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SUMMARY 


This paper presents the results of work on the ejector problem 
carried out by the Research Department, United Aircraft Cor- 
poration, over the past few years. Equations are developed 
which, with the aid of empirical coefficients, enable one to esti- 
mate the performance of steady flow ejectors pumping a fluid 
across a small pressure rise or drop. Experimental data are pre- 
sented which support the theory and provide information on 
some variables not covered by the theoretical treatment. An 
analogy between steady and intermittent flow is presented which 
makes it possible to estimate the performance of ejectors used 
in cooling aircraft engines. A method of applying the results to 
the design of ejectors is suggested. 

The principal uses for an ejector in aircraft engine installations 
appear to be in those cases in which cooling would otherwise 
have to be obtained by the use of a fan or stalled cowl flaps, or 
in which the exhaust thrust cou'd not otherwise be used because 
of flame, noise, or temperature. The results indicate that little 
gain in thrust can be obtained in an installation in which there is 
a pressure drop through the ejector so that pumping is not neces- 


sary. 


INTRODUCTION 


AS EJECTOR IS A DEVICE in which the kinetic energy 
of one fluid is used to pump a second fluid. It 
consists essentially of a duct through which is dis- 
charged a high velocity jet. The secondary fluid is 
drawn from a region of low pressure, mixes with the 
high velocity jet in the duct, and is discharged to a 
region of higher pressure. The ejector may, of course, 
be operated with no pressure rise or a pressure drop 
through it. In this case, the ejector may still increase 
the net thrust of the installation but it will not act as a 
jet pump. 

For airplane installations in which the propeller or 
flight-induced cooling is inadequate (such as in heli- 
copters or airplanes climbing at slow speeds), some 
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form of pump must be used to maintain the required 
cooling airflow through the engine. Cowl flaps may be 
employed as a pumping device. In this case, however, 
the pressure rise is maintained by stalling the cowl 
flaps, which results in an increase in form drag. The 
two alternatives at present are fans and ejectors. Al- 
though an ejector is a far less efficient pump than a fan, 
it need not take power from the engine and it has the 
advantage that the propulsive efficiency of the exhaust 
jet may actually be improved by mixing it with the large 
mass of slow-moving cooling air. 


THEORETICAL TREATMENT 


The object of this analysis is to develop equations 
for the pressure rise of the secondary stream and for 
the momentum from the ejector system in terms of the 
quantity and state of the primary and secondary 
streams and the dimensions of the ejector. The effects 
of mixing and velocity distribution are accounted for 
by empirical coefficients. Incompressible flow is as- 
sumed in order to reduce the number of variables re- 
quired in presenting an overall picture of ejector per- 
formance and to give equations that can be solved ex- 
plicitly for given dimensions. The effects of compres- 
sible flow are discussed qualitatively in a later section. 


This analysis differs from that of Flugel' and Gosline 
and O'Brien? in that the increase in exhaust momentum 
is calculated and the effect on performance of nonuni- 
form velocity distribution at a section is considered. 


When the velocity distribution at any given section 
-of a fluid stream is not uniform, the mean velocities for 
the continuity equation, the momentum equation, and 
the energy equation are all different, because one is an 
arithmetic mean, the second a root mean square, and 
the third a root mean cube velocity. The empirical 
coefficients Cy and Cy are used to relate the mean veloc- 
ities for momentum and velocity head to that for flow, 
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V. Thus, the equations for flow, momentum, and 
velocity head at a given section are: 


Flow = W = pgAV (1) 
Momentum = F = W(VCy)/g (2) 
Velocity head = p(VCy)?/2 (3) 
where 
W = weight flow 
p = mass density 
g = gravitational constant 
V = velocity 
A = cross-sectional area 
F = momentum at a section 
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Fic. 1. Schematic diagram of ejector. 


Fig. 1 is a schematic diagram of an ejector, indicating 
the sections considered in this analysis. The following 
notation is used: 


NOTATION 


Superscripts: 
’ refers to primary jet 
’’ refers to secondary fluid 
Terms without ‘superscripts refer to both streams or the 
mixed stream 


Subscripts: 
A refers to conditions at zero velocity before the ejector 
B refers to conditions before ntixing in the throat of the ejector 
C refers to conditions at end of mixing section 
D refers to conditions at exit of diffuser 


The pressure changes in the entrance section and the 
diffuser sections are found from Bernoulli’s equation 
and the definition of efficiency: 





"(Ve"Cun")? 
P,’ — P; = p (Va Can") (4) 
2nap 
VoCuc)? VpCup)? 
, a eon k ene _ pf D =] (5) 
2 
where 
P = pressure 


4 ideal pressure drop ae 
4p actual pressure drop Comeenny i ee 
trance section for a given change in veloc- 
ity 
actual pressure rise . 








= — -— in the diffuser for a given 
— ideal pressure rise _— 


change in velocity 


1946 


The pressure rise in the mixing section is found from 
the momentum equation and the pressure drop due to 
friction. The friction pressure drop is based arbitrarily 
on the secondary velocity at section B. 
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where the friction factor, f, is defined by the equation 
Friction pressure drop = (fpV?/2)(L/D) 
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In order to reduce the number of variables, the fol- 
lowing dimensionless ratios are introduced: 


length of mixing section wall 
diameter of mixing section 


Area of diffuser exit 





Area of mixing section 
Area of mixing section 
Effective area of primary jet 





_ Secondary flow 
5 Primary flow 
_ Secondary density 
7 Primary density 








For incompressible fluids or for the same gas at different 
temperatures in the primary and secondary streams, the 
density of the mixed stream is found from the following 
relation: 
wi WwW’ + WwW" ge. Bop 

°* (W7e) + W]e") p+ (u/0) 
Introducing the dimensionless ratios, and combining 
Eqs. (1), (4), (5), (6), and (7), the following relationships 
for the pressure rise from the various sections to the dis- 
charge pressure are obtained: 


Pp — Pe _ nevCuc® (1 + #) [1+ aaa F a (<2) 
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e 
The momentum of the combined stream at the ejec- 
tor exit is found from Eq. (2), 


F=W VpCup/g (11) 


The ratio of this momentum to the momentum of the 
primary jet is found from Eqs. 1, 2, 7, and 11: 
1 1 ‘a 
F/F' = (1 + #)( + u/e) Cun (12) 
arC yp’ 

Not allof the coefficients in Eqs.(10)and(12)are worth 
considering; in fact, a fair correlation between test re- 
sults and theory can be found by assuming the coef- 
ficients naz”, Cue’, Cue”, Cuz’, and Cyz” equal to unity, 
neglecting the skin friction, and considering Cyp and 

Cue” Cun \? 
the group Cyc — - nee -() | (hereafter 


- aCe 


called K) to be functions of the diffuser area ratio only. 
Eqs. (10) and (12) then become 
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The experimental flow and thrust coefficients (K) 
and (Cyp) have been calculated from data on com- 
pressed air tests of single nozzle ejectors having a cylin- 
drical mixing section five diameters in length and hav- 
ing 8° (included angle) conical diffusers. Fig. 2 gives 
the average values of these coefficients as a function of 
the diffuser area ratio, a. Fig. 3 shows a comparison 
of the theoretical curves based on the coefficients given 
in Fig. 2, and cross-plots of the actual experimental 
data. These curves give the pressure rise ratio for the 
secondary stream and the thrust ratio as a function of 
the mixing section area ratio, \, for a constant flow ratio, 
nu, of 2. Substantially the same degree of correlation 
was found for other flow ratios. This appears to justify 
the direct use of the theoretical equations in predicting 
the performance of this type of ejector. While better 
values of the empirical thrust and flow coefficients, 
Cup and K, were observed in some tests, these are 
probably as good as could be expected in an actual in- 
stallation, so that their use for preliminary estimates 
seems justified. 


EFFECT OF PRIMARY VARIABLES ON EJECTOR 
PERFORMANCE 


For a given primary jet, ejectors can be designed 
which pump either a little air across a high pressure 
rise or more air across a small pressure rise with an in- 
crease in momentum. The effect of varying the pri- 
mary jet in a given ejector is more difficult to visualize, 
because a change in the primary jet affects several of 
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Fic. 2. Empirical thrust and flow coefficients. 
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Fic. 3. Comparison of data with theory. Flow ratio, 1» = 2; 
density ratio, s = 0.92; a = diffuser area ratio. 


the dimensionless groups of Eqs. (13) and (14) simultane- 
ously. Nor can the effect of changes in the density 
ratio be seen at a glance; its effect on the flow or thrust 
must be found by solving the equations for the particu- 
lar case at hand. 

For the purpose of illustrating the design parameters 
and showing the characteristic performance of steady 
flow ejectors, a series of curves has been calculated 
based on the experimental coefficients given in Fig. 2 
for a density ratio of 3.5, this being representative of 
an ejector using steady flow engine exhaust gas as the 
actuating jet. Since the flow and the pressure rise are 
usually specified in ejector design, plots of pressure rise 
vs. mixing-section area ratio are presented for flow 
ratios of 2, 6, and 10 as shown in Fig. 4. For a typical 
engine, a pressure rise of 10 in. of water corresponds to 
a pressure rise parameter of 0.03. As would be ex- 
pected, the curves show that in order to get more pres- 
sure rise for the same flow and mixing section area, 
more diffuser must be used. The effectiveness of add- 
ing a diffuser falls off as more and more is used. It 
should also be noted from these curves that for any 
given flow ratio and diffuser area ratio, there is a mix- 
ing section area that gives the greatest possible rise. 
If the maximum pressure rise is not required, then a 
range of mixing sections can be used. The use of a 
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Fic. 4. Theoretical pumping capacity. Density ratio, ¢ = 3.5; 
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density ratio, o = 3.5. 


smaller mixing section would give more thrust and re- 
sult in a smaller ejector, but, on the other hand, the 
mixing section area ratio is less critical if the larger mix- 
ing section is used, so that it is easier to meet the de- 
sign point in that region. 

The points of maximum pressure rise for a given flow 
and diffuser area ratio can be found analytically by dif- 
ferentiating the flow equation, Eq. (13), with respect to 
the mixing section area ratio, A, equating to zero, and 
eliminating the flow coefficient, K, by Eq.(13). The 
result is: 
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The second term is negligible. The simplified equation 
determines the locus of the points of maximum pressure 
rise for a given flow and diffuser area ratio, or the points 
of minimum diffuser area for a given flow and pressure 
rise, as shown in Fig. 4. The use of this equation will 
be described in the design section. 

For most aircraft installations the momentum of the 
ejector stream is of interest. Reference to Eq. (14) 
shows that the smaller the diffuser and mixing section 
area ratios, the greater the momentum for a given flow. 
But since these area ratios must be adjusted to main- 
tain the required pressure rise, there is an optimum dif- 
fuser area ratio for a given flow and pressure rise. Fig. 
5 shows the ratio of total to primary thrust vs. diffuser 
area ratio at different pressure rises and flow ratios of 
2, 6, and 10, respectively. It will be noted from these 
curves that for a given flow ratio the diffuser area ratio 
for maximum thrust is greater at the higher pressure 
rises. The thrust ratio is greatest with high flows and 
a pressure drop across the ejector. Fig. 6 is a cross-plot 
of the maximum thrust points for a flow ratio of 6, 
showing the increase in thrust as the pressure rise de- 
creases and finally becomes a pressure drop. The use- 
fulness of an ejector as a thrust-augmenting device can 
only be judged by comparing it with other systems. 
The dotted lines on the right-hand side of Fig. 6 show 
the drag that might be associated with pumping the 
cooling air out by means of cowl flaps or the loss in pro- 
peller thrust that would result from using a cooling fan. 
The exhaust thrust less this drag is the net thrust of the 
alternate system and is less than the total thrust of 
the ejector, shown by the solid line, so that in this re- 
gion the ejector will give better performance than other 
systems. When there is a pressure drop through the 
ejector, however, the net thrust of the alternate system 
is that of the exhaust thrust plus the momentum that 
the cooling air would have if discharged separately. 
This momentum, neglecting losses, would be: 


F” = W"V,"/g (16) 
The velocity may be found from the available pressure 


drop of the secondary stream; hence, the momentum 
is: 





F" = (W"/g) V2(P” — Pp)/p" (17) 


Expressed in terms of the properties of the primary jet, 
this momentum is: 

F” 5 

o 


Feo (ee Fo (18) 


FTA — 
Fig. 6 shows that with an increasing pressure drop 
through the ejector, the sum of the exhaust thrust and 
this secondary momentum soon equals the ejector 
thrust. 
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For a typical engine exhaust and a flow ratio of 6, 
the ejector yields no gain in thrust when the pressure 
drop through the ejector is greater than about 5 in. 
of water. This characteristic of ejectors limits their 
usefulness as thrust augmenters to a small range. 
They are nonetheless useful in reducing the drag as- 
sociated with pumping the cooling air, if that is neces- 
sary. 


EFFECT OF COMPRESSIBILITY 


The experimental coefficients in the steady flow 
theory were derived from compressed air tests in which 
the primary jet velocity was approximately 700 or 800 
ft. per sec. Steady flow tests with superheated steam 
as the primary fluid were, however, run with primary 
velocities from 1,200 to 2,800 ft. per sec. The latter 
tests showed that the design of an ejector is not affected 
by the primary jet velocity. For supersonic velocities 
of the primary jet the performance of the ejector was 
not appreciably affected by removing the divergent 
section of the nozzle. In designing an ejector, how- 
ever, the area of the fully expanded stream, rather than 
the area of the nozzle, must be used. The temperature 
of the primary jet should also be calculated assuming 
complete expansion of the jet to the pressure in the ejec- 
tor throat. 

The high-velocity steam tests showed that the per- 
formance of an ejector fell off sharply when the flow 
mean velocity of the mixed streams reached a Mach 
Number around 0.6. After this velocity was reached, 
practically no increase in secondary airflow or total 
ejector thrust could be obtained by further increasing 
the primary flow and primary velocity. It is, therefore, 
necessary to check the velocity in the mixing section 
when designing an ejector for a given condition. If the 
calculated velocity yields a greater Mach Number than 
0.6, the ejector should be redesigned with a larger mixing 
section area. The theoretical curves presented in this 
report were based on experimental data in which the 
Mach Number of the mixed streams was not over 0.3. 

The effect of compressibility was also checked ana- 
lytically by assuming a primary jet Mach Number of 
unity, zero pressure rise, and a mixing section area ratio 
of 20. The thrust for a given flow ratio was about 2 
per cent less for compressible flow, and the flow ratio 
for a given diffuser area ratio was about 5 per cent 
greater for compressible flow. The agreement was 
even better for ejectors without diffusers, so that the 
assumption of incompressible flow seems reasonable. 


EFFECT OF ALTITUDE 


Data obtained from sea-level tests of ejectors can- 
not be applied in designing an ejector to cool an engine 
at altitude without considering the effect of the change 
in densities. A decrease in the density ratio by itself 
results in a loss in ejector performance. For a given 


primary flow, however, the jet thrust increases with al- 
titude because of the reduced back pressure on the noz- 
zle. The increase in jet thrust has a greater effect upon 
the performance than has the change in density ratio. 
Therefore, if the velocity of the steady flow primary 
jet is subsonic, the overall performance will improve 
with an increase in altitude. 

If this jet velocity is sonic at any given altitude, a 
further increase in altitude may be accompanied by a 
decrease in performance in spite of the increased jet 
thrust. The reason for this is that the jet expands to a 
greater area after leaving the nozzle, giving in effect a 
smaller mixing section area ratio. 

The quantitative effect of altitude on the perform- 
ance of an ejector will depend upon the design condi- 
tions and the particular ejector chosen for these condi- 
tions. Fig. 7 shows the calculated effect of altitude on 
the performance of a steady flow ejector designed for 
an arbitrary set of conditions. The variation in flow 
and thrust ratio with altitude is calculated for a con- 
stant actual pressure rise. The flow and thrust ratio 
in this case increase with altitude until the jet velocity 
becomes sonic. The thrust -remains constant and the 
flow decreases with further increase in altitude. 
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Calculated effect of altitude. Constant pressure rise. 
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Fic. 7. 


DIVERGENT MIXING SECTIONS 


The ejector with a straight mixing section is not only 
the easiest to calculate but it also appears to be theoret- 
ically the best in the absence of losses. No cases have 
been observed, either in theory or in experiment, in 
which a divergent or convergent mixing section led to 
more thrust or a smaller ejector than a straight mixing 
section. One exception might be in the case of high 
Mach Numbers of the secondary stream. In any case, 
once mixing is complete, the addition of a diffuser sec- 
tion has advantages that have already been discussed. 


INTERMITTENT FLOW EJECTORS 


If the exhaust from all cylinders of an engine were 
collected in a manifold and then discharged through a 
nozzle, the flow would be steady enough to use the 
steady flow theory. However, this system would not 
utilize the impulses of kinetic energy from each cylinder, 








564 JOURNAL OF THE AERONAUTICAL 


so that other arrangements are generally used which 
involve intermittent flow, and ejectors must be designed 
to operate with intermittent flow jets. 

An analogy can be made between steady flow and in- 
termittent flow by assuming the flow is steady in the 
latter case and using the primary area that would give 
the same average primary flow and momentum as the 
engine exhaust. The effective primary area can then 
be found from the continuity equation: 


A’ = W'/p'gV’ (19) 


and the mean velocity calculated from the primary 
thrust, which is measured by a balance or a thrust 
target: 


F’ = W’'V'/g 20) 
Combining these, 

A’ = (W’)?/p'g?F’ (21) 
The principal objection to making this analogy is 
that it disregards the effect of the geometrical rela- 
tionship between the area of the ejector and the primary 
nozzle. During the period of the cycle when the pri- 
mary nozzle is actually exhausting the area ratio is 
smaller than that assumed in using the calculated steady 
flow primary area. The actual pumping capacity is 

consequently poorer than predicted by the theory. 





ae 


f ee ee ae 
a: »* | MEASURED | 
+ 20: \ t +—_$AA 4-1 































































































20) =a ’ 3 | 
| | | | | | | \ | j 
= a ms ! ie) 
% je | | “CALCULATED; | NIN 
~'5 ss 7 ee uw fin 15; i 
= wird | i KIS, 
& PRES Te Ce 
=e [ZN ~}] iS 10 wy. | iN a 
> MEASURED ' CALCULATED | \ | [ 
Ss x! FROM MEASURED 
at 2 5, FLOW \ 
” = els x 
' 
i 
| | 
r) 


te) : 
“10 “5 ie) 5 4 “5 1e) 5 
PRESSURE RISE, P,.—R-, INCHES OF WATER 


Fic. 8. Comparison of theory with single-cylinder engine ejector 
data. 


SINGLE-CYLINDER EJECTORS 


Fig. 8 shows the experimental results of typical single- 
cyclinder ejector tests compared to the calculated per- 
formance using the flow and thrust coefficients, A and 
Cup, derived from steady flow compressed air tests. 
The results compare favorably only when there is a 
fairly high pressure drop through the ejector. The dis- 
crepancy is greatest when the ejector is pumping against 
a pressure rise and the return flow of the secondary 
fluid is decreasing the quantity of fluid pumped. Also, 
the nonuniform distribution in time of the velocity of 
the secondary fluid results in mean thrust per unit flow 
greater than would be obtained in steady flow of the 
secondary fluid. This variation between theory and 
test was consistent for all ejectors without diffusers. 
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Fic. 9. Correlation of flow data—single-cylinder 
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The variation was more erratic for ejectors with dif- 
fusers, but the difference between theory and test in no 
case appears to be greater for ejectors with diffusers 
than for ejectors without diffusers. It is evident that 
the empirical coefficients determined for steady flow 
ejectors cannot be used to predict the performance of 
intermittent flow ejectors, so that other means for 
handling the data must be found. 

Figs. 9 and 10 show a correlation of data obtained 
from tests of ejectors without diffusers actuated by a 
single cylinder. These plots seem to correlate fairly 
well data taken over a wide range of engine powers, 
speeds, fuel-air ratios, and exhaust stack restriction. 
Fig. 9 gives the pressure rise parameter as a function 
of the calculated area ratio for various flow ratios. 
The density ratio is not constant for these tests (vary- 
ing from 3.5 to 4.5), which accounts for some of the 
scatter in the data. 

Fig. 10 is a plot of the ratio of the total to primary 
momentum vs. the effective area ratio for various values 
of the flow ratio. The scatter due to variation in den- 
sity ratio may be as much as 7 per cent for these data, 
but there are some isolated points that fall as much as 
20 per cent off the average curve. These variations are 
presumably due to changes in the time distribution of 
the secondary velocity, which was not accounted for in 
any way by the analogy. 

These plots should not be used in designing ejectors 
in which ‘‘Siamesed stacks’ are employed. An ejec- 
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tor actuated by the exhaust from a number of cylinders 
having nonoverlapping exhaust periods will pump a 
considerably greater flow for a given mixing section 
area ratio than indicated in Fig. 9. The ejector thrust 
in this case cannot be compared with that given in Fig. 
10 because the mixing area required to maintain a 
given flow is not the same as that for a single-cylinder 
ejector. 

There are insufficient data to plot curves similar to 
those in Figs. 9 and 10 for ejectors with diffusers. 
Moreover, these curves do not show the effect of 
changes in density ratio and cannot, therefore, be used 
in estimating performance at altitude. 

Even if the theory were developed to completely ac- 
count for all the variables in intermittent flow, it would 
be of little value at present since the additional data 
which would be required are not available. It is, there- 
fore, absolutely necessary at present to employ an 
analogy such as given above in spite of its shortcom- 
ings. 


EJECTORS WITH SIAMESED STACKS 


In most engine installations it will be found advan- 
tageous to discharge the exhaust from several cylinders 
into one ejector. The greatest thrust is obtained by 
having a separate stack for each cylinder, but the piping 
can be reduced at a small loss in thrust by Siamesing 
into a single outlet two or three stacks from cylinders 
whose exhaust periods do not overlap. The outlet 
area should be the same as that for one individual stack. 
In this case the primary flow will be more continuous 
and should approximate more closely that calculated 
from steady flow theory. However, there are available 
no experimental data from engine ejectors with Sia- 
mesed stacks which are complete enough to check the 
theory. 

Figs. 11 and 12 show a comparison of the theory with 
results of intermittent flow compressed air tests. Three 
intermittent flow nozzles were used. The nozzles dis- 
charged for one-third of the cycle with a velocity varia- 
tion representative of that in the exhaust flow from an 
engine cyclinder. In the first case, three separate 
ejectors were used so that intermittent flow was ob- 
tained through each. The pumping performance of 
these ejectors was markedly lower and the thrust higher 
than that given by theory. This is in agreement with 
the comparison made between the theory and the re- 
sults of single-cylinder engine tests. In the second case, 
the nozzles were Siamesed into a single nozzle and 
used in a single ejector having an area equal to the com- 
bined area of the three intermittent flow ejectors. The 
performance in this case was in good agreement with 
the theoretical calculated performance. The use of 
three individual stacks instead of the Siamesed stacks 
in a single ejector gave the same performance in the 
compressed air tests, but slightly better performance 
when tried on a six-cylinder Franklin engine. 
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These intermittent flow compressed air tests provide 
a good enough qualitative check to warrant the use of 
the theory when Siamesed stacks or a number of indi- 
vidual stacks from nonoverlapping cylinders are used 
in an ejector so that reasonably continuous flow is 
maintained. Even if this condition is fulfilled, the re- 
sults calculated using the empirical thrust and flow 
coefficients given in Fig. 2 should be considered some- 
what optimistic with respect to the flow and conserva- 
tive with respect to the thrust. Coefficients based on 
engine tests of the type of ejector being designed should 
be used whenever the data become available. 


RESTRICTED EXHAUSTS 


The exhaust momentum in any type of system can 
be increased by restricting the exit area of the exhaust 
pipe. The amount of restriction that can be made 
without a decrease in engine power depends on the type 
of exhaust system and the exhaust port design of the 
engine. When individual stacks are used, a loss in en- 
gine power will result if the stack area is made less 
than the area of the exhaust valve passage at the maxi- 
mum opening of the valve. A restriction in the exhaust 











566 JOURNAL OF THE AERONAUTICAL SCIENCES—NOVEMBER, 1946 


stacks may, however, result in an overall gain in thrust 
at high speed, if the increase in jet thrust more than off- 
‘sets the decrease in propeller thrust due to loss of engine 


power. 


MIxINnG LENGTH 


As the mixing length of an ejector is increased from 
zero, the velocity distribution becomes more uniform 
and the performance of an ejector improves at first. 
On further increase, however, wall friction in the mixing 
duct begins to become important and cuts down the 
performance. The optimum length results from a com- 
promise between these two factors, and has been found 
by various investigators to be from 4 to 8 L/D for 
single-jet ejectors. Fig. 13, from tests on a seven- 
nozzle ejector, shows the interaction of these two fac- 
tors and indicates that there is a critical mixing length 
below which the performance falls off abruptly. 
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For ejectors with several jets, the critical mixing 
length can be approximated by assuming that the re- 
gion of mixed stream forms a cone of fixed angle whose 
apex is at the primary jet. This leads to the rule that 
the critical mixing section length will vary directly with 
the mean radial spread of .the mixed region to the wall 
or to another jet, provided that the nozzles are distrib- 
uted uniformly over the mixing section area. Fig. 
14 shows results of compressed air tests which support 
this rule for various shapes of mixing sections as well 
as various numbers of nozzles. The diameters of non- 
circular shapes are taken to be the diameters of circles 
with equal area. The mean radial spread of the mixed 
region for multiple-nozzle ejectors is assumed to be the 
same as that for a circular single-nozzle ejector whose 
area is equal to that of the whole ejector divided by the 
number of nozzles. Preliminary tests showed that 
annular nozzles will also shorten the mixing section if 
secondary air is allowed to flow through the annulus. 
All these data were obtained with zero pressure rise 
across the ejector. It should be emphasized, however, 
that no advantage will be obtained from distributing 
the nozzles in the ejector area unless there is essentially 
a continuous flow from each nozzle. The pressure rise 
also affects the optimum mixing length. Fig. 15 shows 
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Fic. 14. Critical length of mixing section. 








that for the same primary jet, the pressure rise in the 
mixing section takes place in a much shorter distance 
when less rise is required, as indicated by the leveling 
of the curves. Fig. 13 shows that the critical mixing 
length actually is less with a diffuser, for which less pres- 
sure rise occurs in the mixing section even with the 
same pressure change across the entire ejector. 

The performance of an ejector without a diffuser is 
not appreciably affected by breaking up the primary 
jet into several scattered jets, provided that the mixing 
section is above the critical length. The improvement 
in performance due to the smaller friction loss of the 
shorter mixing section used is relatively small. How- 
ever, the improved velocity distribution that results 
does increase the efficiency of the diffuser, thereby 
considerably improving the performance of ejectors 
with diffusers. Fig. 16 shows the comparison between 
single- and multiple-nozzle ejectors with diffusers of 
varying area ratio. 


ENTRANCE CONFIGURATION 


Tests have indicated that the entrance radius of the 
mixing section can be as low as one-sixth of the diameter 
without affecting the performance. The nozzle should 
be flush with the entrance or slightly outside the en- 


trance. With steady flow nozzles no appreciable loss 
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in performance will be incurred by having the nozzle 
extend into the duct, other than that which may result 
from reduction of the effective mixing length. Since 
the distance from the exhaust stacks to the ejector exit 
is usually critical in ejector installations, there is no 
point in withdrawing the nozzles far out of the duct. 

In general, the same considerations hold for intermit- 
tent flow ejectors as for steady flow ejectors. It is 
important that the exhaust be correctly distributed in 
the ejector area. Care should be taken that the stacks 
exhaust straight in the ducts and do not strike the duct 
walls. It is extremely important that the stacks do 
not extend into the duct unless three nonoverlapping 
cylinder exhausts are Siamesed together for each nozzle. 
Otherwise, during any given period of the cycle the mix- 
ing area is reduced by the blocking effect of the stacks 
that are not firing. 


DIFFUSER ANGLE 


The experimental coefficients used in the theoretical 
calculations were based on ejectors having 8° conical 
diffusers. The efficiency of such a diffuser is not con- 
stant for different diameters and different area ratios. 
However, in lieu of a more thorough analysis of dif- 
fusers and adequate experimental data, it is suggested 
that 8° conical diffusers be used or, in the case of non- 
circular cross sections, that diffusers with equivalent 
cone angle be used (i.e., maintain the same change in 
area for a given length that a cone with an included 
angle of 8° would have). Experience has shown this to 
be a fairly good criterion for designing a diffuser, al- 
though it is not fundamentally correct. 

For a given overall length somewhat more flow can be 
obtained at the expense of a loss in thrust by using a 
greater diffuser angle. This would be practical only in 
a case where the overall length was limited and it was 
desired to obtain the maximum flow without regard to 
the thrust. 

It has also been found that the optimum diffuser 
angle is greater than 8° for extremely small diffuser 
area ratios. The data are meager, however, and no 
reliable recommendations can be made for such cases. 


DESIGN OF AN INSTALLATION 


The foregoing equations and figures should enable 
the designer to make a rough check on the suitability of 
an ejector for his particular needs. The final design 
will, for the time being, have to be verified by flight 
tests, preferably of several configurations, which seem 
to be near the optimum. 

In any type of installation it is first necessary to find 
the state of the secondary fluid, the required pressure 
rise, and the mean state of the primary fluid as described 
by the intermittent flow analogy (see Eqs. (19), (20), and 
(21). Strictly speaking, the stateof the primary fluid 
should be obtained at the pressure of the ejector throat, 
but for most airplane installations the state at atmos- 
pheric pressure is close enough. It should be remem- 
bered, however, that in the case of a supersonic primary 
jet, the primary area should be calculated as if the jet 
had expanded through a converging diverging nozzle, 
rather than using the throat area as the primary area. 
The pressure rise of the secondary fluid is usually taken 
to atmospheric pressure, mainly because of lack of bet- 
ter knowledge of the static pressure at the ejector exit. 
The calculated thrust will then be roughly correct, 
even though the calculations may indicate a larger dif- 
fuser than will actually be needed. Having determined 
the states of the primary and secondary fluids, the de- 
signer can calculate the required flow ratio, density 
ratio, and pressure rise coefficient. 


Since an adjustable diffuser is mechanically difficult 
to build, this will be one of the limitations in design. 
The overall length and the ejector area will be two other 
limitations imposed by the airplane designer and these 
must be balanced against the desired performance. If 
the thrust is not important but the diffuser area ratio 
is, then the mixing section area for the minimum diffuser 
can be found from Eq. 15. The required diffuser area 
ratio can be estimated by solving the flow equation 
(Eq. 13) for K and finding the area ratio from the empir- 
ical plot of flow coefficient, A, vs. diffuser area ratio, 
Fig. 2. If length is a more critical item than diffuser 
area ratio, a reduction may sometimes be made by using 
a slightly smaller mixing section area and a larger dif- 
fuser area ratio. The mixing length may be reduced 
by the use of multiple or annular nozzles. 


If thrust is most important, the diffuser area ratio 
for maximum thrust may be estimated from Fig. 5 
and the mixing section area found from Fig. 4 at the 
given pressure rise coefficient, flow ratio, and diffuser 
area ratio. If the density ratio differs much from 3.5, 
or if more precision is desired, the optimum diffuser 
must be found by assuming several values of mixing 
section area ratio, solving the flow equation (Eq. 13) 
for the flow coefficient, A, finding the corresponding dif- 
fuser area ratio, a, and then solving the thrust equation 
(Eq. 14) to find which assumed area ratio gave maxi- 
mum thrust. For more accurate work, the state of the 
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primary fluid should be taken at the pressure in the 
ejector throat rather than at atmospheric pressure. 

In the usual airplane installation, however, if an 
ejector is required at all it should be designed to oper- 
ate with the maximum allowable diffuser in the most 
critical cooling condition. Then the diffuser can be 
closed down at high speeds to prevent overcooling and 
to increase the thrust. This procedure will seldom re- 
sult in a diffuser at high speed that is greater than the 
optimum for thrust. 


OUTSTANDING PROBLEMS IN EJECTOR DESIGN 


There are several uncertainties regarding the per- 
formance of ejectors which have not been adequately 
dealt with. First, there are practically no data on the 
effect of ejectors on low pressure areas into which they 
may be discharged, so that for the time being the as- 
sumption of atmospheric discharge pressure is used. 
In the second place, the mixing process has not been 
considered in detail except by the limited work of 
Ledgett.* Certainly the optimum mixing length varies 
with Reynold’s Number, the steadiness of the flow, the 
pressure rise on the mixing section, and other quantities. 
It is still necessary, however, because of lack of better 
knowledge, to assume that the mixing cone spreads out 
at a constant angle. A third uncertainty may be ob- 


served by studying the plot of flow vs. pressure rise in 
Fig. 11. Note that the points are actually waving back 
and forth across the mean line. This waviness is not 
experimental scatter but was observed to some extent 
at high and low Mach Numbers, with steady and inter- 
mittent flow ejectors, and with both air-air and steam- 
air ejectors. A fourth problem is that of intermittent 
flow. Although the data usually correlate, as demon- 
strated in Figs. 9 and 10, occasionally some peculiar 
resonance effects will appear, upsetting all estimates of 
performance. 

These problems, along with lack of accurate and 
complete knowledge of engine exhaust thrust, make it 
impossible to design ejectors precisely. At the present 
time, therefore, it seems sufficient to estimate the 
ejector performance roughly from the foregoing plots 
and equations, and then to modify the final design after 
flight-test data are available. 
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Determination and Presentation of 
Compressive Stress-Strain Data 


for Thin Sheet Metal 


WALTER RAMBERG* anv JAMES A. MILLER? 
National Bureau of Standards 


ABSTRACT 


Various methods of obtaining the compressive stress-strain 
curve of thin sheet metal are described, and their advantages and 
limitations are briefly discussed. The presentation of stress- 
strain data in dimensionless form is outlined, and the advantages 
of such a presentation in correcting for the effect of variations in 
yield strength and in Young’s modulus are brought out. Atten- 
tion is called to a classification of stress-strain curves by shape 
factors following a three-parameter formula proposed by Osgood. 


INTRODUCTION 


BS Iyesn CURRENT INTEREST in methods of obtaining and 
presenting stress-strain curves of thin sheet metal 
in compression dates from the introduction of stressed 
skin structures in aircraft 10 or 15 years ago. The first 
stressed-skin structures were relatively inefficient be- 
cause buckling in compression could be expected well 
inside the elastic range and, therefore, there was little 
need to know compressive properties beyond the elas- 
tic range. As the structures became more efficient, it 
became of increasing importance to know not only 
Young’s modulus, which describes the stress-strain 
curve inside the elastic range, but also the stress-strain 
curve beyond the elastic range. In the beginning the 
materials testing engineers were able to pacify the 
designers who were asking for this information by mak- 
ing tensile tests preferably in a transverse direction and 
telling them that the stress-strain curve in longitudinal 
compression would be approximated by the tensile 
stress-strain curve in a transverse direction. This was 
good enough for some commonly used aluminum al- 
loys, such as 17S-T and 24S-T,! but it soon proved en- 
tirely inadequate for other alloys, particularly stainless 
steels' and magnesium alloys.? In the case of these 
alloys no simple and reliable method has been found for 
determining longitudinal compressive properties from 
any tensile properties. In the end the materials engi- 
neers had to give in and make compressive tests of these 
alloys to determine their compressive properties. 
Several methods of testing have been developed and 
used. All of them have their limitations and none of 
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them is in sufficiently general use to have earned for it- 
self inclusion in the A.S.T.M. Standards. 

A brief survey of the various known methods seems 
advisable at this time for the benefit of structural test 
engineers who are faced with the problem of making 
compressive tests without having had the time to look 
into the various methods. 

The National Bureau of Standards is favorably situ- 
ated for assembling such a survey, since it has been 
concerned with developing methods of compressive 
testing of thin sheet metal for the Bureau of Aeronau- 
tics, Navy Department, and for the National Advisory 
Committee for Aeronautics since 1933. The N.B.S. has 
also been concerned with methods of presenting these 
data in a'rational form which would reduce the test- 
ing of specimens of similar materials to a minimum. 

In reviewing the various methods for making com- 
pressive tests of thin sheet metal one must admit from 
the outset that the experience of the National Bureau of 
Standards is limited insofar as it does not cover all of . 
the ingenious methods that have been developed. In 
several cases, therefore, the authors had to base their 
review on articles in the technical literature rather 
than on their own experience. 


METHODS OF DETERMINING THE COMPRESSIVE STRESS- 
STRAIN CURVES OF THIN SHEET METAL 


All methods of making compressive tests of thin 
sheet metal have one problem in common—namely, the 
problem of overcoming the tendency of the specimen to 
become increasingly misaligned with increasing load 
and finally to fail by instability. An ideal solution of 
the problem requires the provision of simple and easily 
applied supports that prevent butkling without intro- 
ducing friction between supports and specimen and 
without introducing restraint to lateral expansion of 
the specimen as the axial load is increased. 

It will be noted in this connection that the tensile 
test is favored relative to the compressive test insofar 
as alignment and stability improve with increasing load 
instead of becoming worse. The testing machine may 
add to the difficulties of the compressive test by becom- 
ing unstable’ or otherwise permitting a lateral displace- 
ment or a rotation of one head relative to the other. 
In order to obtain suitable loading conditions, various 
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Pack Method 0 
The earliest successful method for testing thin sheet T 
in compression of which the authors have any knowledge | el 
is the pack method developed by Aitchison and Tucker- al 


man of the National Bureau of Standards.’ In the 
pack method the leaf to which the gages are attached 
is supported on each side by an equal number of leaves 
from the same sheet which are identical with the middle 
leaf except for being slightly narrower to ensure that 
the strain gages are attached only to the middle leaf. 
The ends of the ‘“‘pack’’ consisting of the middle leaf 
and the supporting leaves are ground flat and parallel, 
and the pack is then tested as a single stocky compres- 
sive specimen (Fig. 3). Lateral buckling of the outer 
leaves of the pack is prevented by three rows of steel 
pins on each side, '/; in. between centers, which are sup- hardened and ground steel bearing blocks with an in- 
ported by a rigid frame connected to the lower platen termediate shim of plaster of Paris which is allowed to 
of the testing machine. The pack is loaded through harden under an initial load of about 500 Ibs. Strains 
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Fic. 3. Pack ready for test. 
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are measured by a pair of strain gages attached to the 
middle leaf of the pack. 

A number of modifications of the pack method have 
been made to simplify it and to extend its use to thinner 
sheet or stronger materials. Among these are the use 
of clamps, which are placed near the ends of the pack 
before machining the ends and which are left on during 
the test, and the use of a spherically seated bearing 
block at one end of the specimen. Various sizes of spec- 
imens have been used. In general, the length is at 
least twice the width plus the gage length. The Alu- 
minum Research Laboratories have found it unneces- 
sary to make the middle leaf of the pack wider than the 
supporting leaves. Aitchison® found that thinner ma- 
terial could be tested when the leaves were bonded to- 
gether by shellac applied in liquid solution at room 
temperature, dried, and then baked in a clamp at a 
temperature of about 65°C. Using this last technique 
Aitchison was able to determine compressive stress- 
strain curves beyond the yield strength on packs of 
0.020-in. sheet material with a yield strength of more 
than 215,000 Ibs. per sq.in. 

Comparative tests showed that the yield strength of 
a material determined by the pack method was very 
nearly the same as that determined by a test on a solid 
specimen. 


Single Strip Supported by Rollers 


Several schemes were soon tried to improve on the 
pack test by reducing the time required for machining 
the specimen, the sample of material required, and the 
maximum load required of the testing machine. One 
of the most successful of these tests is the test of a single 
leaf or strip supported by rollers, which was suggested 


by W. P. Montgomery and developed further by R. L. - 


Templin and his associates of the Aluminum Research 
Laboratories.’ 

Thé fixture is shown in Fig. 4. The specimen con- 
sists of a piece °/; by 2°/s in. cut from the sheet metal. 
It is slipped between two sets of hardened steel rollers 
0.093 in. in diameter and spaced 0.1 in. between centers. 
The rollers are clamped against the specimen by tight- 
ening the screws. A portion of the specimen, usually 
about 0.03 in. long, projects above the fixture to provide 





Fig. 4. Specimen in Montgomery-Templin fixture. 


clearance between the upper bearing block and the fix- 
ture at the maximum strain. Strain is usually meas- 
ured by a pair of gages attached to the edges of the spec- 
imen. Sometimes an autographic gage is attached to 
one of the edges of the specimen. 

Successful compressive tests, with results in close 
agreement with those obtained by the pack method and 
by tests on solid specimens, have been made with the 
Montgomery-Templin fixture on specimens of 24S-T . 
aluminum alloy down to 0.020 in. in thickness.’ For 
thicknesses below 0.020 in., difficulties have been ex- 
perienced because of buckling at the top. This some- 
times causes damage to the roller supports. The diffi- 
culty has sometimes been overcome by cementing to- 
gether a number of leaves as in the “modified pack 
method.’ Kotanchik, Woods, and Weinberger® state 
that the roller supports must be constructed with pre- 
cision, and they report that unreliable test results were 
obtained because of an inaccurately constructed sup- 
port. Indentation of specimens by the rollers has also 
been observed. Modifications in the method include 
the use of a specimen !/2 by 2!/2 in., a cylindrical seat 
for one of the steel blocks, and a single screw working 
against an equalizing bar for adjustment. 

The time saved in preparation of the specimen as 
compared to the pack test is considerable. The saving 
is greater the thinner the material, since a greater num- 
ber of single strips can be machined at one time. Some 
time can also be saved in setting up the specimen for 
test. 


N.B.S. Solid-Guide Method 


The work with solid guides at the National Bureau of © 
Standards was started by C. S. Aitchison in 1939; in 
1942 he substituted lubricated brass plates for the roll- 
ers of a Montgomery-Templin fixture.° Further de- 
velopment of the method resulted in the design of the 
fixture with tool steel guides, Fig. 5, which is described 
briefly in reference 10 and more fully in reference 11. 

The lubricated faces of the steel guides have a length 
of 2.20 in. and a width of 0.46 in. The specimen is 
made about 0.04 in. wider, but in no case is it more than 





Fic. 5. Specimen between N.B.S. lubricated steel guides. 
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four thicknesses wider than the guides. The initial 
unsupported length is made sufficiently small to pre- 
vent buckling at the ends before the yield strength is 
reached; it is made sufficiently large to provide clear- 
ance between the upper bearing block and the guides 
at the maximum strain. The initial unsupported length 
for specimens thicker than 0.02 in. is about 0.05 in.; 
for specimens from 0.01 to 0.02 in. it is about 0.03 in. 
A thin film of heavy-duty lubricant for automobile 
wheel bearings* is spread on the specimen and the 
specimen is placed between the two guides. The adjust- 
ing screws 7 are tightened and retightened with inter- 
vening sliding of the specimen until the guides will not 
slide down under their own weight but willslidesmoothly 
under an additional weight of 1 Ib. or less. After a lit- 
tle practice it is possible to tell by feel when the desired 
condition is obtained. 

The four holding screws S which transmit the tight 
ening force to the guides are of sufficient length and suf- 
ficiently small section to prevent the setting up of ex- 
cessive pressures as the specimen expands laterally at 
high compressive strain. Strains are measured with a 
pair of 1-in. Tuckerman optical strain gages mounted 
on the projecting edges of the specimen. 

The solid-guide method has enabled us to test sheet 
thinner than 0.006 in., with a 0.2 per cent offset yield 
strength of about 160,000 Ibs. per sq.in., and sheet 0.020 
in. thick with a yield strength of about 240,000 Ibs. per 
sq.in. 


* Texaco Marfak No. 3 was the only grease of this type tried. 
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Fic. 6. Compressive data on 0.032-in. 24S-RT aluminum-alloy sheet in direction of rolling, tested by five methods. 


Tests on specimens with thicknesses of from 0.01 to 
0.128 in. have been made by inexperienced personnel 
following routine procedure, with results that closely 
agreed with those obtained by experienced personnel. 
For thinner specimens more care is required, particu- 
larly with regard to the thickness of the lubricant film. 


Tests made with solid guides, with the Montgomery- 
Templin fixture, and by the pack method on various 
sheet metals were reported in reference 11. The re- 
sults obtained by the various methods were in close 
agreement. For example, Fig. 6 shows compressive 
stress-strain data for longitudinal specimens from a 
sheet of 0.032-in. 24S-RT aluminum alloy which were 
obtained by five methods at the National Bureau of 
Standards. The tests were made by the pack method, 
with the Montgomery-Templin fixture, with lubricated 
brass guides, with lubricated steel guides, all with the 
subpress, and with lubricated steel guides without the 
subpress. The differences between the stress-strain 
data are magnified by plotting as a stress-deviation 
curve, the abscissa being the deviation from an average 
straight line through the points in the elastic range for 
the pack and Montgomery-Templin tests. The stress- 
strain data are in close agreement. There is no indi- 
cation of systematic differences between the results 
obtained by the various methods. 


The results obtained with the solid-guide method were 
found to be as consistent as those obtained in tensile 
tests. The time required for the two kinds of test is 
about the same. 
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Consolidated Vultee Solid-Guide Method 


One of the earliest published methods for obtaining 
the compressive stress-strain curve of thin sheet metal 
from tests of single specimens between lubricated solid 
guides is due to K. R. Jackman of Consolidated Vultee 
Aircraft Corporation." 

In this method the specimen consists of a piece 1 by 
3!/g in. with two tabs on one side (Fig. 2), punched 
from the sheet. The ends are finished by milling or 
grinding. The tabs serve as points of attachment for a 
2-in. extensometer, which may be of any type with a 2- 
in. gage length. The specimen is guided between 
greased plates of either metal or wood. Reference 12 
contains no data on the minimum thickness of speci- 
men that can be tested successfully or on the accuracy 
of the test. However, it is emphasized that the method 
has an advantage for industrial use in saving time. It 
is stated that an experienced testing-machine oper- 
ator can carry through a test on a specimen in 
from 8 to 10 min. if he uses a Southwark-Templin 
autographic strain recorder to record the stress-strain 
curve. 

The Consolidated Vultee method has not been used 
at the National Bureau of Standards where precision is 
more important than the saving of time, because errors 
may be introduced by reading strains on tabs on one 
edge of the specimen and errors may be caused by 
changing the properties of the material along the edges 
in the punching operation. 


Welter’s Solid-Guide Method 


Professor Georges Welter of the Ecole Polytechnique 
de Montreal used lubricated guides to support com- 
pressive specimens of 0.032-in. 17S-T aluminum-alloy 
sheet.'* Few details of the method are given. Three 
different devices are shown with cutouts to permit the 
attachment of gages to the sides of the specimen. The 
compressive load is introduced into the specimen 
through blocks that are fixed to the specimen either 
mechanically or with Wood’s metal (Fig. 2). The 
specimens are */, in. wide at the ends and !/2 in. wide 
in the reduced section. 

No data obtained by other methods on the same 
sheet are given from which the reader may judge the 
accuracy of the method. However, autographic curves 
obtained with a nonaveraging 2-in. extensometer are 
not nearly so smooth as those shown for the tensile 
tests. 

An adequate study of the method seems necessary to 
ascertain that the results are not in error because of 
nonuniformity in strain distribution, frictional forces 
between guides and specimen, and buckling in the 
unsupported portions of the specimen. It would be 
desirable to make additional tests to determine the 
applicability of the method to sheet stronger and also 
thinner than that tested. 





Fic. 7. ARMCO single-strip compression jig. 


ARMCO Solid-Guide Method 


LaTour and Wolford of the American Rolling Mill 
Company Research Laboratories have developed a par- 
ticularly compact fixture for testing single strips be- 
tween lubricated steel guides. The procedure is de- 
scribed in detail in reference 14. Fig. 7 shows the spec- 
imen F clamped in the fixture between the movable 
guide B and the flat guiding surface on the rigid frame 
A. A light coat of graphite grease or oil is applied to 
the specimen as lubricant. The guides are clamped 
on the specimen by tightening thumb screws C with 
A positioning ring, D, is used to 
bearing 


moderate pressure. 
center the fixture midway between the 
blocks. 


Strain gages with a l-in. gage length are attached to 
both edges of the specimen, where the guides are cut 
away, to measure strain under load. An approximately 
uniform load distribution is obtained by loading the 
specimen in a testing machine between bearing blocks 
resting on wet blotting paper. A preload of 25 to 100 
Ibs. is applied and the thickness of the blotting paper is 
adjusted by rubbing away at the appropriate side until 
the strain gages indicate approximately the same strain 
on the two edges of the specimen. The fixture is also 
used in conjunction with an autographic compressometer 
that measures the change in distance between the 
bearing blocks. A correction for ‘‘foundation effect’’ 
is given. 


Compressive tests were made with the fixture on 
specimens 0.035, 0.054, and 0.100 in. thick, cut from a 
bar of SAE 4330 steel. The values of initial modulus 
of elasticity and of yield strength from data obtained 
with gages at the edges of the specimen were all within 
2 per cent of the values obtained on a solid cylinder of 
the material, with the exception of the modulus of the 
0.100-in. specimen. As is to be expected, the data ob- 
tained with the autographic recorder were not so con- 
sistent. The yield strength, offset = 0.2 per cent, has 
been obtained by this method on sheet as thin as 0.01 
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in. with a yield strength of over 140,000 Ibs. per sq.in. 
and on 0.02-in. sheet with a yield strength of about 
180,000 Ibs. per sq.in. 

Comparison of the ARMCO method with the N.B.S. 
solid-guide method shows two principal differences. 
In the ARMCO method there is no check to determine 
that the initial friction is small. Also, the arrangement 
for clamping the guides is more rigid. This may lead 
to the development of excessive lateral pressure at high 
compressive strains, particularly with a thick speci- 
men, which may in turn produce relatively high fric- 
tional forces between the guides and the specimen. This 
may account for the high value of modulus (32,000,000 
Ibs. per sq.in.) for the 0.100-in. specimen that is re- 
ported in Table I of reference 14. It may also account 
for the relatively low strains for a given stress above 
70,000 Ibs. per sq.in. that are reported for single strips 
in Fig. 7 of reference 14 as compared to the strains 
measured on an unsupported solid cylinder from the 
same material. 

Two further differences between the two methods 
are the use of a positioning ring and of a shim of wet 
blotting paper in the ARMCO method. Our experi- 
ence with the N.B.S. solid-guide method indicates that 
good results can be obtained without these expedi- 
ents. 


N.A.C.A. Grooved-Guide Method 


Kotanchik, Woods, and Weinberger of the N.A.C.A. 
made a study* of compressive tests in which the speci- 
men was guided between lubricated brass plates with 
opposite and with offset vertical grooves, Fig. 2. The 
brass plates were held in a modified Montgomery- 
Templin fixture. For comparison they also made tests 
on specimens guided between rollers and between 
flat brass guides, as in reference 9, and between flat 
wooden plates. 

Tests made on 0.064-in. aluminum-alloy sheets 
showed good agreement between the results obtained 
with the various types of support as long as the initial 
lateral pressure between the guides and the specimen 
did not exceed 45 Ibs. 

In the tests with the offset-grooved guides the values 
of Young’s modulus were in close agreement with those 
obtained in the test with rollers for the whole range of 
initial lateral pressure (22 to 380 Ibs.). However, the 
indicated yield strengths increased with pressure above 
45 lbs. This indicates that an appreciable fraction of 
the load was by-passed as a frictional force between 
the guides and the specimen. Such a frictional force 
may be expected to arise in a relatively rigid retaining 
clamp, such as the Montgomery-Templin fixture, since 
the specimen tries to expand laterally at high compres- 
sive strains. This brings out again the advisability of 
introducing a certain amount of flexibility in the clamps 
to prevent the building up of excessive lateral pressures 
at high compressive strains. 


Van den Broek’s Solid-Guide Method 


Van den Broek has recently described a method of 
testing single strips of thin sheet in compression between 
solid guides in which a simple and ingenious means is 
employed to keep the lateral pressure nearly independ- 
ent of the load on the specimen." 

A sketch of Van den Broek’s fixture is reproduced in 
Fig. 2. The specimen, S, is clamped between the 3- 
by 1-in. flat surfaces of two cast-iron guides, G. No 
lubricant is used between the guides and the specimen. 
The clamping pressure is applied through a pair of 1-in. 
square bars, B, which are 40 in. long and which are held 
together at the ends by C-clamps, C. The C-clamps 
are tightened by trial and error to develop the minimum 
clamping pressure that will prevent buckling of the 
specimen before the yield strength is reached. The bars, 
B, are of sufficient flexibility so that the variation in 
clamping pressure with increasing compressive strain is 
negligible during a given test. 

In the absence of lubricant between guides and speci- 
men an appreciable portion of the load applied by the 
testing machine may be by-passed through the guides 
by friction. A correction to each stress-strain curve is 
made by assuming that the frictional error develops 
very fast during the early stages of the test and then 
remains constant during the test. From the illustra- 
tion given it would seem that such a correction could 
not be made accurately for Alclad and other sheet 
where the stress-strain graph begins to curve at low 
stresses. It is not clear from reference 15 what steps 
are taken to prevent stressing of the specimen by the 
weight of the fixture, which is estimated to be about 25 
Ibs. 

No comparisons are given with stress-strain curves 
obtained by other methods. The method in its present 
state of development does not seem well adapted to 
routine testing. ; 

Van den Broek’s method has not been tried at the 
National Bureau of Standards. If it should be tried, an 
attempt would probably be made to reduce the friction 
between the guides and the specimen with a lubricant 
as in the N.B.S. fixture. 


A.T.S.C. Solid-Guide Method 


Major B. C. Madden of the Air Technical Service 
Command at Wright Field has developed a method of 
testing single-sheet specimens in compression between 
lubricated steel guides, which is described in reference 
16. 

In the A.T.S.C. method the specimen, S, Fig. 8, is 
1.800 + 0.010 in. wide and usually 5.500 + 0.005 in. 
long. It is supported against buckling by being clamped 
between the ground faces of two guide plates, which 
are spaced by shims to leave a total clearance of 0.0005 
to 0.005 in. between the specimen and the plates. This 
clearance is filled with lead soap base grease. Each 
plate has two 0.12- by 0.80-in. apertures to provide ac- 
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Fic. 8. A.T.S.C. fixture. 


cess for l-in. gages to the center of both faces of the 
specimen. The dimensions of these windows are made 
small to prevent buckling of the specimen before the 
yield strength has been attained. Strain is usually 
measured with a pair of 1-in. Tuckerman strain gages. 
The load is applied through a spherically seated bearing 
block acting on the top end of the specimen through a 
0.040- by 1- by 2-in. piece of relatively soft metal. 
The amount of friction between guiding plates and 
specimen is limited by requiring that the specimen 
should slide under a load of 5 Ibs. at the beginning of 
the test. With this precaution it appears from tests 
reported in reference 16 that friction at the maximum 
load is less than § Ibs. 

Reference 16 shows close agreement between stress- 
strain curves obtained by the A.A.F. method and by 
the pack method on hard-rolled 0.035-in. manganese 
steel. 


ao 
or 


The minimum thickness of specimen that can be 
tested to its yield strength is limited by the tendency to 
develop a buckle at the apertures for the strain gages. 
It is estimated in reference 16 that the minimum thick- 
ness is about 0.02 in. 

Since the specimen is much wider than the specimens 
of other methods, relative tilting of the heads causing 
differences in strain at the two edges would cause greater 
inequalities of strain, for the same tilt, than for other 
methods. No measurements of this inequality can be 
made with gages at the center of the width. Difficulties 
in obtaining proper contact which are encountered for 
thin specimens might be avoided by the use of a cylin- 
drically seated bearing block in place of the spherically 
seated block. 


PRESENTATION OF STRESS-STRAIN CURVES 


Compressive stress-strain curves must usually be 
known with greater accuracy than tensile stress-strain 
curves because of the dependence of buckling loads at 
high stresses on the shape of the knee of the stress- 
strain curve. In estimating buckling loads it is impor- 
tant to know the tangent modulus vs. stress curve—i.e., 
the slope of the stress-strain curve rather than the 
stress-strain curve itself. 

A rough approximation to the stress-strain curve can 
sometimes be constructed from a knowledge of Young’s 
modulus and yield strength by drawing a modulus line 
to the yield strength and continuing with a horizontal 
line. More is certainly needed for even the roughest 
approximation to the tangent modulus vs. stress curve. 
With new materials especially it will be necessary to 
obtain a complete stress-strain curve and tangent modu- 
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Fic. 10. Dimensionless tangent modulus versus stress graphs, 
0.064-in. R3801-T aluminum-alloy sheet. 


lus vs. stress curve by an accurate method. However, 
it would be unfortunate if such complete stress-strain 
curves were to be required for each new batch of a given 
material. 

Fortunately, it turns out! that the stress-strain 
curves of many materials, particularly materials of the 
same nominal composition and treatment, are affine. 
This means that the stress-strain curves can be brought 
into coincidence by multiplying the stresses by a suit- 
able factor and the strains by another suitable factor. 
These factors are preferably chosen to reduce the ordi- 
nary stress-strain curve to a dimensionless stress-strain 
curve in which the effects of variations in Young’s modu- 
lus and in yield strength are eliminated. 

Fig. 9 shows the dimensionless stress-strain graphs 
for two longitudinal and two transverse specimens of 
0.064-in. R301-T aluminum-alloy sheet. The coordi- 
nates o, ¢ of the stress-strain curve have been replaced 
by the dimensionless coordinates o/a9.7, €£/oo.7, where 
E is Young’s modulus and where o.7 denotes the secant 
yield strength defined as the ordinate of the intersec- 
tion with the stress-strain curve of a line through the 
origin with a slope 0.7E. The conventional 0.2 per cent 
offset yield strength could not be used for the reduction 
since it will not reduce affinely related stress-strain 
curves to a common curve. The secant yield strength 
for a slope of 0.7E was chosen because examination 
of the stress-strain curves of reference 1 indicated that 
the secant yield strength for this slope comes near to the 
conventional 0.2 per cent offset yield strength for a 
number of high-strength sheet materials. 

A dimensionless curve of tangent modulus, E’, vs. 
stress may be derived directly from the dimensionless 
stress-strain curve by plotting the slope, E’/E, of the 
curve against o/oo.; (Fig. 10). From this curve, in 


turn, one may compute a dimensionless curve of the 
reduced modulus for a rectangular section: 
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Fic. 11. Dimensionless reduced modulus versus stress graphs, 
0.064-in. R301-T aluminum-alloy sheet. 


E/E = 4E'/E(1+WVE'/E)? (1) 
This is shown in Fig. 11. 


EXAMPLE FOR USE OF DIMENSIONLESS TANGENT- 
MOoDULUS AND REDUCED-MODULUS CURVES 


An example will be given to illustrate the applica- 
tion of these dimensionless tangent-modulus and re- 
duced-modulus curves to estimate the effect of varia- 
tions in Young’s modulus and in yield strength for ma- 
terials with affinely related stress-strain curves. 

A structure fabricated from 0.064-in. R301-T alu- 
minum-alloy sheet buckles under compressive loads 
acting in the longitudinal direction at a stress of 


o = 55,000 Ibs. per sq.in. 


Compressive tests of the sheet give the following prop- 
erties for the longitudinal direction 


E = 10.8 X 10° lbs. per sq.in. 
go.7 = 62,000 Ibs. per sq.in. 


The buckling stress is to be estimated for similar mate- 
rial having ‘‘minimum properties”’ of 


E = 10.4 X 108 lbs. per sq.in. 
go.7°= 57,000 Ibs. per sq.in. 


‘ 


It is assumed that the stress-strain curves in longitu- 
dinal compression of the test material and the hypo- 
thetical material with minimum properties are affinely 
related. The buckling stress will be computed under 
two assumptions (see, for instance, reference 18, pp. 
384-390). For the first computation it is assumed that 
the buckling stress is proportional to the tangent modu- 
lus and, for the Second case, that it is proportional to the 
reduced modulus. In the first case then 


o/E’ = const. 


and in the second case 
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Fic. 12. Auxiliary curves for estimating changes in buckling 
stress with changes in o/oo.7, 0.064-in. R301-T aluminum-alloy 
sheet in longitudinal compression, 


o/E, = const. 
The test material will be denoted without an asterisk 
and the hypothetical material with minimum proper- 
ties by an asterisk. 
For the test material 
a/oy.7 = 55,000/62,000 = 0.89 


Carrying through the computation for the first case, 
one obtains from the dimensionless curve of tangent 
modulus in longitudinal compression, Fig. 10, 
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E’/E = 0.48 
Hence, the tangent modulus for the test material is 
E’ = 0.48 X 10.8 X 10° = 5.18 X 108 Ibs. per sq.in. 
and the fixed ratio 
o/E’ = o*/E’* = 55,000/5.18 X 10®° = 0.01062 
For the hypothetical material 
oy.7*/E* = 57,000/10.4 X 10° = 0.00548 
Combining this ratio with the previous ratio gives the 
following equation for o*: 
oo7*/E* E'*/E* _ 0.00548 


— a = 0.516 
o*/E’* 








~ @*/ay2* 0.01062 


This equation may be solved for o*/oo.;* by drawing 
a line (shown dotted) through the origin of Fig. 10 
with a slope of 0.516 and determining the abscissa of 
intersection with the modulus curve for longitudinal 
compression. Or it may be solved very simply with the 
auxiliary curve for (E’/E)/(o/oo.7) in Fig. 12. 
The following value is obtained 
o*, 00.7" = 0.895 
The estimated buckling stress is therefore 


o* = 57,000 X 0.895 = 51,000 Ibs. per sq.in. 


For the second case, one obtains analogously from 


Figs. 11 and 12 


E,/E = 0.67 
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Dimensionless compressive stress-strain curves having markedly different shapes. 
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Fic. 14. 


E, = 0.67 X 10.8 X 10° = 7.24 X 10 Ibs. per sq.in. 
o/E, = o*/E,* = 55,000/7.24 X 10° = 0.00760 
oo.7*/E* = 57,000/10.4 K 10® = 0.00548 

oo.7*/E*  E,*/E* _ 0.00548 


o*/E,* 


=——— = 0,721 
o*/a0.7* 0.00760 
o*/oo.7* = 0.898 


o* = 57,000 X 0.898 = 51,200 Ibs. per sq.in. 


which happens to be only slightly larger than the esti- 
mated buckling stress for the first case. 

“jit should be mentioned that the above procedure is 
just one of several ways of correcting buckling stress for 
changes in properties, all of which are based on the as- 
sumption of an affine relationship between the stress- 
strain curves. Karol outlines a convenient procedure of 
this type in reference 19 and actually gives charts of 
correction factors for Alclad 24S-T and Alclad 75S-T 
aluminum-alloy sheet. 


CLASSIFICATION OF STRESS-STRAIN CURVES 


As mentioned above, the primary purpose of the di- 
mensionless plotting is to reduce the stress-strain curves 
to a form in which the effects of variations in Young's 
modulus and in yield strength are eliminated, so that 
only the deviations from affinity—that is, the changes 
in shape of the stress-strain curve—remain. 

This suggests a classification of stress-strain curves 
by shapes that may vary between the extremes, illus- 
trated in Fig. 13, of gradual curvature beginning at low 
stresses, and of practically no curvature up to the yield 
strength. An examination of stress-strain curves of 


|.5 2.0 €E/oz 


Dimensionless stress-strain curves described by Eq. (3). 


aluminum-alloy and of alloy-steel sheets with widely 
different shapes, given in reference 1, shows that most of 
them can be described by Osgood’s formula :!” 


« = (¢/E) + K(o/E)" (2) 


in which K and are constants that are adjusted to 
give an optimum fit to the stress-strain curve. The 
meaning of these constants is made clear by writing Eq. 
(2) in terms of the dimensionless variables o/05.7, 
e/oo.7. It is shown in reference 17 that this leads to 


eE/o0.7 7 (go, 00.7) + of ‘*(a, ‘O0.7)” (3) 
Comparison of Eq. (3) witli Eq. (2) shows that 
K = 3/; (E/o0.7)"™ 


i.e., K depends on the Young’s modulus and the secant 
yield strength of the material. Only the shape parame- 
ter m remains in the dimensionless Eq. (3). Fig. 14 
shows the effect of variations in m on the shape of the 
stress-strain curve. 

The applicability of Eqs. (2) and (3) may be tested 
by plotting the stress-deviation curve (deviation = e — 
a/E) on a log-log scale. This should lead to a straight 
line if Osgood’s formula holds, since 


log [e — (¢/E)] = log K + mn log o/E (4) 


The values of K and of m may be determined from the 
ordinate at a given abscissa and from the slope of this 
straight line. 

Approximate values of m can be obtained quickly, for 
stress-strain curves described by Osgood’s formula, by 
determining the ratio of two secant yield strengths as 
explained in reference 17. Fig. 15 shows the relation 
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between n and o».7/0.35, Where o.3,; denotes the secant 
yield strength for a line with slope 0.85E. 

A comparison of measured stress-strain and tangent- 
modulus curves with curves obtained from Osgood’s 
formula” shows that the formula gives an adequate 
and simple description for many stress-strain curves. 
It follows that ” is a convenient parameter for classify- 
ing such stress-strain curves. 

Osgood’s formula is used in reference 20 to derive a 
nomogram for estimating the effect of changes in 
Young’s modulus and in yield strength on the buckling 
load of a structure. This provides a particularly rapid 
solution for materials with affinely related stress-strain 
curves described by Osgood’s formula. 

It should be pointed out that the formula will not 
be adequate for clad materials for which the stress- 
strain curves show a change in slope when the clad- 
ding yields at low stress, followed by a linear por- 
tion up to the region in which the core begins to yield. 
In such cases a useful analytical description of the 
stress-strain curve can sometimes be obtained by divid- 
ing the stress-strain curve into two portions and de- 
scribing each portion independently by Eq. (2). . 


CONCLUSIONS 


Accurate values of the compressive stress-strain 
curve up to stresses of 60,000 Ibs. per sq.in. may be ob- 
tained on sheet with a thickness not less than 0.02 in. 
by any one of the following methods: pack method, 
Montgomery-Templin method, N.B.S.  solid-guide 
method, N.A.C.A. grooved-guide method, ARMCO 
solid-guide method, and A.T.S.C. solid-guide method. 

The stress-strain curves for much stronger materials 
and for thinner gages may be obtained by several of 
these methods if proper attention is given to the prep- 
aration of the specimen and to details of the testing 
procedure. 
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The effects of small variations in Young’s modulus 
and in yield strength on the stress-strain curve and on 
the buckling stress in the plastic range can be estimated 
frequently from dimensionless curves of stress vs. 
strain and of tangent modulus and reduced modulus vs. 
stress. 

A large class of stress-strain curves, which, however, 
does not include all clad aluminum alloys, can be de- 
scribed by a simple formula proposed by Osgood which 
requires a knowledge of only one parameter in addition 
to Young’s modulus and a secant yield strength. A 
convenient set of three such parameters is Young’s 
modulus and the secant yield strengths corresponding 
to slopes of 0.7E and 0.85E. 
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Letter to the Editor 


Dear Sir: equation 4.16 of von Karman and Biot’s book seem to be an Ex 

The following is with reference to Alexander Flax’s letter of | ¢Xtrapolation on the von Karman and Biot conclusion which is: 7 
July 8, 1946, appearing in the October, 1946, JouRNAL and com- “ Therefore, to reduce the stress in a long pipe materially by a empii 
menting on my paper, ‘“‘Notes on the Design of Pressurized Air-  S€ries of stiffener rings the spacing must be of the order of magni- Sion | 


‘ gow J /atfa 9 
tude of 2] = 2 v/rt vV/3. coeffi 


craft” (June, 1946, issue of the JOURNAL). 
(3) In reply to Mr. Flax’s comment No. 3, ‘‘Table 1 and Fig. 2 


(1) Mr. Flax’s comments on the assumptions concerning the streay 
altitude variation of load factor do not agree with some of the of the paper for the stress coefficients of a rectangular membrane efficic 
basic philosophy underlying aircraft design criteria—namely: appear to be in error,’’ the footnote on page 303 of the writer’s tho ve 

(a) The maximum speed at which aircraft are flown is Paper, as well as Fig. 1, shows that the dimensions (a) and (b) with. 
legally limited by placard and not by the terminal velocity that re twice those used by Timoshenko. A simple numerical com- tropic 
the airplane-may attain in a dive: Such placards are used to Putation will show that the coefficients that he has listed are ateacth 
limit maximum horizontal speed, as well as diving speed; altered considerably by the fact that Timoshenko’s a is !/2 that —_ 
hence, aerodynamics and engine performance do not form the used by the writer*—Mr. Flax has also omitted any statements result 
complete basis for load factor determination. It should also tegarding the value of Poisson’s ratio used in his computations. for ca 


(4) Regarding comment No. 4, the writer is well aware of the of th 


be realized that discussions of load factors cannot be conducted 
literature on the effects of internal pressure on buckling of skin 


purely on an academic basis, since in the case of the maneuver trente 
factors for commercial aircraft the load factor attained de- Under shear and compressive loads. It is pointed out, however, 
pends mainly on the discretion of the pilot. There is probably that most of the literature is ‘‘restricted” and therefore should 
no transport aircraft in existence today which cannot be torn N0t be quoted. 
apart in the air in the hands of a reckless pilot. (5) Mr. Flax’s comment Nos. 5 and 6 seem to have little bear- | I | 
(b) The maneuver factor prescribed by Civil Air Regulations ing on the illustrative example presented in the writer’s paper, 
for large aircraft is generally smaller than the gust factor and, since Mr. Flax has not seen fit to present a comparative set of prob 
consequently, is secondary for design purposes. In this regard figures demonstrating how the “critical structural altitude”’ is Since 
the arbitrary load factor reduction formula suggested by the altered by the assumptions he proposes. It is agreed, however, d 
writer is conservative, since it applies only a density correction that a great deal of additional work can and should be done to an 
on the gust factor and does not include the effects of the actual determine effects of pressurization on such items as stringer loads prese 
réduction of effective gust velocity (see discussion in paper on introduced by tension field and similar buckling formations. the f 
gust loads). The writer’s paper was written as an introduction to the problem knov 
(2) Mr. Flax’s comments on hoop tension are interesting, but | 4d was not intended to be a comprehensive treatise. mad 
‘ sly academic point of view. » has convenien 
from a purely academic point of vie ’ He has co veniently Oh ieeeen as 
ignored the writer’s comments (page 303 of June, 1946, JoURNAL) B ; . ; 
; neh ~ ell Aircraft Corporation load 
on skin continuity. Perhaps Mr. Flax’s comments would be 
legitimate for the design of a pressurized tai] cone without cutouts MELE. ics, Pecans : ; wat the s 
. 3 é : * For comparison purposes Mr. Flax should have multiplied his coefficients 
and without skin seams. Further, his statements based on py 22/3 = 1.588. In 
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Additional Developments of Greene’s 


Equation to the Study of Compressible 
Airflow 


GEORGE 
Const rlidated J ‘ullee 


SUMMARY 


Extension of Greene’s ‘‘Attenuation Method for Compressible 
Flow Systems” to include the supercritical region leads to an 
empirical interpretation of supercritical conditions. The correc- 
tion that is then applied to the airfoil incompressible pressure 
coefficient predicts supercritical flow conditions on the down- 
stream side of the peak negative incompressible pressure co- 
efficient. The calculated position of the shock that accompanies 
the return to subsonic conditions is shown to progress downstreain 
with increasing Mach Number. With a suitable choice of a poly- 
tropic flow process the calculated chordwise pressure distribution 
checks with high-speed test data and this automatically brings 
the lift and pitching moment values into agreement with test 
The energy loss in the shock yields a reasonable means 
Each surface 


results. 
for calculating the increase in drag at high speeds. 
of the airfoil and each item in the airplane configuration are 
treated separately in estimating the drag. 


INTRODUCTION 


_— PROBLEM of accurately predicting the behavior 
of an airfoil at high speeds is one of the serious 
problems confronting the aerodynamics engineer today. 
Since the density of air is variable with the temperature 
and pressure, it invalidates incompressible theory at 
present-day speeds. The effects of compressibility on 
the forces produced by airfoils at high speeds must be 
known if accurate predictions of performance are to be 
made. The moments of these forces are important 
factors in stability and control problems. The total air 
load and its distribution are extremely important to 
the structures engineer. 

In the light of our present knowledge, the most 
practical way of accounting for compressibility is to ex- 
press the relationship of the pressure coefficients for 
compressible and incompressible flow systems in some 
rational manner. The wealth of incompressible pres- 
sure distribution data and theory are then useful. The 
corrections that are applied to pressure coefficients 
reflect in forces and moments upon integration of the 
corrected pressure distributions. If it be remembered 
that the fundamental corrections are applied to pres- 
sure, then more reliable corrections, when they become 
available, will be applied to all aerodynamic en- 
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NOTATION 
dy = slope of lift curve for section (low Mach Number) 
a = speed of sound 
A = aspect ratio 
Ay =(W1— Gy — 1/(1 — Cr) 
b = span 
B? = M? + [2/(n — 1)](1 — F*""?) 
c = chord 
Cp = drag coefficient 
Cr = lift coefficient 
Cy = pitching moment coefficient 
C, = pressure coefficient Ap/q 
F = (1 + '/onM?Cpy)'’" = density ratio 


K = Cp er, Cro 
M = Mach Number 
Mr, = critical free stream Mach Number 
m = mass of air affected per sec. 
= polytropic index 
“EE = kinetic energy per sec. 


= wing area 
= density 


n 
K 
r = radius 
¢ 
p 
n =m/paMS 


Subscripts 


0 = incompressible 

M = compressible 

l = limiting local conditions 

2 = local conditions just before shock 
3 = local conditions just after shock 
4 = local conditions 

none = free stream conditions 


GLAUERT CORRECTION 


The commonly used correction for compressibility is 


that of Glauert V1 — M? (reference 1). It was de- 
veloped for application to two-dimensional pressure 
data. However, it is found that the same factor appears 
in the forces and moments. For this reason and also 
because it is easy to use, the Glauert factor is often 
applied in places where it is not sufficiently accurate or 
where it does not apply. For example, the slope of the 
lift curve for both the airfoil section and for a complete 
wing are almost always corrected by Glauert’s factor. 
If the wing lift curve slope is properly corrected by 
application of the factor to section data the resulting 
slope is given by 


a. 
da V1—M? + (a)/rA) 
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KARMAN-TSIEN CORRECTION 


The Karman-Tsien correction? is becoming more 
widely used and gives better agreement with experi- 
ment than does the Glauert correction. It applies to 
pressure distribution for the section and care should be 
exercised when it is applied to wings of finite aspect 
ratio. To obtain the corrected force coefficients one 
must integrate the corrected pressure distribution at 
each Mach Number selected. It gives excellent agree- 
ment with experiment at subcritical Mach Numbers. 


M? Cc. 
Cou = Gall VIB ae) # | 


GREENE’S ATTENUATION METHOD 





The attenuation method of Greene* is superior to 
either of the previously mentioned corrections. The 
author’s reasons for believing that Greene’s equation is 
the best available are: 

(1) It is sufficiently accurate for most engineering 
purposes. 

(2) It permits the use of any desired polytropic 
process including a variable index. 

(3) Itis easily adapted to the correction of pressure 
distribution for bodies of revolution. 

(4) It gives an empirical interpretation for pressure 
distribuiion at supercritical Mach Numbers. 

(5) By judicious choice of the polytropic index it 
can be made to agree with experiment. 

The author’s interpretation of the attenuation 
method and its practical applications are the principal 
subject matter of this paper. Agreement with test 
data in the supercritical region is the prime considera- 
tion. 


INTERPRETATION 


The word ‘‘critical’’ is used with several possible 
meanings. The author suggests that ‘‘critical’’ be re- 
tained as a generally descriptive term and that “‘sonic’”’ 
be used to describe conditions in which the local speed 
of sound is reached. Similarly, “‘limiting’’ is used to 
describe conditions at which shock losses begin to ap- 
pear, and ‘‘absolute’’ describes the minimum possible 
static pressure (theoretically zero) conditions. Then 
sonic Mach Number would mean the free stream Mach 
Number at which the speed of sound is reached locally, 
and limiting Mach Number would designate the Mach 
Number above which shock losses occur. In the con- 
sideration of drag, critical usually indicates the Mach 
Number at which the rise in drag begins. 

The equation for Greene’s attenuation method may 
be written 


dCyo/dCpm = M(FB — M)/A,FB? (1) 
The local Mach Number is given by 
M, = B/F*- 9” 








1946 


When M, = O is solved for Cy it yields the well-known 
equation for the pressure coefficient at the stagnation 
point 


D) n/(n — 1) 
Coun ah +(n- p/2ae| - i} (3) 


When the local speed is equal to the free stream speed 
it is found that Cj is zero by definition. When the 
local Mach Number is unity, Eq. (2) yields 


2 ([2+ (n—1)Mey/e-» \ 
ail n+1 | 7 . (4) 


The limiting conditions without shock losses are de- 
fined by FB = M. It can be shown that the local 
stream tube cross section is the same as in the free 
stream. 

The values for C,y for this case at selected Mach 
Numbers are found in the process of integration 
(see Fig. 4a, reference 3) or by Eq. (25) of reference 3, 
When the density ratio is zero the absolute conditions 
are indicated and 


Cya(abs.) = 2, ‘nM? (5) 








Cym(sonic) = 


Numerical integration of Eq. (1) with a variable poly- 
tropic flow process is accomplished easily by the dif- 
ference table method.‘ The author’s interpretation re- 
quires that the curves be plotted with Cy, and C)y 
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tropic index that brings the theory into agreement with 
experiment varies with the type of airfoil. The NACA 
4412 airfoil appears to check when = 1.7 — 0.3 Gy. 
Low-drag sections appear to require a variation of ” 
similar to that illustrated in Fig. 4. It is probable that 
the relative humidity of the wind tunnel is an important 
factor in the type of polytropic process required. Con- 
densation shocks cause a reduction in the peak negative 
C,a observed as compared to theory. 
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NUMERICAL APPLICATIONS 


The curves in Figs. 2a and 2b are typical results from 
high-speed wind-tunnel data. It is significant that the 
curves deviate from the origin by increasing amounts 
as the Mach Number increases, and that following the 
shock the curve passes close to the origin. The local 
speed equal to the free stream speed occurs at different 
chordwise stations as the Mach Number is changed. 

The author is aware of the observation that in the 
supersonic region the airfoil shape is of fundamental 
importance and that C,,, in this region is not necessarily 





Fic. 2a. Cpm vs. Cpo. Test data. 





NACA 4412. 
Fig. 3. 


Fic. 2b. Cymvs. Coo. Test data. 
a = 9; M = 0.756; reference 5, 
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conditions. 


a function of Cj». However, Cy,» is a function of the air- 
foil shape. 

The author’s empirical interpretation for supercritical 
conditions is illustrated in Fig. 3. The curve on the left 
(A) Cy vs. Cyo is obtained by integration of Eq. (1). 
The solid curve on the right (B) is the given incompres- 
sible pressure distribution. A new scale for Cy is made 
for the curve on the left so that the required maximum 
negative C,o matches the limiting conditions. This 
scale is illustrated above the curve. Near the trailing 
edge the values of C,,, which correspond to each Cyo 
are found for the subsonic region using the regular Cyo 
scale. Beginning at the stagnation point the chord- 
wise flow is followed, finding the values of C,, for each 
C,o on the new scale. When C, passes the peak nega- 
tive point the larger negative value of C,y is used. 
The shock is drawn at the chordwise station so that it 
terminates on the local sonic speed. 

If the shock is severe, the energy losses raise the tem- 
perature of the air and reduce the effective free stream 
Mach Number. In cases where severe shocks are en- 
countered the region following the shock is revised, using 
a reduction of approximately 0.03 in the free stream 
Mach Number. This moves the shock forward slightly 
to the new C,,, (sonic) position. 


ADIABATIC vs. POLYTROPIC FLow 


An adiabatic flow process gives unsatisfactory agree- 
ment with experience. This is true for both the sub- 
critical and supercritical cases. Since the K4arman- 
Tsien correction gives excellent agreement with test 
data in the subsonic region, a value of — 1 is selected for 
n to apply where the flow is subsonic. (When —1 is 
used for n, Greene’s equation is substantially equivalent 
to Karman-Tsien.) In the region close to the stagna- 
tion point an adiabatic (m = 1.4) process is indicated 
since there is no way to change the entropy. 
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By checking the peak negative C,,, from all avail- 
able data, the envelope was assumed to represent the 
practical values of Cj, (abs.). The corresponding 
values for m at test Mach Numbers fell on a smooth 
curve which could be represented by 

2 2 
n= ——— TF — _— (6) 
M?*C,m(abs.) 1.44 — M 

Fig. 4 illustrates the currently recommended values of 
n for low-drag sections. A value of 1.4 is always used 
for m at the limiting conditions. A comparison of the 
adiabatic and variable polytropic curves is shown in 
Fig. 5. 


COMPARSION WITH TEST RESULTS 


The curves in Fig. 6 show a comparison of the pre- 
dicted pressure distribution for a 4412 section and the 
measured pressure distribution from reference 5. Fig. 
7 illustrates the comparison of the calculated lift and 
pitching moment with experiment at high speed. When 
the pressure peak is near the leading edge the estimated 
pressure distribution for supercritical flow conditions 
does not agree with test data. It appears probable 
that local separation caused by centrifugal force is the 
explanation. The conditions at which the normal 
pressure gradient just balances the centrifugal force 
are expressed by 


dC,/dr = 2psVs2/rpV? . (7) 


CORRECTION OF PRESSURE ON 
THREE-DIMENSIONAL BODIES 


Greene (reference 3, Eq. 27) has adapted his equation 
to bodies of revolution. Unfortunately, this adaptation 
requires a new integration for each fineness ratio. In 
most engineering work the accuracy required and the 
time available do not indicate such refinements. Fre- 
quently, the bodies used are more nearly ellipsoids 
than bodies of revolution. An approximate correction 
can be obtained by using an appropriate fraction of the 
two-dimensional correction 


b 


Coa —_ v4 
4 a+b 


“ Com (2-dim.) (8) 


a 
Cro + a+b 


where a and 8 are appropriate major and minor diam- 
eters. If a body of revolution is corrected in this manner 
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Fic. 6. Comparison of test and calculated pressure distribution. 
NACA 4412. a = —2°. 
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Fic. 7. Calculated and measured airfoil characteristics at high 
speed. 
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ADDITIONAL DEVELOPMENTS OF GREENE’S EQUATION 


it amounts to using one-half of the two-dimensional 
correction. 


ESTIMATION OF DRAG 


The drag that is added because of the shock can be 
estimated by the method presented below. For drag 
estimates each of the two surfaces of the airfoil should 
be separately calculated, using the appropriate critical 
Mach Numbers. In this case “‘critical’’ refers to ‘‘limit- 
ing Mach Number”’ or the Mach Number at which the 
drag begins to rise sharply. The drag coefficient for a 
complete airplane can be found by use of the same 
formula by calculating the contribution of each compo- 
nent of the airplane separately. Since each item of the 
airplane configuration is separately calculated, it is easy 
to estimate the effect of any particular item on high- 
speed performance. 


Formula for Drag 


The kinetic energy change in the shock is assumed to 
be represented by the increase in drag. Then 


VAD = KE, — KE; = ACpS'/»a°?M3 = 
1/sm(a2?M,? — a37.M;3") (9) 


*/gm (a2°M2? — as°M3*) _ 


KE, — KE 
ACp = ———* 





S'/spa°M*® ———«S"/spa*M*a M 
m a2*Mz? — as°M3" 


pa MS  @?M? 





(10) 


Borrowing an idea from Glauert, it is assumed that the 
square of the limiting velocity is equal to the product of 
the velocities immediately before and after the shock. 


ay?>M,? = a2M2a3M; (11) 

In order to have symmetry in the equation, let 
a2M, = (M/M,.)aM, (12) 
a;M; = (M.,,/M)aM, (13) 


Then the increase in drag coefficient becomes 
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Now m/paM is the cross-section area of stream tubes 
(in the free stream) which pass through the shock 
zone. For a wing, this area is some fraction, f, of the 
product of the chord and the span, or 

m/paMS = fcb/paMcb = » (15) 


The drag coefficient at Mach Number then may be ex- 
pressed in terms of the low-speed drag coefficient. 


a2M,?}/M? Me? 
Ot Qi + - > 
—_ | TT owe (S M2 )| 


2+ (n — 1)M? M? x 
2+ (n — 1)M,? M 


(z Ma tia 
M2? M? sii 


The values of K and » are found from high-speed tests 
and are approximately 


K= 





a Cook E + | 


1.20 
1.0M?2/Me? 


u] 


Substitution of these values in Eq. (16) gives: 
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Fic. 8. Drag correction—high speed. 


TABLE 1 


m a,"M,? M \ M...\? ; 
scp = MeL (MY _ (Me) ] ay 
paMS a2? L\M.. M 
Part Ma. M=0 M=0.46 
Wing 
U.S. (inner) 0.65 0.0025 0.0028 
Wing 
U.S. (outer) 0.68 0.0020 0.0022 
Wing 
L.S. (inner) 0.70 0.0025 0.0027 
Wing 
L.S. (outer) 0.70 0.0020 0.0022 
Fuselage 0.75 0.0045 0.0047 
Tail 0.70 0.0022 0.0024 
Power Plant 0.65 0.0015 0.0017 
Misc. 0.60 0.0010 0.0012 
Total se 0.0182 0.0199 


M=065 M=07 M=075 M=08 M=085 M=09 
0.0030 0.0061 0.0088 0.0110 0.0128 0.0144 
0.0023 0.0032 0.0051 0.0066 0.0081 0.0091 
0.0028 0.0030 0.0053 0.0072 0.0088 0.0103 


0.0023 0.0024 0.0043 0.0058 0.0070 0.0082 
0.0049 0.0051 0.0054 0.0085 0.0112 0.0135 
0.0025 0.0026 0.0047 0.0064 0.0077 0.0090 
0.0018 0.0037 0.0053 0.0066 0.0077 0.0087 
0.0028 0.0043 0.0054 0.0064 0.0072 0.0079 


0.0224 0.0304 0.0443 0.0585 0.0765 0.0811 
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Fic. 9. Cp vs. M. Assumed airplane. 


G 5+M? Mi? (MM? My? 
——™ = 1.20 E + otal - — a) (17) 
Coo 5+ Me Ma?\Me? MM? 


Eq. (17) is valid only when the critical Mach Number 
is exceeded. It is observed that no drag increase occurs 
below M/M,,. = 0.6. Fig. 8 contains the results of 
Eq. (17) in a usable form. For more accurate evalua- 
tion of the values of K and 7, a systematic series of 
tests on symmetrical sections at zero lift are required. 
The critical speed for both upper and lower surfaces 
would then be the same. 

The drag of a complete airplane can be estimated by 
applying Eq. (17) to each part. Table 1 illustrates an 
example for an assumed airplane at zero lift and Fig. 9 





represents the final result. The critical Mach Number 
for each surface of the wing is required at each lift 
coefficient. The induced drag is assumed to be given 
by C,?/7A at all Mach Numbers. 


CONCLUSIONS 


(1) In the absence of high-speed wind-tunnel tests 
Greene’s attenuation method may be used to estimate 
pressure distribution at supercritical Mach Numbers. 

(2) By checking out the polytropic index selected 
against test data for similar sections, the attenuation 
method is made extremely useful for obtaining details 
not included in the test data. 

(3) Reasonably accurate preliminary estimates of 
performance, air loads, and weights for new designs 
can be made from data, including corrections for com- 
pressibility, by the attenuation method. 

(4) The calculation of the drag for each item sepa- 
rately is useful for estimating the effects of changes in 
airplane configuration on the high-speed performance. 
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Tests on Elbows of a Special Design 


JOHN J. HARPER* 
Georgia School of Technology 


SUMMARY 


This paper presents the results of a series of tests on rectangular 
elbows of the type proposed by Szczeniowski.! Six elbows were 
tested, including one conventional guide vane elbow designed in 
accordance with the recommendations in R. & M. No. 1773.% 
Four of the special elbows had an aspect ratio (A.R.) of 0.67, 
while the fifth had an A.R. of 1.5. 

The results show that under certain design conditions the 
special elbows show fewer losses than do conventional elbows 
with guide vanes. Increasing the aspect ratio of the elbow 
reduces the losses. 


INTRODUCTION 


i acntge THE USE of special elbows without guide 
vanes, a reduction of losses in many applications 
may be possible. This type elbow might be used to 
advantage on the carburetor air intake of airplane 
engines, as there would be no guide vanes to ice up. 
The manufacture of these special elbows should be 
relatively easy. It is also possible that elbows of this 
type could be extended to wind tunnels, where elimina- 
tion of corner guide vanes would mean a saving in 
power required. 

In view of these possibilities, a test program was in- 
stituted consisting of velocity surveys on five of the 
rectangular two-dimensional elbows and one conven- 
tional elbow with guide vanes. All of the above were 
90° bends. 


THEORY 


The mathematical equations for the streamlines of 
these elbows were developed by Szcezeniowski! from per- 
fect fluid theory as follows. The elbows are shaped in 
such a way as to fulfill the equation of continuity in a 
two-dimensional flow: 


d*y/dx? + (0%Y/dy") = 0 


where ¥ = constant and is the stream function and 
where x and y are rectangular coordinates. 

For the 90° rectangular elbows considered, Szczeni- 
owski! gives the following complex potential fulfilling 
the equation of continuity: 


ei the OP ee eitD(w au) + ef -D@ au) (1) 


wheres = x+iy,i= -Y—l,w=o+ iy. The sym- 
bol ¢ designates the velocity potential; y, the stream 
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Station. 


function; a, any linear parameter; and u, any constant 
freestream velocity. . 

It is desired that the total velocity, V, at any point, 
should not exceed the freestream or boundary value u, 
or V <u. Likewise no stream line should have singu- 
lar points. Solution of Eq. (1) gives 














‘ ' cot ¥ + cot kd tanh 
% I tan au au au + 
cot £ cot ¥ — tanh 4 
au au au 
: »Z + LN 92(cosh 22 + cos 2¥) (2) 
- au au au 
cot £. + cot £ tanh ¢ 
y 1 = au au au 
- = 5 tan _ 
. cot = cot — tanh bd 
au y au 
—_ 
v ( ¢ v 
=| —-—t+ZLN ¥2| cosh 2 — + cos 2— 
2 au \ au au (3) 
0 glcosh 2(@/au) + cos 20/au) 
u cosh 2(¢/au) + sin 2(~/au) 


In the above, ¢ is a variable parameter for an arbi- 
trary streamline ¥, which is a constant. By selecting 
values of ¥/au and assuming values of ¢/au, the stream- 
lines can be plotted. The following streamlines com- 
posed the walls of the elbows tested (Fig. 1): y/au = 





54/12 and —7/8, ¥/au = 5/12 and w/12, y/au = 
y Fau 
Y15530 CUT OFF 
2a wr 





























fra tk 
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Fic. 1. Diagram of various streamlines. 
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Fic. 2. 


Photograph of elbow with side removed. 


8 and .-(/8; A.R. = 1.5.) 


(Streamlines 





Fic. 3. Photograph of three of elbows tested. 





Fic. 4. Photograph of three of elbows tested. 





Fic. 5. 


Photograph of test equipment. 


117/24 and 0, and ¥/au = 32/8 and 7/8. These 
streamlines were plotted by means of Eqs. (2) and (3). 
From Eq. (4), the following maximum velocities were 


found to be: ¥/au = —7/4, Vaz. = ©; Y/au = 
—/8, Var. = 2.415u; ¥/au = 0, Viner. = 1.4144; 
y/au = m/12, Var. = 1.15u; and y/au = r/8, 
Vimar. = tb. 


Some of the properties of these streamlines as given 
by Szezeniowski are as follows: 

(1) In cutting off the elbows at */,7a from the origin 
(Fig. 1), the error made in the channel width is less 
than 1 per cent of the width at infinity. 

(2) Thestreamline Y/au = 1/2 (Fig. 1) is the outer- 
most boundary of the field since it does not carry 
around the bend before reaching infinity. 


(3) The streamline ¥/au = 37/8 is the outermost 
streamline without an inflexion point. 
(4) The streamline ¥/au = 7/8 is the innermost 


streamline along whose length the velocity never ex- 
ceeds the entrance velocity 1. 

(5) The streamline —7/4 is the innermost boundary 
of the field as the velocity becomes infinite at the 45° 
point on the streamline. 


APPARATUS AND TESTS 


Fig. 1 shows the possible streamlines, for all the 
special elbows, which were plotted from Eqs. (2) and 
(3). White pine was used to construct the elbows, 
which had sides of !/4-in. plywood. The interiors were 
given two coats of wood filler, rubbed down with steel 
wool, and finally given four coats of automobile lacquer. 
Each coat of lacquer was carefully rubbed down with 
steel wool. Fig. 2 is a photograph of one of the special 
elbows with one of the sides removed. 

All elbows had entrance and exit dimensions of 4.5 by 
3 in. Of the five special elbows, four had an A.R. of 
0.67, and the fifth had an A:R. of 1.5. The aspect ratio 
may be defined as the depth of the stream divided by 
the width of the stream channel (width being the 
distance between boundary streamlines)—i.e., A.R. = 
3/4.5 = 0.67. For the fifth elbow, the reverse is true, 
or A.R. = 4.5/3 = 1.5. It should be understood that 
in order to keep the entrance and exit areas of all 
elbows constant it was necessary to use a different 
length for each elbow. This means that the value of 
the linear parameter a, which appears in Eqs. (2) and 
(3), must be different for each elbow. From Fig. 1 it 
is seen that the distance between a pair of streamlines is 
(r/4)a and the cutoff point measured from the Y axis is 
3/gra. Thus any pair of streamlines used has a certain 
value of a which is necessarily different from that of any 
other pair. 

The conventional elbow had five aluminum-alloy 
guide vanes and was designed according to the method 
recommended by Patterson.” Figs. 3 and 4 show the 
set of six elbows labeled with their corresponding 
streamlines. 
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Fic: 6. Lengthwise distribution for elbow of streamlines 111/24 
and 0 at different box pressures. 


The test equipment consisted of a 4 by 4 ft. plywood 
pressure box fitted with a convergent nozzle having 
exit dimensions of 4.5 by 3in. Air was supplied by a 5 
hp. centrifugal blower equipped with a throttle valve. 
A standard pitot tube was used in conjunction with 
water-filled, reservoir-type manometers. Fig. 5 is a 
photograph of the test equipment. 

In order to ensure good distribution at the elbow 
entrance, preliminary velocity surveys were made on 
the nozzle in the plane of exit. One survey was made 
lengthwise down the centerline of the exit, and others 
were taken at three equal intervals across the short 
dimension (% in.) of the nozzle. These velocity 
surveys were made at three different box pressures 
(100, 200, 300 mm. water) to determine the effect of 
change in velocity on the nozzle distribution. The 
velocity distributions obtained were all geometrically 
similar; therefore, the nozzle was considered satis- 
factory for testing the elbows. The same procedure 
was followed for one of the elbows, but there was no 
apparent change in the similarity of the distribution at 
different box pressures (Fig. 6). 

The final pitot surveys were identical for each elbow 
tested. Static pressure in the box was maintained at 
300 mm. water. A total of four surveys on each elbow 
was made: one down the centerline in the exit plane 
and three crosswise surveys at equal intervals across 
the short dimension (3 in.) of the exit. 


RESULTS AND DISCUSSION 


Figs. 7 through 14 show the velocity distributions 
obtained with various streamline combinations. These 
distributions were plotted isometrically, and the inside 
of the bend was marked in each case. 

Since all surveys on the nozzle were made with no 
bend in place, it cannot be assumed that the nozzle 
losses or velocity remain unchanged with an elbow on, 
since the flow resistance is increased and the velocity 
through the nozzle decreases. Consequently, to get 
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Velocity distribution for nozzle at box pressure of 300 
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Velocity distribution for elbow of streamlines 57/12 
and —7/8. 
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Velocity distribution for elbow of streamlines 117/24 
and 0. 
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Fic. 10. Velocity distribution for elbow of streamlines 57/12 


and 7/12. 








Fic. 11. Velocity distribution for elbow of streamlines 37/8 and 
4/8. 
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Velocity distribution for elbow of streamlines 37/8 
and 7/8; A.R. = 1.5. 
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the true elbow losses, it is necessary to eliminate the 
losses chargeable to the nozzle. All results have been 
expressed as a per cent of the total velocity head 
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Fic. 13. Velocity distribution for conventional elbow with five 


guide vanes. 
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Fic. 14. Velocity distribution for elbow of streamlines 57/12 


and —7/8 with two guide vanes. 


available at the elbow entrance. From the three nozzle 
surveys made at 100, 200, and 300 mm. water box 
pressure, a curve of pressure loss was obtained for the 
nozzle versus quantity, cu.ft. per sec.; this log plot is 
shown in Fig. 17. The surveys made at the elbow 
exits include the nozzle losses, and the total pressure 
loss may be expressed as 


AP rota = Hs — */2pV? 


where Hs = box static pressure and 1/29? = dynamic 
head measured at elbow exit. Therefore, for a given 
quantity of flow through the nozzle-elbow combination, 
the pressure loss, AP, for the nozzle is obtained from 
the graph (Fig. 17). Then 


AP nen ~ AP rotat es AP nesses 


and the per cent loss expressed in terms of the available 
head at the elbow entrance becomes 
BP se mm. H,O x 100 


loss = “ 
mean dynamic head for nozzle (mm. H,0) 
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TABLE 1 Dag Gas 9 yy pot ~~ 
Elbow (1) (2) (3) (4) (5) (6) (7) ~— (8) 
AP Total AP 
Mean Pressure Elbow Total 

Dynamic 5 Drop APN Loss Head 

Head, g, Velocity, Quantity, Through — Loss Through Available Per 

at Exit Ft. per Cu.Ft. Elbow + Through Elbow at Elbow Cent 

Mm. H,O Sec. per Sec. Nozzle Nozzle Alone Entrance Loss 
a 202. 80 195.4 18.3 97.20 4.25 92.95 295.75 31.45 
la ~ ~ 
Sp 0 242.52 213.5 19.9 57.48 4.70 52.78 295.3 17.90 
co 247.80 216 20.15 52.20 4.80 47.40 295.2 16.06 
3x . - — 
38 251.00 217.2 20.35 49.00 4.90 44.10 295.1 14.94 
Conventional—5 guide vanes 252.28 218 20.45 47.72 4.95 42.77 295.05 14.5 
#2 AR. = 15 264.00 22 20.85 36.00 5.00 31.00 295 10.5 
53 T : ae me i m= ve 
19” sa es guide vanes 235.60 210.5 19.7 64.40 4.65 59.75 295.35 20.15 
Conventional bend—no vanes 167.00 177.2 16.6 133.00 3.80 129.20 297 . 20 43.5 
(1) From velocity surveys (4) 300 — (1) (7) 300 — (6) 
(2) V = 13.7 Vq mm. H,0 (5) From Fig. 17 (8) (6)/(7) X 100 
(3) Q = AV = 0.0936 X (2) (6) (4) — (5) 


Thus the losses chargeable to the nozzle itself are 
eliminated in the final results. 

Szezeniowski has pointed out that an elbow composed 
of the stream lines 37/8 and 7/8 is probably the most 
reasonable one, since the use of these streamlines tends 
to suppress any vortex formation in the throat of the 
elbow. As can be seen from the isometric plot, all the 
“shorter” elbows exhibit a marked separation of flow 
at the inner wall. Since the elbow of 37/8 and 7/8 is 
the “‘longest’’ one of the group, it should have better 
flow characteristics; the isometric plot and the results 
confirm this. 

The effect of aspect ratio was studied to a limited 
extent, in that the bend of streamlines 37/8 and 7/8 
had its aspect ratio increased from 0.67 to 1.5. In- 
creasing the aspect ratio reduces the losses. The tests 
on the conventional bend were included for comparison 
(Fig. 13). Likewise Fig. 14 is included; this is the plot 
of velocity distribution for the ‘“‘shortest’’ special 
elbow with two guide vanes installed. These vanes 
were calculated as two streamlines so placed as to give 
equal theoretical volumes of fluid flowing between 
them. These vanes did decrease the losses but, not as 
much as expected. As shown in Fig. 15, there is some 
divergence between the inner wall and the first vane. 
This may account for the comparatively poor results; 
however, lengthening the vanes on the exit side would 
probably prove helpful, since the vanes were cut off 
short. Fig. 16 is included for comparison, since it is a 
plot of the velocity distribution for the conventional 
bend with guide vanes removed. 

Table 1 lists the results obtained for the various 
elbows. From these results, a few suggestions for 
future study can be made. It is quite probable that 
bends with guide vanes and still higher aspect ratios 
would have lower losses than the elbows tested, since 
the special elbows were calcualted from two-dimen- 





Photograph of elbow of streamlines of 51/12 and —7/8 
showing divergence of guide vanes. 


Fic. 15. 
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Velocity distribution for conventional bend with 
guide vanes removed, 


Fic. 16. 


sional flow and any skin friction on the walls will be 
detrimental to the realization of the calculated flow. 
Also, there are an infinite number of possible streamline 
combinations which will form the walls of the elbows, 
and it is likely that combinations other than those 
tested would give better results. One other possibility 
is the use of a pair of ‘‘long’’ streamlines (such as 37/8 
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Fic. 17. Pressure drop through nozzle. 


and 7/8), cut off at points nearer the origin. Tests on 
bends of angles other than 90° would be desirable. As 
stated previously, the ‘‘longer’’ bends give better 
results (7/8 and 37/8 being the ‘“‘longest’’), and elbows 
which fulfill the conditions of no inflection points; 
that is, Var. Z u, have proved to give the best re- 
sults. 

The results show conclusively that the very ‘‘short”’ 
special elbows have poor flow characteristics, as evi- 
denced by Figs. 8 and 9. 

In the case of the ‘‘shortest’’ elbow tested, the air at 
the inner wall experiences a reversal of flow which, 
with the accompanying turbulence, greatly increases 


the losses. With the ‘‘longer’’ elbows, where the air is 
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turned more gradually, lower losses would be expected, 
and, as stated before, this was found to be true. 

A complete test program should include all the 
variations mentioned in the foregoing discussion and, in 
addition, tests at low velocities would be desirable. 


CONCLUSIONS 


Based solely on the results of these tests, the following 
conclusions can be made: 

(1) The “short” elbows of special design have 
higher losses than conventional bends of the same size. 

(2) Lower losses are obtained with “‘longer’’ special 
elbows having carefully selected streamlines. 

(3) Increasing the aspect ratio of a bend from 0.67 
to 1.5 reduces the losses. 

(4) Addition of guide vanes to the “‘short’’ elbows, 
as done in these tests, does not appear to be worth while. 

(5) Still lower losses are indicated by proper com- 
bination of point of cut off of streamlines and aspect 
ratio. 

(6) Further investigation of bends having greater 
and lesser angles than 90° with various streamlines 
would appear to be in order. 
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The Effect of Radius of Curvature and Pre- 
liminary Artificial Eccentricities on Buckling 
Loads of Curved Thin Aluminum-Alloy 


Sheets for Monocoque Constructions 


GEORGES WELTER* 
Ecole Polytechnique de Montréal 


ABSTRACT 


In two previous papers, ' ? the results of tests on thin curved 
aluminum-alloy sheets in compression were presented. The 
present paper, whose publication was restricted during the war, 
shows that sheets of aluminum alloy tested under special condi- 
tions offer highly increased buckling and ultimate strengths 
under compression loads By decreasing the radius of cold- 
rolled curved sheets 40 to 50 per cent, a much greater ultimate 
buckling load can be obtained if these panels are tested under 
the usual radius of curvature (24 in.). Buckling and ultimate 
loads are distinctly higher, the usual buckling loads being ex- 
ceeded by about 100 per cent. 

Local eccentricities of the curved panels, generally assumed to 
be exceedingly detrimental to the bearing capacity in compres- 
sion, did not decrease the buckling load of the panels. Even 
holes produced by shooting bullets through the panels did not 
weaken materially their normal buckling strength or their ulti- 


mate resistance. 


INTRODUCTION 


IT A PREVIOUS INVESTIGATION, !: ? it was observed that 
if the radius of curvature of curved thin aluminum- 
alloy sheets is decreased from 24 in. to 10 in., and tested 
with a radius of 24 in. by opening the plate elastically 
to the greater radius, a maximum buckling load was 
reached. This preliminary observation was based only 
on an insufficient number of specimens, so that this 
point had to be verified and a few series of tests with 
panels having a radius of curvature of 6 and 8 in., re- 
spectively, increasing in steps of 2 to 4 in. up to 24 in., 
were made. Also, completely flat panels as well as 
panels with a negative curvature of minus 24 in. were 
loaded with a curvature of plus 24-in. radius. 


TESTS 


Two complete series of tests were carried out in order 
to obtain reliable results. The first series was con- 
ducted with 17ST sheet “KB” of 0.036-in. thickness 
(24-in. width and 96-in. length). Another series, com- 
prising tests on 18 panels, was made with cold-rolled 
pure aluminum sheet 2S */,H of 0.032-in. thickness, 
designated by “L.’”’ The width of the panels was, in 
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all cases, 6 in. and the length between the grips was 9 
in. 


Some results of tests on duralumin sheets are given in 
Fig. 1. The radius of curvature of these sheets varied 
between plus 24 in. and minus 24 in. The records of 
this series are so conclusive that it is hardly necessary 
to explain in detail these results. Similar curves are 
given in the report of the main investigation. More in- 
formation regarding the test setup, the test apparatus 
used, the uniformity of loading achieved, the degree of 
edge restraint imposed, and the precautions taken to 
assure uniform condition throughout the series of tests 
may be found in a previous article by the author.’ Of 
the first two panels with a radius of 8 in. and tested 
with 24-in. radius, one panel buckled at 3,080 Ibs. and 
the ultimate load was lower than the buckling load. 
The second panel (KB 4) showed no buckling and the 
stress-strain curve went up to the maximum load of 
3,850 Ibs. without any inflection before attaining the 
ultimate value. The same result was obtained with 
three panels of 12-in. radius (KB 2, 5, and 13): maxi- 
mum loads of 3,800 to 4,090 lbs. were reached without 


the characteristic buckling effect at lower loads. For 
higher radii, however, the buckling reappeared. Thus, 


for instance, with panel KB 9 of 14-in. radius, the buc- 
kling took place at 3,000 Ibs., while the maximum load 
of 3,420 lbs. was appreciably lower than with the 
panel of 12-in. radius. Furthermore, somewhat lower 
buckling loads were recorded for panels KB 16 and 21 
with radii of 18 and 24 in., respectively. The buckling 
took place at 2,510 Ibs. and 2,630 Ibs., while 3,500 Ibs. 
and 3,370 Ibs. were necessary to make these panels col- 
lapse. The flat plate KB 3 (radius ©) tested with a 
radius of 24 in. had a still smaller buckling load (1,650 
Ibs.); it showed, however, about the same maximum 
load (3,410 Ibs.) as the preceding panels. The lowest 
loads under which the panels buckled were registered 
with panels of negative curvature (minus 24-in. radius) ; 
panel KB 6 withstood only 960 Ibs. and panel KB 20 
only 940 Ibs. before buckling. This last panel was 
rolled first with a radius of minus 24 in.; then it was 
rolled in the opposite direction to a radius of plus 24 
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Fic. 1. Duralumin 17ST. Influence of radius of curvature. 
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Fic. 2a. Aluminum. Influence of radius of curvature. 
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Fic. 3. Influence of radius of curvature. 


in. and back again to minus 24-in. radius. The maxi- 
mum load of two panels treated this way was 3,270 and 
3,480 Ibs. 

It will be seen from these results that the radius of 
curvature of the panels being smaller than that of the 
steel grips with a radius of 24 in., has a pronounced 
influence on the buckling load as well as on the ultimate 
load of the sheets. According to these results, buckling 
does not take place with panels between about 8- and 
14-in. radius. No buckling can be observed and the 
collapse load corresponding to the buckling load is 
about 20 per cent higher than that of panels with 
smaller or greater radii. These maximum loads are those 
experienced in the main investigation’ (pages 368, 
369 and Fig. 19, reported as ‘‘maximum special’’). 

To confirm the observations of this interesting series, 
the same kind of tests was repeated with */, hard- 
rolled pure aluminum panels 2S of the same thickness 
(0.032 in.) and with curvatures of plus 6- to plus 24-in. 
radius. 

As shown in Figs. 2, this series of tests is also convinc- 
ing with regard to a certain increase of the maximum 
load for panels with smaller radii. Thus, for instance, 
panels with radii smaller than about 12 to 14 in. showed 
no buckling and collapse occurred at loads between 
1,580 and 1,782 Ibs. (L 1,9, 2, 10,11). Beginning with 
a radius of 12 in. the panels started buckling at loads as 
high as 1,566 Ibs. For a radius of 14 in., the buckling 
in one case occurred at the high load of 1,822 Ibs. 
(L12). From here the buckling load decreased syste- 
matically with increasing radius, reached 1,294 to 1,432 
lbs. for a radius of 24 in., and was as low as 884 to 946 
Ibs. for a radius of minus 24 in. 
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Buckling loads and ultimate loads of 17ST and 2S */,H sheets (0.032-in. thickness). 


It is interesting to note that, between 12- and about 
20-in. radius, the ultimate load at the point of collapse 
is appreciably lower than the buckling load: as much 
as several hundred pounds. However, at radii of 24 in. 
or more, the stress-strain records again became normal— 
that is, the ultimate load is, as in the usual tests, higher 
than the buckling load. 

In general, it may be said that these results corre- 
spond fairly well with those of the preceding series with 
regard to the influence of the radius of curvature on the 
behavior of these thin light-alloy panels. Furthermore, 
it is of interest to state that the buckling load for */, 
hard cold-rolled and stretched pure aluminum, with a 
yield point of about 35 per cent, that of 24ST alloy, 
has for panels with 12- and 14-in. radius a buckling load 
that is about as high as that of 24ST alloy. On the 
other hand, however, the collapse load is only about 
half that of 24ST, as may be seen in Fig. 2a, representing 
normal buckling test results with 24ST sheets of the 
same thickness (0.032 in.). Even for a radius of 24 in., 
the buckling load of 1,294 to 1,432 Ibs. is not much 
smaller than that of 24ST sheet, especially if we take into 
consideration the fairly wide variations measured dur- 
ing the main tests. In any case, it -an be seen from 
these investigations that something important happens 
with panels rolled to a radius between 8 and 14 in. and 
afterward opened for testing to a radius of 24in. The 
average results of both these series are illustrated in 
Fig. 3 which shows the pronounced effect of their radius 
of curvature on the test results. 

To investigate this region somewhat more precisely, 
it seems necessary to make a few more series of supple- 
mentary tests in this field. The absence of buckling 
for radii between 8 and 14 in. and the important in- 
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crease in ultimate load of these panels may be of prac- 
tical interest for highly stressed skin constructions. 


EFFECT OF PRELIMINARY ARTIFICIAL ECCENTRICITY 


In the bibliography on the theory of buckling of 
curved aluminum-alloy sheets in compression, the 
opinion is generally expressed that the maximum buc- 
kling load of panels is extremely sensitive to slight differ- 
ences in initial eccentricities of the sheets. The wide 
variations of the buckling loads for individual tests are 
mainly explained by these initial eccentricities. 

In order to study the influence of the initial eccen- 
tricity of sheets on the buckling load, a few tests on 
24ST curved panels (radius of curvature 24 in.), pro- 
vided with initial local eccentricities of different values 
were made. Thus, several indentations with a 10-mm. 
steel ball were made in the upper third of the free length 
of a panel. The pressures were of 500, 1,000, 1,500, 
and 2,000 kg. and the eccentricities of the indentations 
varied between about one to three times the thickness 
(0.032 in.) of the sheet. The buckling of this panel with 
nine indentations of different eccentricities, as shown 
in Fig. 4 (H 36), placed under compression load, 
was not influenced at all by these artificial eccen- 
tricities. The buckling took place in the lower part of 
the panel, as is usual with panels free of eccentricities, 
and the buckling had the normal appearance of that in 
panels tested without eccentricities. The load-deflec- 
tion diagram of this panel was also normal and the max- 
imum buckling load of 1,530 Ibs., as well as the ultimate 
load of 2,860 Ibs., was not different from other results 
obtained with this material under the usual test con- 
ditions. 

However, in order to determine if indentations with 
a steel ball of larger diameter would not give somewhat 
different results, this test was repeated with panel Ic 7. 
Two indentations with a ball of 1-in. diameter under 
loads of 500 kg. were made in the upper and the lower 
parts of the panel, while four other indentations with a 
loadof 1,000 kg.and the same ball were made on the right 
and left sides in the upper and lower parts of this panel 
(Fig. 4). Again, the compression test showed that 
buckling was not at all influenced by this series of arti- 
ficial eccentricities of a few times the thickness of the 
sheet. On the contrary, the direction in which buc- 
kling took place seems to indicate that it avoids these 
eccentricities. The buckling appeared right in the 
middle of the sheet — that is, in that part of the sheet 
which was free of any eccentricity. The lower edge of 
the buckling seems even to stop along and at some dis- 
tance from the main direction of these indentations. 

Therefore, it may be said that artificial eccentrici- 
ties of this kind are not the cause of premature buc- 
kling, and the divergence of buckling results of thin 
sheets under compression loads cannot be attributed 
mainly to this cause. Moreover, the load deflection 
diagram does not show any abnormal shape and the 
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Fic. 4. Influence of eccentricity. 


buckling load of 2,060 Ibs. is comparatively high for 
this sheet (see Fig. 4). 


As these tests had given ithe desired information 
about local artificial eccentricities, complementary 
tests with a somewhat different kind of eccentricity 
were executed. As shown in Fig. 5, three indentations 
of maximum effect were made in panel H 26. With 
a 0.22-in. rifle, three bullets were shot through this 
panel (two in the upper part and one in the lower part 
of the panel, Fig. 4, H 26). The edge of the holes was 
distinctly deformed plastically, so that part of the mass 
of the sheet was forced outside the plane of the panel. 
Despite this rather crude lccal deformation of the sheet, 
buckling did not start from or spread across these initial 
eccentricities. Independently of these holes produced 
dynamically, buckling took place in the usual manner, 
not in the least influenced by this kind of deformation. 
Consequently, the stress-strain record showed no peculli- 
arity whatever and it can be characterized as normal 
for this sheet. 


It may be seen from these tests that, even with a 
pronounced artificial eccentricity produced by ball 
indentations of 10-mm. and up to 1-in. diameter, 
which under loads of up to 1,000 kg. result in eccen- 
tricities of several times the thickness of the sheets, 
there could not be observed any kind of influence on 
the mechanism of the buckling of these thin alloy sheets. 
If the sensitivity of the sheets to slightest differences 
in initial eccentricities is as great as supposed by 
different investigators, it would have been proved by 


(Continued on page 604) 
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Higher Modes of Vibration by a 
Method of Sweeping 


L. BESKIN* ano R. M. ROSENBERG{ 


Consolidated Vultee 


SUMMARY 


A method is presented whereby the frequencies and deflection 
curves in bending and torsional vibrations of continuous beams 
of variable section can be found. Emphasis is placed on the de- 
termination of modes higher than the first. 

This method was developed in an effort to use the Stodola iter- 
ation process in the determination of higher modes of vibration. 
One of the most important advantages of the method presented 
here is that its application does not require unfamiliar tools of 
mathematics such as matrix techniques. Therefore, actual cal- 
culations may be reduced to a routine that may be entrusted to 
computers having a high-school background only. 

This paper, containing a numerical example worked out in de- 
tail, was written in such a manner as to emphasize application 
rather than veil it. 

In applying the method, it is necessary to determine the first 
(n — 1) modes. The mth mode is then found in the following 
manner: 

(1) A “reasonable’’ estimate (as defined in the body of the pa- 

per) of the mth mode is made. 

(2) The preceding (w — 1) modes are cleared from the esti- 

mated mth mode by a method of sweeping. 

(3) The swept mode is operated on with Stodola’s iteration 

method leading to convergence to the mth mode. 

For great accuracy, alternate sweeping operations and Stodola 
iterations are continued leading to an indefinitely close approxi- 
mation of the mth mode. 

This method has been used by the authors with excellent results 
and presents the best compromise between simplicity, accuracy, 
and reduction of labor known to them. 


INTRODUCTION 


I THE FIELDS of aeronautical and mechanical engi- 
neering the need frequently exists for computing the 
frequencies and deflection curves of higher modes of 
vibration of continuous systems. In systems of con- 
stant section, straightforward solutions to this problem 
exist. In systems of variable section, more or less labo- 
rious iteration methods are known which permit the 
calculation of approximations to these higher modes. 
Usually, these iteration methods rely on either of two 
principles: (1) The Holzer method and its variations 
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make use of the fact that, in the absence of damping, 
the force in bending vibrations (or the torque in tor- 
sional vibrations) required to maintain vibrations in 
the principal modes, is zero. (2) Most other methods 
make use of the orthogonality relation between normal 
modes. In application, the Holzer method and its 
variations require an estimate of the frequency of vi- 
brations of the mode; the methods using orthogonality 
restraints require an estimate of the deflection curve of 
the mode. It is believed that a good estimate of the 
deflection curve of the mode requires far less experience 
and can be much more easily visualized than a good 
guess of the frequency. In this respect, methods relying 
on the orthogonality principle are, therefore, more eas- 
ily applied by novice and expert alike. However, those 
of the latter methods now in general use suffer from 
certain defects. Some, although theoretically sound, 
give numerically unsatisfactory results; others require 
a knowledge of matrix methods to carry out the itera- 
tions. Therefore, the present method was developed 
in an effort to provide a technique that would be so sim- 
ple to use that a high-school graduate could apply it 
without learning special mathematical techniques, 
that would give excellent numerical results carrying 
only a small number of significant figures (usually five), 
and that would be sufficiently rapid to be no more te- 
dious than methods now in general use. 

A relatively large number of contributions!: ?: * to 
the field of vibration analysis have recently appeared 
in the literature. These have, generally, provided 
worth-while labor saving devices or greater accuracy of 
results. They have, however, generally relied on the 
principles enumerated above for the calculation of 
higher modes of vibration. 

In this respect the present paper is no exception. 
A method of ‘‘sweeping”’ is applied to an assumed deflec- 
tion curve in order to clear from it the lower modes. 
Stodola’s method‘ is then applied to the swept curve 
and results in the convergence of the Stodola process to 
the nth mode if the preceding (x — 1) modes have been 
swept from the function. 


In applying the method discussed here to the mth 
mode a “reasonable’’ estimate of the deflection curve 
of the nth mode is made. A ‘‘reasonable’’ estimate dif- 
fers from a “‘possible’’ assumption! in that it not only 
satisfies the boundary conditions, as does the latter, but 
it also possesses the known general character of the 
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mode. Thus, a reasonable deflection curve of the nth 
mode of a cantilever beam in bending vibrations pos- 
sesses (n — 1) nodes and has zero deflection and slope at 
the root. The components of the preceding (n — 1) 
modes (as determined by previous computations) are 
cleared or ‘‘swept’’ from the assumed curve, whereupon 
it is subjected to the Stodola process. 

If the preceding (m — 1) modes were exactly known, 

they could be completely swept out, and, if the nu- 
merical computations were made with complete accu- 
racy, repeated operation on the swept curve with the 
Stodola process would lead to continued convergence 
to the mth mode. This procedure is, unfortunately, 
not possible because all modes up to (m — 1)st were, 
themselves, determined by a converging iteration proc- 
ess. They can, therefore, be swept from the assumed 
curve only within the accuracy with which they have 
been determined. It must be assumed, then, that re- 
siduals of the lower modes are left even after sweeping, 
and continued application of the Stodola process will 
amplify these residuals to the extent that repeated 
Stodola iterations, done upon the swept curve, will 
eventually converge to the first mode. This difficulty 
can be overcome if each Stodola iteration is followed by 
a sweeping procedure. If this is done, the method pre- 
sented here will continue indefinitely to converge to 
the desired mode. It can be shown that the sweeping 
operation need not be applied in its entirety to each 
iteration. For instance, when calculating the fourth 
mode, if the preceding three modes have been deter- 
mined with equal accuracy, the first mode should be 
swept out after each iteration; but it is sufficient to 
sweep out the second mode after every other iteration 
and the third mode after every fifth. Practice has 
shown that the frequency can be calculated within less 
than 1 per cent of that resulting from three iterations 
and four complete sweeping operations if the assumed 
curve is swept once, and the frequency is calculated 
from the results of a single Stodola iteration. This 
minimum of work is sufficient in problems involving 
forced vibrations and resonance. In three-dimensional 
flutter analysis, however, accurate knowledge of the 
deflection curve is required, and the location of the 
nodes may determine stability or instability against 
flutter. In these cases it is recommended that the suc- 
cession of sweeping operations and Stodola iterations 
be continued until successive values for the deflection 
curve and frequency are obtained which do not change 
sensibly after additional iterations. 

The nomenclature used here is given in Table 1. 


SToDOLA’sS ITERATION METHOD 


Stodola’s iteration method may be described as fol- 
lows: A system of length /, having a mass per unit length 
dm, is assumed to deflect in a shape 7,'(x), henceforth 
identified as j,'. If the frequency of the system isw = 
1, the inertia loading per unit span is ¥,! dm. Loading 
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the system with j,' dm, a new deflection curve 7, is 
obtained. The ratio 7;'/j,? is a first approximation of 
the square of the frequency. This process may be re- 
peated by loading the system with ¥,? dm and will con- 
verge to the first mode as shown below.* 

The assumed deflection curve may be developed in 
the orthogonal series 


J} = Y, + a2¥2+...4,Y,+... (A) 


When the frequency » = 1, the inertia loading be- 
comes 


Jil dm = Y,dm + a.Y2dm+...a,Y,dm+... (B) 


For the deflection curve of the nth mode, the inertia 
loading is (w,)*Y, dm. Therefore, the loading a, Y, dm 
gives a deflection a, Y,/(w,)?, so that the second de- 
flection curve of the process becomes 


a Y a2V2 a,f, 


es ab 
” (@)? ” *** (@,)? 


(w1)? 


This differs from the first (assumed) deflection curve in 
that each term is divided by the square of its frequency. 
The gth deflection curve obtained from the (¢ — 1)st 
iteration is 


1 2¢ 2¢ 
eon) oats + (3) att] 
1 2 3 


(D) 





+... (C) 


Since it can be proved that w; < w: and w;< w,, etc., all 
terms higher than the first decrease with increasing q. 
When g— = all terms higher than the first disappear. 


DETERMINATION OF THE SECOND MODE 


Let an assumed deflection curve of the second mode 
be represented by the orthogonal series 


ja! = Yi +a2¥2+...4,Y,+... (1) 
Inasmuch as the orthogonality relation 
JS'Y,,Y, dm =0,m+n (2) 


* This proof is given in reference 5. It is repeated here be- 
cause the same proof will be shown later to apply to convergence 
of the described process to any mode. 
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HIGHER MODES 


holds for principal modes, it can be shown that the co- 
efficient of the mth term is 


a, = Soin Y, dm/ So'( Y,)? dm (3) 
Substituting Eq. (3) for m = 1 into Eq. (1) results in 


So'52! Vi dm 


BS a: Sie ’ 
Je Se(¥a? dm ties 


Yi+a¥2+. 
The first mode is now cleared from the assumed curve 
and the swept function becomes 


'5.1V, dm 
yo! = Je} — ee Y, = a2V2 +.. ‘7, ee 
(5) 
The swept curve, obtained by taking the difference as 
indicated in the middle portion of Eq. (5), is represented 
by the series of the right-hand side which contains no 
first mode impurities. In carrying out the sweeping 
process numercially, the final deflection curve, obtained 
by applying the conventional Stodola process to the 
system, is used for Yj. 

If the series in Eq. (5) ts divided through by a, it 
evidently becomes similar to the series of Eq. (A). 
Therefore the proof of convergence to the first mode of 
Eq. (A) holds identically for Eq. (5), except that re- 
peated application of Stodola’s procedure to yz causes 
convergence to the second mode. 

One of the numerical limitations of the process is 
that the accuracy with which a higher mode can be de- 
termined is a function of the errors made in the previ- 
ous modes. 

Let 


Y, = w'+ ANY 
where y,° is the deflection curve of the first mode ob- 


tained after (¢ — 1) Stodola iterations. The quantity 
swept from the curve is then, 


o = LOWY: — AY:) dee 
aac SON: — AY,)* dm 
whereas, for complete sweeping, a quantity 


ay, = (fo'yt¥:1 dm/ So'( V1)? dm) Y; 


(Y, — AY,) 
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TABLE 2 
Beam Parameters and Estimated Reasonable Deflection Curves 








tion* It wt ai} Fo} Fs! av 
10 0.0001 0.0020 1.000 1.000 1.000 1.000 
9 0.0005 0.0235 0.780 0.485 0.262 0.300 


8 0.0014 0.0385 0.591 0.100 —0.275 —0.195 
7 0.0026 0.0550 0.4832 -—0.125 —0.475 —0.22% 
6 0.0044 0.0690 0.309 —0.270 —0.312 —0.066 


5 0.00725 0.0843 0.202 —0.320 —0.020 0.140 

4 0.0105 0.0990 0.135 —0.265 0.226 0.060 

3 0.0148 0.1150 0.080 —0.165 0.275 —0.136. 
2 0.0225 0.1262 0.038 —0.090 0.160 —0.140 

1 0.0246 0.1342 0.011 —0.030 0.050 —0.060 

0 0.0246 0.1342 0 0 0 0 











 * Ax = 3.15 in, 
T (In.)*. 
t Lbs. per in. 











TABLE 3 
Computed Deflection Curves—First Mode 
Station ji} ji? ar* j4 
10 1.00000 1.00000 1.00000 1.00000 
9 0.78000 0.83895 0.84194 0.84216 
8 0.59100 0.68043 0.68600 0.68642 
7 0.43200 0.53199 0.53907 0.53962 
6 0.30900 0.39946 0.40680 0.40737 
5 0.20200 0.28594 0.29251 0.29303 
4 0.13500 0.19163 0.19683 0.19724 
3 0.08000 0.11693 0.12051 0.12080 
2 0.03800 0.06121 0.06324 0.06341 
1 0.01100 0.02155 0.02232 0.02238 
0 0 0 0 0 
w;? = 1,392 cycles per min. Z(ji* — Fi)? = 0.048827654 
wi} = 1,392 cycles per min. =(7i4 — F172)? = 0.000269215 
wi = 1,391.5 cycles per min. =(ji4 — 7,5)? = 0.000001407 


should have been eliminated. Therefore, a first mode 
residual 
 - LOY dm S'5i(Va — AV:) dm 
’ So(¥1)? dm . So(¥i1 — AY:)? dm 
(Yi-— AY) (6) 





is still contained in the function y2', even after sweeping. 
It can be shown that the major portion of the residual is 


a = (Yi Jon ay dm— 
AY, Joi" V1 dm)/ So'(¥1)? dm 
Operation on 2 with the Stodola process will amplify 


this residual because, as was shown in the section ‘‘Sto- 
dola’s Iteration Method,’ the Stodola process always 




















TABLE 4 
Computed Deflection Curves—Second Mode 
Station Fe! y2} 52? 2? 528 y2* 524 y24 
10 1.000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 

9 0.485 0.53705 0.60459 0.60465 0.60974 0.60986 0.61060 0.61068 
8 0.100 0.18545 0.24010 0.24022 0.24645 0.24667 0.24775 0.24790 
7 —0.125 —0.02815 —0.03262 —0.03247 —0.02919 —0.02891 —0.02823 —0.02804 
6 —0.270 —0.17129 —0.19123 —0.19108 —0.19162 —0.19132 —0.19151 —0.19132 
5 —0.320 —0.23067 —0. 24690 —0.24676 —0.25090 —0.25062 —0.25165 —0.25147 
4 —0.265 —0.19764 —0.23114 —0.23103 —0.23783 —0.23761 —0.23909 —0.23894 
3 —0.165 —0.12335 —0.17421 —0.17413 —0.18143 —0.18128 —0.18271 —0.18261 
2 —0.090 —0.06764 —0.10405 —0.10400 —0.10939 —0.10931 —0.11031 —0.11025 
1 —0.030 —0.02237 —0.04058 —0.04056 —0.04302 —0.04299 —0.04342 —0.04340 
0 0 0 0 0 0 0 0 0 

we? = 4,940 cycles per min. Z(y.4 — yo')? = 0.017630543 

w2? = 4,936 cycles per min. Z(y.4 — y2)? = 0.0003188165 

wet = 4,936 cycles per min. X(y.4 — y.3)? = 0.0000082542 





ee 
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TABLE 5 et 
Computed Deflection Curves—Third Mode 
Station Fs! ys! Fs? ys? 938 ¥s8 3" y3* = 
10 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 —- 
9 0.26200 0.38104 0.38181 0.38281 0.38063 0.38131 0.38042 0.38107 ~ 
8 —0.27500 —0.08438 —0.10265 —0.10105 —0.10367 —0.10260 —0. 10364 —0. 10263 Biaces 
7 —0.47500 —0.28615 —0.27778 —0.27640 —0.27455 —0.27371 —0.27403 —0.27327 
6 —0.31200 —0.20207 —0.18728 —0.18673 —0. 184386 —0.18419 —0.18416 —0.18404 
5 —0.02000 —0.00648 —0.00101 —0.00135 —0.00205 —0.00253 —0.00243 —0.00296 mer, 
4 0.22600 0.16601 0.14197 0.14109 0.13847 0.13760 0.13806 0.13717 9 
3 0.27500 0.19927 0.18329 0.18234 0.18230 0.18143 0.18244 0.18157 8 
2 0.16000 0.11455 0.13867 0.13799 0.14128 0.14066 0.14181 0.14119 } 
1 0.05000 0.03505 0.06252 0.06222 0.06521 0.06495 0.06563 0.06536 5 
0 0 0 0 0 0 0 0 0 : 
w3? = 11,170 cycles per min. Z(yst — ys')? = 0.003619254 ; 
w3? = 11,161 cycles per min. =(y3* — ys”)? = 0.000061208 0 
w3* = 11,161 cycles per min. =(y3* — ys)? = 0.000001113 — 
ze = - < 
TABLE 6 4 
Computed Deflection Curves—Fourth Mode | 
Station ia ys! je? ys" 533 3 ja4 ys* 
10 1.00000 1.00000 1.00000 1. 00000 1.00000 1.00000 1.00000 1.00000 an 
9 0.30000 0.22292 0.18475 0.18538 0.16756 0.16841 0.16146 0.16033 | 
8 —0.19500 —0.28717 —0.28106 —0.28027 —0.27813 —0.27707 —0.27312 —0.27486 | 
7 —0.22800 —0.20652 —0.17033 —0.17003 —0. 14676 —0.14632 —0.13365 —0.138516 — 
6 —0.06600 0.056380 0.12258 0.12242 0.13047 0.138035 0.13202 0.13114 
5 0.14000 0.32571 0.22162 0.22169 0.19094 0.19084 0.17946 0.17907 ai5% 
4 0.06000 0.13897 0.06748 0.06779 0.04866 0.04895 0.04204 0.04189 Sta 
3 —0. 13600 —0.20880 —(.12747 —0.12679 —0.11068 —0.11005 —0.10447 —0.10456 = 
2 —0.14000 —0.22881 —0.17457 —0.17392 —0.15109 —0.15050 —0.14028 —0.14034 
1 — 0.06000 —0.10137 —0.10148 —0.10115 —0.08973 —0.08942 —0.08299 —0.08301 
0 0 0 0 0 0 0 0 0 : 
ggg 19,973 cycles per min. Z(yet — ya!)? = 0.064720144 10 
ws? = 19,749 cycles per min. Z(ys4 — 947)? = 0.0063868492 9 
ws? = 19,712 cycles per min. Z(yit — 943)? = 0.0005581705 P 
a Ta a a ee ea el 5 ies sistahopade manera? eee é 
Computed Deflection Curves—Fourth Mode, with Partial Sweeping 5 
Station 5a) yt ie" y2* 53* °° j** ys! 4 
10 1.00000 1.00000 1.00000 1.00000 ‘1.00000 1.00000 1.00000 1.00000 3 
9 0.30000 0.22292 0.18475 0.18522 0.17015 0.16717 0.15484 0.15613 
8 —0.19500 —0.28171 —0.28106 —0.28036 —0.27442 —0.27811 —0.27995 —0.27812 . 
7 —0. 22800 —0.20652 —0.17033 —0. 16982 —0. 14432 —0.14558 —0. 13380 —0.13253 1 
6 —0.06600 0.05630 0.12258 0.12279 0.13071 0.13218 0.13660 0.13711 
5 0.14000 0.32571 0.22161 0.22187 0.18970 0.19197 0.18231 0.18252 0 
4 0.06000 0.138897 0.06748 0.06757 0.04705 0.04844 0.03978 0.04007 
3 —0.13600 —0.20880 —0.12747 —0.12729 —0.11193 =—0.11177 —0.110438 —0.10999 - 
2 —0.14000 —0.22881 —0.17457 —0.17439 —0.15178 —0.15220 —0.14608 —0. 14568 
1 —0.06000 —0.10137 —0.10148 —0.10139 —0.08997 —0.09032 —0.08605 —0.08585 
0 0 0 0 0 0 0 0 0 - ae 
wy? = 19,973 cycles per min. X(yst — ys)? = 0.0647201444 
w4?* = 19,748 cycles per min. D(yat — yu2*)? = 0.0064259824 ws? = 
ws** = 19,709 cycles per min. Z(ya4 — yu*)? = 0.0006223984 
wy = 19,712 cycles per min, X(t — yst*)? = 0.0001521064 
bra 
amplifies the lowest mode present in the series. There- In the general case, the swept function of the nth mode tior 
fore, unless the first mode has been computed to high after the (¢ — 1)st Stodola iteration is given by the 
accuracy, it will be necessary to sweep after each ap- - wit! 
plication of Stodola’s procedure leading, for the second ; " So5,°Y, dm Vv (8) sect 
7 Re Fue “IW \2 dm of 7) 
mode, from 2% to y2". _ si — So(i 2 dm rect 
isa 
DETERMINATION OF Mopes HIGHER THAN THE SECOND sideé 
; PRACTICAL APPLICATION 
From the preceding development it is clear that, Th ed f lati ; 
when determining the third mode, the first and second 4 - eae a. a _ ning vibration I 
; 7 g ass ¢ astic chi - ) 
modes are swept from the deflection curve by the oper- wre es of a system are the mass and elastic characteris 
dian tics along the span. After these are known it is con- con 
; ' venient to divide the length into ten or more equal ope 
TT <1V i j - : 
=H) — S05s' Vi dm + So5s Y2dm (7) increments, called sections. Stations are located in the Alt! 
rs Ne . a . . : 
So(V1)? dm So( V2)? dm center of each section. When calculating bending vi- two 
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TABLE 8 
Third Mode—First Sweeping Operation 
@ D ® @ ® ® ® ® @ D i 
Sta- = ere 
tion ¥i* Y.* Sat w w(¥i)2 w(¥:2)?2 wY¥i5s! wY¥253! a'F; a2! ¥2 ya as (normal) 
- (—0.17754) (—0.16183) } g 
@®@x@ @©x@’@x@x@@x@x® x@®@ x@ @®-@-@ @/1.33937 
o 2 1.00000 1. 0.00200 0.0020000 0.0020000  0.0020000  0.0020000 —0.17754 —0.16183 1.33937 1.00000 
9 0.84216 0.61068 0.26200 0.02350 0.0166672 0.0087639 0.0051852  0.0037600 —0.14952 -—0.09883 0.51035 0.38104 
8 0.68642 0.24790 —0.27500 0.03850 0.0181400 0.0023660 -—0.0072674 —0.0026246 -—0.12187 -—0.04012 -—0.11301 —0.08438 
7 0.53962 —0.02804 —0.47500 0.05500 0.0160154 0. 432 —0.0140975 —0.0007325 -—0.09580 0.00454 —0.38374 —0.28651 
6 0.40737 —0.19132 —0.31200 0.06900 0.0114508 0.0025256 -—0.0087700 0.0041187 -—0.07232 0.03096 —0.27064 —0.20272 
5 0.29303 —0.25147 —0.02000 0.08430 0.0072384 0.0053309 —0.0004940  0.0004240 -—0.05202 0.04070 -—0.00868 —0.00648 
4 0.19724 —0.23894 0.22600 0.09900 0.0038515 0.0056521 0.0044131 —0.0053460 -—0.03502 0.03867 0.22285 0.16601 
3 0.12080 —0.18261 0.27500 0.11500 0.0016782 0.0038348 0.0038203 -—0.0057750 -—0.02145 0.02955 0.26690 0.19927 
2 0.06341 —0.11025 0.16000 0.12620 0.0005074 0.0015340  0.0012803 —0.0022262 -—0.01126 0.01784 0.15342 0.11455 
1 0.02238 —0.04340 0.05000 0.13420 0.0000672 0.9002528  0.0001502 —0.0002912 -—0.00397 0.00702 0.04695 0.03505 
0 0 0 0 0. 13420 0 0 0 0 0 0 0 0 
4.17243 0.812: 0.89100 0.076161 0.0323033 —0.0137798 0.052278 = 0.740770. 13150 1.76327 1.31648 
ai! = Sw¥151/ zw (¥)? = —0.0137798/0.0776161 = —0.177538 ie we a 
a2! = Sw¥25!/Zw(¥2)? = —0.0052278/0.0323033 = —0.161835 
* Obtained from preceding calculations. 


t+ Reasonable estimate of third mode. 









































TABLE 9 
Third Mode—First Stodola Iteration 
i 2 (3) 4) 6) (6) G 8) 9) 40 @ 13 
Sta- 1 dy (os 3) 1G 3) (432) 2/ w(43*)? 
tion I w ya! wys! S M M/I 3.15 dx 3?/9.922 9.922 / norm. 98. 4 56 98.456 
Table : 2 T able 2 ‘Table 8 
Column Column Column ; bi \ . . ‘ 
1 2 12 2) x @) {= {= Ee /D t=@ f= 9) /29.91263 (6) x @ 9)? 2) x @ 
‘ 10 0 00010 0.00200 1.00000 0.002000 29 .91263 1.00000 894.7654 1.78953 
0.002000 18.49157 
i) 0.00050 0.02350 0.38104 0.008954 0.002000 4.00000 11.42106 0.38181 0.008000 130.4406 3.06535 
0.010954 14.49157 
8 0.00140 0.03850 —0.08438 —0.003249 0.012954 9.25286 —3.07051 —0.10265 0.119862 9.4280 0.36298 
0.007705 5. 23871 
7 0.00260 0.05500 —0.28651 —0.015758 0.020659 7.94577 —8.30922 —0.27778 0.164152 69.0431 3.79737 
— 0.008053 — 2.70706 
6 0.00440 0.06900 —0.20207 —0.013943 0.012606 2.86500 —5.60216 —0.18728 0.036116 31.3842 2.16551 
—0.021996 — 5.57206 
5 0.00725 0.08430 —0.00648 —0.000546 —0.009390 —1.29517 —0.03010 —0.00101 0.012162 0.0009 0.00008 
—0.022542 — 4.27689 
4 0.01050 0.09900 0.16601 0.016435 —0.031932 —3.04114 4.24679 0.14197 0.097110 18.0352 1.78548 
i — 0.006107 —1.23575 
3 0.01480 0.11500 0.19927 0.022916 —0.038039 —2.57020 5.48254 0.18329 0.097768 30.0582 3.45669 
0.016809 33445 
2 0.02250 0.12620 0.11455 0.014456 —0.021230 —0.94356 4.14809 0.13867 0.020032 17.2067 2.17149 
0.031265 2.2780 
1 0.02460 0.13420 0.03505 0.004704 0.010035 0.40793 1.87008 0.06252 0.004094 3.4972 0.46932 
0.035969 1.87008 
0 0.02460 0.13420 0 0 0.046004 1.87008 0 0 0 086031 0 0 
0 64532 19 .06380 
60 eb Tn I(d2y/dx*)*dx 
w©=>5 — cycles per min., where g = 386.4 in./sec.2? and E = 10.3 X 10 Ibs. per sq.in. (here). 
= Se wy? is 
e 60 386. 4X 10.3 X 108 X 0.645327 Fs F 
oo? = — ——— = 11,170 cycles per min. 


6.2832 19.06380 X 98.456 


brations, the mass and elastic parameters of each sta- 
tion are paired. When calculating torsional vibrations, 
the mass moment of inertia at each station is patred 
with the torsional elasticity parameter at the end of the 
section in which the station is located. When these 
recommendations are followed, the length increment 
is a constant all along the span and may be taken out- 
side the summation sign—1.e., 


S f(x) dx = Yf(x) Av = AxD f(x) 


In performing the computations it has been found 
convenient to carry out Stodola iterations and sweeping 
operations on separate tables (see Tables 8 and 9). 
Although this requires working back and forth between 
two tables, it has the advantage of greater simplicity 


because each table contains repetitions of the same proc- 
ess. For the Stodola iterations a tabular method is 
preferred to the graphical one, originally devised by 
Stodola.* The method presented here is due to Burgess’ 
except that, because of a simplification, the number of 
operations has been reduced by one (see Table 9). 

In practice it has been found convenient to ‘‘normal- 
ze’ the deflection curves after each Stodola iteration 
and after each sweeping operation. Although this 
leads to two additional operations that are not essential, 
it is preferred because errors are more quickly discov- 
ered. The term “normalizing condition” is conven- 
tionally defined by /,'Y, dm = M. However, re- 
cently it has come into practice to refer to a mode as 
normalized if a reference amplitude (say the tip ampli- 
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tude of a cantilever beam) is equal to unity.® In this 
paper, the term “normalizing” is used in the latter 
sense—i.e., a mode is normalized by dividing the am- 
plitudes at all stations by the reference ampli- 
tude. 

In calculating frequencies from the deflection curves, 
Raleigh’s principle of equating maximum kinetic and 
strain energies is applied. This can be shown to be iden- 
tical to the weighted average method used by Boukidis 
and Ruggiero.! 


NUMERICAL EXAMPLE 


For a numerical example, the first four bending vi- 
bration modes of a cantilever beam have been com- 
puted by the method presented here. In all modes a 
reasonable deflection curve has been assumed and three 
Stodola iterations have been carried out. In all modes 
higher than the fundamental, the first Stodola itera- 
tion was preceded by a complete sweeping of all modes 
lower than the one to be calculated, and each Stodola 
iteration was again followed by a complete sweeping 
procedure. In addition to these computations, certain 
partial sweeping operations have been carried out for 
the fourth mode; these are compared with results 
from the complete procedure. For each mode the re- 
sults of each Stodola iteration have been used to com- 
pute frequencies identified by superscripts. Therefore, 
for the mth mode three frequencies have been calculated, 
identified as w,’, w,*, and w,4. The last is, of course, 
the most accurate value, and the others may be com- 
pared with the last to estimate the degree of accuracy 
obtainable with an abbreviated procedure. 

In addition to the frequencies, the deflection curves 
have also been compared. In the comparison of the 
deflection curves the principle of least squares has been 
applied. 

Let Ay, = Y, — y,% be the deviation of the computed 
values from the correct values. If the problem at hand 
were one of comparing observed values of a phenomenon 
with computed values, the Gaussian law of errors would 
give the probability of obtaining the observed values, 
as® 


P dai (h/ V a)me~ PCAyn)* 


9 


which isa maximum when }> (Ay,)? is a minimum. 
Therefore }>(Y, — y,*)? may be considered to furnish 
an evaluation of the accuracy of the computed deflection 
curves. In computing this quantity, the final computed 
mode has been defined as the pure mode—i.e., 


Y, =n 


and the quantities to be compared with one another 
are 


2\%a* — %n*)* 
Do (y_4 — 957)? 
200! — ¥a")* 
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.The results of this investigation are contained in Tables 
2 to 7. In addition, the first sweeping operation and 
Stodola iteration, used to calculate the third mode, are 
shown in detail in Tables 8 and 9 in order to illustrate 
the procedure. 

The fourth mode has also been calculated by the 
following abbreviated procedure: 


(a) Assume a reasonable deflection curve 7,'. 

(b) Sweep out the first, second, and third modes 
leading to y4!. 

(c) Apply the Stodola procedure leading to ¥,’. 

(d) Sweep out the first mode only leading to y4°*. 

(e) Apply the Stodola procedure leading to ¥,**. 

(f) Sweep out the first and second modes only leading 
to y4**. 

(g) Apply the Stodola procedure leading to ¥,**. 

(h) Sweep out the first mode only leading to y,**. 


This shortened procedure leads to excellent results 
for both the frequency and deflection curve, as shown in 
Table 7. 

Figs. 1, 2, 3, and 4 show the assumed and final de- 
flection curves for the first for bending modes. 
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VIBRATION DEFLECTION CURVES 
THIRD BENDING MODE 
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Appendix I 


Boukidis and Ruggiero! have presented a method of 
determining higher modes of vibration which is, in some 
respects, similar to the present one. From the stand- 
point of mathematical elegance their method is superior 
to the one contained in this paper, because all modes 
are calculated from a single assumed curve and because 
the tedious Stodola process is applied to the first mode 
only. In actual application, however, the method of 
reference 1 seems to have certain defects that render the 
method presented here preferable when modes higher 
than the second have to be calculated. 

It is shown in reference 1 that* 


Jot = Hi? — (w1)?H0*? (I-1) 


In the section on Stodola’s method of this paper, it was 
shown that the first term on the right-hand side of Eq. 


(I-1) is 


1 2q wy \74 
2 os E + (2) a2Y2+ (2) a3Ys+.. | 
(w) We w3 
(I-2) 
The second term on the right-hand side of Eq. (I-1) is 


1 2¢+2 
wit! = al y+ (=) a2Y2 + 
1 2 
Voge 


wy 2¢+2 
— as3V3 + o° 
W3 


Neglecting terms of an order higher than the second, 
the quantity of Eq. (I-1) is evidently 


1 1 29 1 2 = 
m= calla) L(G) ont 


Comparing Eq. (I-4) with either of Eqs. (I-2) or (I-3), 
it appears that yo" is the difference between two quan- 
tities that tend toward one another for large values of g 
(a large number of iterations). For proof we need only 
realize that in Eq. (I-4) 





* For reasons of simplicity, the bars will be omitted hereafter. 


(I-5) 


I 
o 


lim (w /wo) 2q 
I> © 
or, comparing the difference of Eq. (I-1) with either one 
of the subtrahends 


lim [Eq. (I-2) — Eq. (I-3)]/Eq. (I-2) = 0 (1-6) 
I o 
This fact explains the recommendation in reference 1 
that as many places should be retained in the calcula- 
tions as the calculating machine can handle. 

It will now be shown that the method of reference 1 
requires a relatively large number of iterations (large q) 
for dependable results. 

Since the second mode is calculated from the second 
term of the series of a function that is made to resemble 
the first mode as closely as can be estimated,' the coef- 
ficient of the second term (or the presence of the second 
mode in the function) will be small. Furthermore, no 
control can be exercised over the magnitude of the co- 
efficients of the various terms. Therefore, the coef- 
ficient of the mth term in the assumed curve may easily 
be larger than that of the mth term (where ” > m). 

The method presented by Boukidis and Ruggiero! 
has been used to calculate the second and third modes of 
the beam treated in the numerical example. To carry 
out these calculations, ten significant figures and, oc- 
casionally, 13 places were retained in the computations. 
This was necessary because retaining five significant 
figures as used in the method presented in this paper 
resulted in only a single significant figure for the third 
mode. 

The deflection curves, calculated by means of refer- 
ence 1, are identified by stars. The results are given 
in Table 10, and the average error is compared with 
that resulting from abbreviated calculations by means 
of the method presented in this paper. It is seen that 
the method of reference 1 gives acceptable results for 
the second mode (comparable to two to three Stodola 
and two to three sweeping operations), whereas results 
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for the third mode are poor. In fact, the average error 
resulting from application of the method of reference 1 
gives in the third mode an average error approximately 
50 times that which would result from using the esti- 
mated and swept curve of the present method but with 
no Stodola iterations. For this reason, fourth mode cal- 
culations by means of the method of reference 1 have 
not been attempted. 
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TABLE 10 





Comparison Between Present Method and That of Reference 1 
Station y24 y2* y3* 43* 
10 1.00000 1.00000 1.00000 1.00000 
g 0.61068 0.60918 0.38107 0.31362 
8 0.24790 0.24570 —0.10263 —0O.12084 
7 —0.02804 -—0.02947 -—0.273827 —0.22531 
6 —0.19132 -—0.19117 —0.18404 —0.12592 
5 —0.25147 —0.24991 —0.00296 0.00377 
4 —0.23894 —0.23668 0.13717 0.11828 
3 —0.18261 —0.18047 0.18157 0.14297 
2 —0.11025 —0.10879 0.14119 0.09722 
1 —0.04340 —0.04278 0.06536 0.03833 
0 0 0 UO 0 
=(y:4 — y2*)? = 0.0000237942 =(y:4 — ys*)? = 0.151152694 
Z(y.4 — y.3)? = 0.0000082542 Z(y3t — y3)? = 0.000001113 
Z(y2! — y.2)? = 0.0003188165 (yt — ys!)? = 0.003619254 


y2* and ys* were found by means of the method of reference 1. 


7 Burgess, C. P., The Frequencies of Cantilever Wings in Beam 
and Torsional Vibrations, N.A.C.A. T.N. No. 746, 1940. 

8 Sokolnikoff, I. S., and Sokolnikoff, E. S., Higher Mathematics 
for Engineers and Physicists, 2nd Ed., p. 537; McGraw-Hill Book 
Company, Inc., New York, 1941. 

® Biezeno, C. B., and Grammel, R., Technische Dynamic, Chap. 
X, p. 16; Julius Springer, Berlin, 1939. 


The Effect of Radius of Curvature and Preliminary Artificial Eccentricities 
on Buckling Loads of Curved Thin Aluminum-Alloy Sheets for 
Monocoque Constructions 


(Continued from page 596) 


these special tests. On the contrary, it was found that 
the panels are not affected at all by these artificial 
eccentricities, and no sensitivity or weakness could be 
detected as a result of this sort of eccentricity. 


RESUME OF RESULTS 


A series of tests on panels of various radii of curvature 
showed that these panels, which had been rolled to a 
smaller radius and tested with the usual radius of 
24 in., had strikingly higher resistant qualities than 
panels of the usual shape. Between 8- and 14-in. 
radius, no buckling took place; the stress-strain dia- 
gram has a high elastic range and the collapse loads 
are distinctly superior to those of other panels with 
greater radii. These results have been confirmed by 
tests on panels made of pure aluminum. 


Furthermore, it was possible to show that initial 
eccentricities of the panel, generally supposed to be the 
main cause of a low buckling load, are not at all serious 
if they remain localized in the sheets. Indentations 
with 10-mm. and 25-mm. steel balls did not affect the 
strength qualities of the panels against compression 
loads. Even holes produced by shooting 0.22-in. 
bullets through the panels did not decrease their normal 
buckling loads or their ultimate loads. 
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1 Welter, Georges, Curved Aluminum-Alloy Sheets in Compres- 
ston for Monocoque Constructions, Journal of the Aeronautical 
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* Welter, Georges, Influence of Different Factors on Buckling 
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No. 4, pp. 204-213, April, 1946. 
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Vibrations of a Helicopter on the Ground 


GABRIEL HORVAY* 
McDonnell Aircraft Corporation 


ABSTRACT 


The present paper serves as an introduction to ground vibra- 
tion analysis. Considerable insight may be gained of the 
phenomenon of ground vibrations, without actually going into 
the analysis of these somewhat complicated self-excited oscilla- 
tions, by limiting one’s attention to the much simpler problem— 
that of forced vibrations of the helicopter on the ground with the 
blade hinges locked. This preliminary study permits the de- 
termination of the natural modes of the helicopter, the ratios of 
the various amplitudes at resonance, and the magnitudes of the 
oscillations and of the blade unbalance forces at which shock- 
strut action sets in. The connection between the present study 
and ground vibration analysis is also taken up briefly, and a nu- 
merical example is given. Finally, the spring and damping con- 
stants of a tire-shock-strut combination are determined, and a 
discussion is given of the various constants of the shock absorber 
which play a role in ground vibrations. 


INTRODUCTION 


| pasevcesaiel CREATED by the interaction of a heli- 
copter’s rocking on the tires with the oscillations of 
the rotor blades about the lag hinges is known as ground 
vibration. In earlier construction of rotorcraft there 
were no lag hinges and no ground vibrations. How- 
ever, rotor blades were subjected to large alternat- 
ing chordwise bending moments, and breakages oc- 
curred. 

To eliminate the large variations in the bending mo- 
ments at the blade roots and to provide smoother per- 
formance, lag hinges came into use. The hinges pro- 
vided much improvement, but the danger of ground 
instability also entered the picture. Many rotorcraft 
were destroyed on the ground because of the new kind 
of instability. The trouble was remedied by the intro- 
duction of lag motion dampers attached to the blades. 
The dampers eliminated ground vibrations, but they 
created fairly large bending moments at the blade roots, 
and the structural strength of the blades has again be- 
come a matter of serious concern. It is therefore im- 
portant to know what is the least amount of blade 
damping that is adequate to prevent instability. Of 
course, this determination may be carried out experi- 
mentally and with safety for any given helicopter, as 
described by Deutsch in his recent paper.! The de- 
termination of the general effect of the various param- 
eters of the helicopter—such as height of helicopter 
c.g., height of rotor, width of wheel base, eccentricity of 
lag hinge, coupling between lateral motion and rolling, 


Presented at the Rotating Wing Aircraft Session, Fourteenth 
Annual Meeting, 1.A.S., New York, January 29-31, 1946. 
* Research Engineer. 


etc.—upon the ground vibration characteristics is a_ 
study that requires a theoretical approach. 

Coleman? was the first to compute the amount of 
blade damping required for a given helicopter on the 
basis of the following simplifying assumptions: (1) 
The rotor has three blades; (2a) H = height of heli- 
copter c.g. above axle = 0; (2b) A = height of heli- 
copter rotor above axle = 0; (3) flapping hinges are 
locked; (4) blades are rigid; (5a, b, c) pitching, roll- 
ing, yawing of the helicopter on the tires can be dis- 
regarded; (5d) vertical bouncing of the helicopter on 
the tires can be disregarded; (5e) gyroscopic effects can 
be disregarded; and (6) elasticity and damping exist 
only between the tires and ground and at the lag 
hinges. 

Conditions 5a, b, d, and e really follow from conditions 
2aandb. Since axle center, c.g., and rotor center coin- 
cide, no excitation of the pitching, rolling, and vertical 
modes takes place; therefore these modes may be omit- 
ted from consideration. 

Coleman found the important result that, when the 
product of blade damping and lateral landing gear damp- 
ing exceeds a certain quantity, stability is obtained. 

The present paper is based on a study, recently under- 
taken at the McDonnell Aircraft Corporation under 
the sponsorship of the Bureau of Aeronautics, U.S 
Navy, which aims at a more accurate determination of 
the required blade damping than was available hereto- 
fore. Therestrictions 2, 3, and 5 of Coleman are dis- 
carded in that study, and assumption 6 is relaxed to in- 
clude pylon bending modes. 

Because of the extended length of the derivations, the 
above study cannot be presented within the scope of 
this paper. The most important result, the criterion 
for coupled lateral-rolling instability is given in the 
section “Comments on Ground Vibrations’ without 
proof (proofs to the statements will be supplied in refer- 
ence 3). Elsewhere the present investigation deals with 
the simplified problem where the lag and flapping hinges 
are locked, the forces of self-excitation are replaced by 
impressed forces, and the helicopter is considered as 
rigid except for the elasticity and the damping provided 
by the landing gears. 


DESCRIPTION OF THE HELICOPTER 


In what follows it will be assumed that the helicopter 
is a twin-rotored side-by-side type (McDonnell-type) 
t The writer wishes to express his appreciation to his colleague, 


A. Gail, and to Lt. Comdr. L. Arnold, U.S.N.R., Bureau of Aero- 
nautics, for helpful discussions. 
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helicopter, with right rotor center at Op, left rotor 
center at O, at a distance 2Z from Og. L is the length 
of the pylonarm. The rotor midpoint (i.e., the midpoint 
of O,Og) at O is assumed to be directly above C, the 
c.g. of the helicopter. 

For a single rotor helicopter Og coincides with O, and 
there is no left rotor. Other changes necessary to 
adapt the present analysis to single rotor helicopters 
will be noted at the appropriate places. 

Fig. la gives a schematic front view, Fig. 1b a sche- 
matic top view, and Fig. lc a schematic side view of the 
helicopter. For the sake of simplicity, the rotor 
blades and the points Og, Oz are omitted from the dia- 
gram, and only the rotor midpoint O is shown. A-A, 
of length 2/,,, is the axle of the main landing gear, F its 
center. C is at height H, O at height 4 above the axle. 
The vertical axis OC pierces the horizontal plane 
through A-A at point S. F is at distance /, forward of 
S, the rear axle R is at distance /, aft of S. RSF can 
be called the longitudinal axis. 

The center of gravity of the helicopter can be con- 
sidered as the center of a coordinate system. Dis- 
placements with respect to C are denoted by u (for- 
ward), v (to starboard side), Z (downward). 

The rigid body rotations of pitching £, rolling 7, 
yawing y are defined as positive when O moves forward 
with respect to C (this is contrary to the usual defini- 
tion), when O moves to the starboard with respect to 
C, and when RSF rotates clockwise about C as seen 
from above, respectively. It will be always assumed 
that the displacements u, v, Z and the rotations &, n, ¥ 
are so small that-their higher powers can be néglected. 

Figs. la and 1b show the vertical axis OCS displaced 
sidewise by v into position 0,C,S, and rotated by 7 into 
position O2C,S:. The longitudinal axis RSF assumes 
the positions R,S,F; and R,S2F2 in the course of these 
displacements. (Note that certain points, like S,, O., 


etc., are indicated by small circles in Fig. 1b but are 


not labeled because of lack of space.) A yaw by angle 
¥ brings R.S2F, into the R3S;F3 position. The axle 
A-A assumes at the same time the A;—A; position. 

Fig. le shows OCS moved forward by u into position 
O:C,S,, downward by Z into position 0.C,S2, and rotated 
through angle £ into position 03C353. 

Let W be the weight, I = W/g the mass of the heli- 
copter (blades included). Let J; be the centroidal 
pitching moment of inertia, J, the centroidal rolling 
moment of inertia, J, the centroidal yawing moment of 
inertia of the helicopter. Let W’, M’, I,’, I,’, 1,’ 
denote the above quantities when the helicopter is 
stripped of its blades. Naturally, in this case the 
moments of inertias must be referred to point C’ (the 
c.g. of the helicopter when stripped of the blades), 
which is at height H’ above the axle. 

Let the rotor blades have the length R. This in- 
cludes the lag hinge link length a. Let dm denote a 
mass element of the blade at distance 7 from the lag 
hinge. Then 


mp = Seiragdm, m= S,25%dm (1a) 


are, respectively, the “‘total mass’’ and the ‘‘outboard 
mass” (i.e., the mass outboard of the lag hinge or, 
briefly, the ‘“‘mass’’) of the rotor blade. 


mpB = Sytr=,(a + r)dm, mgps* = 
Sn,=-,(a+7r)%dm (Ib) 
yR5* rdm, mp? = Jf, 85° r'dm (1c) 
define the c.g. distance and the radius of gyration of 
the total blade with respect to the rotor center and of 


the outboard blade with respect to the lag hinge. 
The important ratio (Coleman’s A) 


¢ = ab/p? (2) 


istobenoted. It may be called the eccentricity param- 
eter of the lag hinge. ¢ plays a fundamental role in 
ground vibrations. Another numerical constant, 


6 = pb?/2p? (3a) 


mb = 


which can be called the interference parameter, likewise 
plays a fundamental role in ground vibrations. Here 


u = 3m/M (3b) 


is the ratio of the outboard blade masses carried by a 
single rotor to the total helicopter mass (including the 
blades). (6 is Coleman’s A; for single rotor helicopters 
and Coleman’s A;/2 for twin rotor helicopters.) 

In place of the centroidal coordinate system u, v, Z 
it will be convenient to introduce a new coordinate 
system X, Y, Z, with origin at S, the intersection point 
of the vertical and longitudinal axes. The horizontal 
displacements in the two coordinate systems are related 
by 

X =u-— fet (4) 
Y=v— Hn 
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while the vertical displacement Z is the same with 
respect to both systems. The pitching, rolling, and 
yawing moments of inertia of the helicopter with respect 
to S are 


L=1,+HM,I,=I1,+ HM,1, (5) 
The numerical constants 
i, = 3mppp*/I,, i = 3M ppp*/T, (6) 


(the gyroscopic coupling parameters of rolling and pitch- 
ing) are of importance in the theory of the single rotor 
helicopter. 

Let Ky, Ky, Kz denote the elasticities of the main 
landing gear in the X, Y and Z directions. Kz is a 
combination of shock absorber elasticity A, (which 
acts vertically) and vertical tire elasticity K, Ay and 
Ky represent the longitudinal and lateral tire elastici- 
ties. Capital letters always refer to a pair of tires, or 
struts, or their combinations; the lower case letters 
ky, ky, Rz, ki, ks will be used to denote the constants of 
a single tire or strut, or a single tire-strut combination. 
Let the damping constants Cx, Cy, Cz C,, C,, Cx, Cy, Cz, 
¢, ¢; be defined similarly. The tail wheel exerts re- 
straint only in the vertical direction. « is the spring 
constant, y the damping constant of the combined rear 
tire and strut. 


NOTATIONS 


Numbers in parentheses refer to the equations in which the 
quantities are introduced. Symbols introduced earlier are not 


repeated. 

£, F = lateral, longitudinal exciting force, Eq. (8) 

mM = exciting moment for yawing 

w = rotor frequency, Eq. (9) 

Q = vibration frequency, Eq. (16) 

D = d/dt, Eq. (11) 

D(v;) = quantity defined in Eq. (29c) 

D = (with various subscripts) displacement of 
shock absorber piston 

d = asin d/dt 

d = quantity defined in Eq. (38c) 

d = depth of the center of oscillations, Eq. (44) 

t = time 

t = coefficient of coupling between lateral mo- 
tion and rolling, Eq. (18b) 

r = height ratio, Eq. (18a) 

vy, Vy, Yt, ¥y, Wy = Natural frequencies of the pure modes 
X, Y, &, », ¥, Eqs. (12) and (14b) (the 
rotational modes are referred to S as 
center) 

Ax, Ay, ete = extinction coefficients of the pure modes 


X, Y, etc., Eq. (12) 
Vz, 1, ¥s, \z, Mt, As = Quantities defined in Eq. (14) 


q = vertical frequency ratio, Eq. (22) 

a = (lateral) frequency ratio, Eq. (22) 

Ay% = percentage critical of the pure Y mode, Eq. 
(15) 

Ex, Ey, Ex, E, = quantities defined in Eq. (11) 

M%, Ve = quantities defined in Eq. (19) 

A =v, + ite = frequency determinant of the lateral-roll- 
ing motion, Eq. (19) 

Aund = frequency determinant of the undamped 


system, Eq. (20) 


"4, V2 = natural frequencies of the coupled lateral- 
rolling mode, Eq. (21) 

Wi, We = critical rotor frequencies, Eq. (26) 

Ai, Ae = extinction coefficients of the »;, ve modes, 
Eq. (38) 

ly, In, Uy, Uy = coefficients in Eq. (23) 

€1, €2, €3, €4, Ai, A2 = coefficients in Eq. (37) 

ay = damping constant of the blade damper 

ry = blade extinction coefficient, Eq. (27) 

Vy = blade natural frequency, Eq. (32) 

D(w, 2) = O = frequency equation of ground vibrations, 
Eq. (25) 

Me, = apparent extinction coefficient of landing 


geai, Eqs. (28) and (29) 
Ly, Ly, Uy, Uy, = coefficients in Eq. (29) 
= gyroscopic splitting of the ¥, mode, Eqs 


Vy—» Yat 
(35) and (36) 
Vi-, "i+, ¥2-, V+ = gyroscopic splitting of the »,, ve modes 
Votat., Hnetat, = static deflections, Eq. (41) 
V = veitical force, Eq. (45) 
Vi = V/2, W, = W/2 = quantities used in the section on 
“Constants of the Shock Absorber”’ 
V*, ¥*, Hn*, £* = limiting values at which shock-strut action 
sets in, Eqs. (46) and (47) 
¥ = coefficient of packing friction of shock ab- 
sorber, Eq. (61c) 
A = piston area, Eq. (53) 
P = period, Eq. (60) 
v = (with various subscripts) air volume in 
strut 
p = (with various subscripts) air pressure in 
strut 


TIRE MODES OF THE HELICOPTER 


Rotor Stationary 


A lateral displacement v of the helicopter c.g. and a 
roll by angle » about the c.g. produce certain tire reac- 
tions. A lateral landing gear spring force Ky(v — Hn) 
opposes the sidéwise displacement v — Hn of the axle 
A-A. The roll » of A-A is supported by the moment 
HKy(v — Hn) of the lateral spring force and is opposed 
by the moment /,Kz/,n of the vertical spring force. 
If £ denotes the lateral exciting force which acts at 
rotor height, then 


v:; Mi + Ky(v — Hn) = £ } (7) 
n: 1,7 — HKy(v — Hn) + Kalen = (h — ELS * 


are the equations of motion of the system. By Eqs. 
(4) and (5) the equations can also be written in the form 


¥: MY + KyY + MH7 = £& } 
n. HMY + 1,4 + Kgly*n = hf 


Selection of the centroidal coordinate system u, v, Z 
implies static coupling between lateral motion and roll. 
Selection of the coordinate system X, Y, Z with origin 
at S replaces the static coupling by dynamic coupling. 
This is a distinct advantage in the present study be- 
cause it permits separation of the lateral and vertical 
spring effects of the landing gears. 

A forward displacement u of the helicopter c.g. and 
a pitching by angle — about the c.g. produces the fore- 
and-aft tire reaction Ky(w — Hé) which opposes the 
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forward motion, the moment HKx(u — He) of the 
fore-and-aft tire reaction which aids, and the moment 
LKzé + Lxl,£ of the vertical landing gear reaction 
which opposes the pitching £ If & denotes the 
longitudinal exciting force that acts at rotor height, 
then 


X: MX + KxX + MHE = § 
t: HMX + 1,€ + (Kal? + «lZ)& = hs 


are the longitudinal-pitching equations of motion 
of the helicopter referred to the S-coordinate sys- 
tem. 

For simultaneous excitation by the forces £ and F 
the equations of motion, Eqs. (8a) and (8b), hold simul- 
taneously, and there is no intercoupling between the 
lateral and the forward motions. 


\ (8b) 


Rotor Revolving, Blade Hinges Locked 


Consider now a single three-bladed rotor with center 
at O, revolving at angular speed w in a clockwise sense 
(sense of MM, Fig. 1). The centrifugal forces of the 
three blades cancel out. The Coriolis forces 


+2w(a + r)dmi cos wt; —2w(a + r)dmé sin wt (9a) 


which act on the blade element dm at the azimuth wt 
as measured from the downwind position, the first 
pointing down, the second up, create a pitching mo- 
ment 


28(a + r)2dm cos wt(—7 cos wt + E sin wt) (9b) 
and a rolling moment 
Qu(a + r)*dm sin wt(—4 cos wt + E sin wt) (9c) 


Integrating over the blade, cf. Eq. (16), and summing 
for the three blades, one obtains —3mgp,’wy and 
+3mgpp*wt for the two impressed gyroscopic moments. 

Taking also landing gear damping into account, the 
equations of motion of the single rotor helicopter on ground, 
with rotor revolving but blade hinges locked now assume 
the form 


X: MEx(D)X + MD*HE = § (10a) 

t: HMD*X + I,E,(D)E + IpugoDn = hF (10b) 

Y: ME,(D)X + MD*Hy = & (10c) 

n: HMD*Y + 1,E,(D)n — I,u,0Dt = h& (10d) 


where y:, uw, were given in Eq. (6) 


D = d/dt \ 


Ex(D) = D? + xD + z°f (11) 


and Ey, E;, E, are obtained from Ey by changing the 
subscripts X to Y, &, and 7. 
The quantities 


Vy? = Kx/M, vy? = Ky/M 
ve? = (Kal? + 11;?)/Ip ¥,? = Kebw®/Ir 


1 om 


are the squared natural frequencies of the pure X, pure 
Y, and (keeping S fixed) pure £ and pure 7 modes. 
(Note that the centroidal rolling frequency, keeping 
C fixed, is : 


Vy, centroid = V (Kaho? + KyH?*)/I, (13) 





It depends on both the Z and the Y elasticities of the 
landing gear.) Similarly, 


hy = Cx/ M. »Ay = Cy/M 


A: = (Czl/ _ yl,?)/TIp, r, a A (12b) 


are the extinction coefficients (or briefly, the “‘damping’’) 
of the above modes. 

For twin rotor helicopters with two identical rotors 
the gyroscopic terms cancel out because of the reversal 
of the sign of w for the second rotor. The equations of 
motion split into two groups: the modes of X, & and 
the modes of Y, 7. Thus, the equations of motion of the 
twin rotor helicopter on ground, with rotors revolving but 
blade hinges locked, are obtained from Eqs. (10) by omit- 
ting the gyroscopic coupling terms and were given earlier, 
except for the damping terms; in Eqs. (8a) and (8b). 
For single rotor helicopters the gyroscopic interactions 
produce a coupling between ~ and n, and the modes X, 
éand Y, 7 are no longer independent. 

When the lag hinges are unlocked, blade oscillations 
in the lag plane take place. These oscillations interact 
with the X, &, Y, 7 vibrations of the helicopter. It is 
found that the equations that govern the helicopter 
vibrations are still Eqs. (10). (Omit the 4, u, terms 
for twin rotor helicopters.) But £ must now be iden- 
fied with the lateral and $ with the longitudinal whip- 
ping forces of the hinged blades. Evidently, the blade 
unbalance forces also create an unbalance yawing 


‘moment 9 and thus excite yawing vibrations y of the 


helicopter. Less directly, vertical oscillations Z of the 
helicopter are also excited. It will be shown in refer- 
ence 3 that the Z motion is usually unimportant, while 
the ¥ motion is unimportant only when /, and L are 
small. 

Considerable simplifications and symmetry are 
achieved in the theory by restricting the degrees of 
freedom to the X, &, Y, 7 motions. This plan is fol- 
lowed in the present paper. 

For the sake of completeness the following quanti- 
ties, used in subsequent sections, are also defined at this 
place. 


vz? = K;/M, v2 = K,/M, v7? = K,/M_ (14a) 


are the squared vertical natural frequencies of the heli- 
copter on the front landing gears, on the front 
shock absorbers alone, and on the front tires alone. 
The extinction coefficients dz, \,, 4, are defined simi- 
larly. 


yy? = (Kxly? + Kyl?)/I, (14b) 


is the squared yawing natural frequency of the heli- 
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copter when the point S (or C) is kept fixed. The ex- 
tinction coefficient d, is defined similarly. 
The notation 


(15) 


will be used to indicate that the Y mode has a 6 per cent 
critical damping. 


- 100) Aya = Ay 2vy = Cyvy 2Ky = 0.06 


COUPLING BETWEEN LATERAL OSCILLATIONS AND 
ROLLING 


Natural Frequencies of the Twin Rotor Helicopter 


Eqs. (10), with the gyroscopic interaction terms yz, 
u, omitted, govern the motion of the twin rotor heli- 
copter. The X, £ equations are of the same form as the 
Y, » equations. Thus, for an understanding of the 
response of the system it suffices to investigate the Y, 
n motion. 

For a harmonic excitation 


L(t) = Le (16a) 
the responses are also harmonic 
Y(t) = Ye™, Hn(t) = Hne™ (16b) 


and Eqs. (10c) and (10d), with u, omitted, reduce to 

Y: (vy? — Q? + 2Ddy) Y — 2°An = L/M (17a) 
Hn: —t9?Y + (v,? — 2? + i2d,)Hn = rt£/M_ (17b) 
Here 


r=h/H (18a) 


is the height ratio of rotor and c.g. (referred to axle as 
base) ; 


t = MH?/I, (18b) 


is the coefficient of coupling between lateral and rolling 
motion. 

The system is in a natural mode if the determinantal 
equation of Eqs. (17) vanishes: 


A= V} + 102 = 0 (19a) 


where 
y= 24(1 —)=— 27(vy? + vy + AyA,) ? vy,” (19b) 


v2 = OQ[A, (vy? — 2%) + Ay(v,? — 22)] (19) 


The roots of the damped equation A = 0 are given 
(accurately, to the extent that second order small 
terms in Ay and X, are negligible) by the roots », and v 
of the undamped equation 


Auna, = (vy? — 0?)(»,? — 2?) — 4 = (1 — 1) X 
(2? — v,?)(Q? — »~*) = 0 (20) 


One obtains 


ee |: +a? + Vil — a’)? + a vy? (21) 
2(1 — 2) 
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where subscript 1 refers to the lower root, subscript 2 
refers to the upper root, and 


a = v,/vy (22a) 


is the lateral frequency ratio, or briefly, the frequency 
ratio of the system. a is related to the vertical frequency 
ratio 


v,/v2 = ¢ = lw t/H = V%,/dz = Xqu/dzu (226) 
by 
a = gVK;/Ky (22c) 


The natural frequencies »; and v2 are independent of 
the height ratio 7. », and v2 are plotted in Fig. 2 for 
a = 2, ~/2,1, V/'/2, '/2, when ¢variesfromOtol. For 
very small coupling, ¢ ~ 0, the frequencies 
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nt = {1 — [tat/(1 — at} Jat ) 
vw? = {1 + [t/(1 — a@?)] boy? (when a < 1) 


(when a = 1) 


yy? = vy?/(1 + V2), wy? = vy?/(1 —Vvi) + (23) 


vy? = {1 — [t/(a? — 1)] yy?, 
vy? = {1 + [ta?/(a? — 1)] boy (when a > 1) } 


,are essentially the uncoupled frequencies vy and », of 
lateral motion and rolling. As the coupling increases, 
v, becomes somewhat lower and v2 much higher. In 
the limit of > © (complete coupling) 


v2—> © (24a) 


and 
1 > ary/V1 + a? (24b) 


is the only natural frequency of the system. 


Comments on Ground Vibrations 


Without going into derivations or proofs, some com- 
ments about the application of the present theory to 
ground vibrations seems to be in order at this point. 

Eqs. (10) govern the oscillations of the helicopter 
also in the case when the rotor blades are free to oscil- 
late about the lag hinges. But the exciting forces F 
and £, which represent the blade unbalance forces, 
are now created by the oscillations X, &, Y, n themselves. 
This fact can be expressed by a second set of equations 
between X, é, Y, » and §, £; the phenomenon is one 
of self-excitation. As a result, the frequency of vibra- 
tion 2, which is the same as the frequency of excitation 
by § and £, cannot be arbitrary; Q is a prescribed 
function of the rotor frequency w. This functional re- 
lation assumes the form of a frequency equation 


D(w, 2) = 0 (25) 


One finds that for rotor speeds near the so-called critical 
speeds 


w, = v,/(1— Ve) (j = 1,2) (26) 


where 1, v2 are the two natural frequencies of the 
coupled lateral-rolling motion, the oscillation frequency 
2 of an undamped helicopter has a negative imaginary 
part. This indicates instability. The magnitude of 
the imaginary part gives the rate at which the oscilla- 
tions build up; the real part gives the frequency of 
oscillations. This is very nearly v,, indicating that the 
gravest instability sets in with a frequency that is 
approximately equal to a natural frequency of the heli- 
copter. ; 
Instability can be eliminated by introducing sufficient 
damping into the system. One finds that if the blade 


extinction coefficient 

de = 6/mp? (27) 
(c is the viscous damping constant of the blade damper 
in units of lb. ft. per rad. per sec.) and the ‘‘apparent”’ 
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landing gear extinction coefficient \, have a sufficiently 
large product 


Node 2 Svy*(1 — vVe)/Ve (28) 
then stability is ensured. Here d,, the apparent ex- 
tinction coefficient of the landing gear, is a linear combina- 
tion of ry and d,: 

Ag = Lydy + L,X, (in the lower mode »;) (29a) 
Ag = Uydy + U,), (in the upper mode v2) (29b) 


where 

Ly = [a°(vy?/m,*) — 1]/D() 

L,= [(vy?/v4?) — 1] /D(n) (29c) 

D(n) = a? + tr? — (1 + tr? — 2Qir)n,?/vy? 

and Uy, U,, D(v2) are obtained from Ly, L,, D(r,) by 
replacing »; by » in the above expressions. The ec- 
centricity parameter e and the interference parameter 
6 were defined earlier in Eqs. (2) and (3). 

For a > 1, ¢ = O (rolling has a higher natural fre- 
quency than the lateral motion and is uncoupled) or 
fora = © (Kz = ~&, pure lateral mode), Eq. (29a) re- 
duces to Ag = Ay, and (28) reduces to the well-known 


stability criterion of Coleman 
Mary S dry?/(1 — We)/We (30a) 
which can also be written in the form 
(¢o/mp*)(Cy/Ky) 2 61 — Ve)/VWe — (30b) 


(Note that Eq. (30) is a special case of Coleman’s more 
general formula when k, = blade centering spring stiff- 
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ness = 0. The effect of blade centering springs is 
omitted in the present analysis for convenience. Blade 
flexibility can be taken into account in an approximate 
way by determining an effective lag hinge position as 
the intersection of the unflexed blade and a mean secant 
of the flexed blade. This increases the value of «.) 
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For a = 0 (Ky = ©, pure rolling mode), formulas 


(28) and (29) reduce to 


MA, > tr%v,2(1 — Ve)/Ve (31a) 
which can also be written in the form 
> i=¥e (31b) 





(wis) )Mata) > 7 


From (31) it is evident that a large rolling moment 
of inertia I,, as provided by twin rotor helicopters with 
long pylon arms, is advantageous in ground vibrations. 
However it must be remembered that in such a case the 
excitation of the yawing mode y is no longer negligible. 
When the distance /; of the c.g. behind the axle center 
is sufficiently short and the brakes are off (Ky = 0), 
then the yawing frequency vy, is sufficiently low so that 
its effect on ground vibrations can be disregarded. 

The curves of Ly and L, are plotted in Figs. 3a—3f vs. 
tand rfor a = 1/+/10, '/2,- V'/z, 1, V2,2. The curves 
of Uy and U, are plotted in Figs. 4a—4d vs. ¢ and r for 
a = 0, 1/2, 1, 2. It is found that a large value of the 
height ratio r is detrimental in the lower and usually 
favorable in the upper mode. Since, however, the 
lowest mode is almost always the most critical mode, 
large r is undesirable in ground vibrations. For a large 
frequency ratio a most of the excitation goes (at »;) into 
the slightly damped lateral mode vy; for small a and 
small coupling ¢ most of the excitation goes into the 
rolling mode yv,, which is greatly damped by shock ab- 
sorber action. Small a and small ¢ are therefore de- 
sirable conditions in ground vibrations. They are best 
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achieved by twin rotor helicopters with long pylon 
arms. Naturally, the interference 6 should be as small, 
the lag hinge eccentricity « as large, and the natural 
frequency »; as small as permitted by other design con- 
siderations. 
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The preceding considerations apply strictly to twin 
rotor helicopters with sufficiently low yawing frequency 
vy, sufficiently small 6, and sufficiently large e, so that 
the required Aygy, A,%, Noy, as determined from Eq. (28) 
do not become excessive (i.e., they stay below, say, 
35 per cent). Here A» must be referred to the blade 
natural frequency 


n= w/e (32) 


and thus has different values in the two modes », and 
vo. If dy, Ay, A» are excessive, then the criterion, Eq. 
(28), breaks down, and the cumbersome frequency 
equation, Eq. (25), must be used to establish the con- 
ditions for stability. This will be done in reference 3. 
However, for most helicopters Eq. (28) is adequate. 

For isotropic single rotor helicopters the same sta- 
bility criterion, Eq. (28), applies provided the expres- 
sions of v,, ¥1, v2 are corrected for the gyroscopic splitting. 
This calculation is explained in the next section, and is 
usually somewhat cumbersome. For most purposes 
the gyroscopic splitting can be neglected; then Eq. 
(28) and Figs. 3 and 4 apply directly to the isotropic 
single rotor helicopter. 


* Natural Frequencies of the Isotropic Single Rotor Heli- 


copter 


The methods of the preceding sections apply im- 
mediately to the cases Kz = © and vy/vy = O or © of 
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the single rotor helicopter. In the present section 
the equally important special case of isotropic sym- 
metry 

er oy eS Vo» Ax = hy, A: = Ay Mg = My (34a) 


is considered. The equations of motion, Eqs. (10), 
now become 


X + iY: [D® + dAyD + vy] (X = tY) + 
D°H in) = (F + 1L)/M 
“* in) = (F = iL)/ (34b) 
H(t = in): tD*(X += 7Y) + [(D? +A,D + 


v,? = pwiD|H(— + in) = tr(F + i£)/M 


and can be brought to the form of Eqs. (17), provided 
one writes 


Vat? V_? a+ we (35) 


v,—-? = v,? — pw 


in place of v,?._ Thus, because of gyroscopic interaction 
between rolling and pitching, v, splits into two frequen- 
cies, v,— and v,+, which are quite close to each other 
since yw, is small. Correspondingly, the natural fre- 
quencies »; and v2, obtained from Eq. (20), also split 
into the two closely adjoining frequencies »;-, + and 
Ve—, Vz+, respectively. Unfortunately, v,— and v,+ de- 
pend on both the rotor frequency w and on the vibration 
frequency 2. For a given rotor speed w the frequency 
equation, Eq. (25), can be solved for Q, and then the 
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values in Eq. (35) can be determined. This procedure 
is indicated in reference 3. For stability calculations 
by the criterion, Eq. (28), it is sufficient to know the 


values of v,+?, v,—? assumed at the critical speeds w; and . 


wo. These are 


2 
Vat 


Vy? + Byn?/(1 — ey 
ry — myn2/(1 — Ve) san 


at the speed w:. Replacing », by ™, one obtains the 
expressions of v,+? and v,—? at the speed w». 

The frequencies v,+ and 4+ can now be determined 
in the following way: Solve Eq. (20) for » using the 
approximation v,+ = v,, and then, from Eq. (36), de- 
termine y,+. Replaced into Eq. (20), this yields a 
value for 1+. One such iteration is usually sufficient 
since yw, is small. The determination of the pairs 
{v,—, Vy— } ; {r+ Vo+ \ and i; Vs— } is similar. 


2 
Va- 


Effective Damping of the Normal Modes 


It is well known that the equation 


Q4 + @,23 + eQ? + 2 + ae = 0 
(a? ~ e32 ~ 0) (87a) 


factors into 
(Q? — A,Q — »?)(Q? — AQ — m7) = 0 (37b) 


where v;? and v2? are the roots of 


Q4 + @2? + eg = 0 (37c) 
and 
2 2 
Ain és + ei" ii~ es + e1r2" (37d) 
ve? — vy? yy? — vp? 


Correspondingly, Eq. (19a) factors into 

A = + tv, = (1 — £)(Q? — 1A,Q — v7) X 

(Q? — 4deQ — v*) (38a) 
where vr,’ and 7” are given by Eq. (21), and 
M1 = lyAy + La de = Uydy + Uyry (38b) 
with 
by = [(1 — é)(Qs? — v7?) — 1 — a? + 

2(1 — t)a*|/d 

1, = (1 — (4? — nt) — 1 — a + 


2(1 — #)]/d 
uy = [(1 — 2)(m? — mn?) +1+ a? — > —(38c) 
2(1 — t)a*|/d 
u, = [(1 — £)(m? — m4?) + 1+ a? — 


2(1 — 4))/d 








d = 2(1 — £)?(* — »?) 


\ and dz are the extinction coefficients of the free vibrations 
of the lower and upper modes v, and v2. The coefficients 
ly, 1, uy, U, are plotted vs. ¢ in Fig. 5 for a = 3/2, 1, 
and 2. 


For ¢ = 1 the frequency equation degenerates into a 
cubic and factors into 


A = —i(A, + Ay)(Q — a2) (Q? — AQ — w?) (39a) 
where 


Ae = (vy? + v,7)/(Ay + A,) 
nm? = a®vy?/(1 + a?) 
Ar = (a*Ay + A,)/(1 + @?)?, 


There is now only one oscillating mode, . This is as- 
sociated with the extinction coefficient, \;. The other 
mode represents an aperiodic motion with the very high 
extinction coefficient, ),. 

It is interesting to note that the coupled system a = 
t = 1, Ay = X, has a 50 per cent smaller effective damp- 
ing than has either the mode vy or the mode »,. From 
Eq. (39) one obtains 


A = Vay + Vad, = ‘dy (40a) 


(39b) 


Part of the decrease is compensated by the decrease in 

frequency; the percentage critical in mode » is 29 

per cent less than the percentage critical for either Ay 
Or v,: 

1/2 = Ay/V Bry = A, /V8r, (40b) 

The effective damping constants \i, \: of the present 


section are not to be confused with the apparent damp- 

ing constant A, of Eq. (29). 

Resonance Magnification. Amplitude Ratio. Center of 
Oscillations 


From Eqs. (17a) and (17b) it follows that, at any forced 
frequency Q, 








El ee Re poke oe 





y(Q 


is th 
to la 


Hyn(S 


is th 
Noti 
nan 


Y(y 


and 
Evic 
nanc 


of r. 


The 
1/, fi 


mod 


with 


belo 
and 
Hn/ 
large 
dam 
high 
tire 
itsel 
reas 


UNI 


data 
H/l, 
12,0 
per | 


med 


Tire 





toa 


9a) 


9b) 


as- 
her 
igh 





toe 


eS ee 





VIBRATIONS OF A HELICOPTER ON 





: ee eee 
( ) 1 v1 + tv JM 


{ } vy” V stat. 


is the vibration amplitude of the c.g. displacement due 
to lateral motion, and 


(41a) 





ra) = {= LN = OUI + ar) 8 
V1 + Ww M 
{ J v92Fnstat. (41b) 


is the amplitude of the c.g. displacement due to roll. 
Noting that Ay’? ~ A,” ~ Oand that 1, = 0 at the reso- 
nance frequency 1, 


Y(1)/YVstat. ~ vy* [tr — 1)n? + »,7] /tv2(r1) 


(resonance) (42a) 


and a similar formula, Eq. (42b), holds for Hn(v1)/Hstat- 
Evidently, a large height ratio r causes a large reso- 
nance magnification Y()/ Ystat.- 

The resonance amplitude ratio Hn/Y is independent 


of r. Placing £ = 0 in Eq. (17a), 
Hn/Y = (vy?/v;?) — 1 (resonance) (43) 
The ratio Hn/Y is plotted vs.¢ for a = 2, ~/2, 1, ~W/*/s, 


1/, for the lower mode » in Fig. 6a and for the upper 
mode vy; in Fig. 6b. 

The center of oscillations, the intersection of O2C,S, 
with OCS, Fig. la, is at depth 


d = Y/n = H/(Hn/Y) (44) 


below the axle in the mode », when /n/Y is positive, 
and is at height d above the axle in the mode ». when 
Hn/Y is negative. It is desirable to have small Y and 
large },.» vibrations. The former limits the slightly 
damped lateral tire amplitude; the latter makes the 
highly damped vertical amplitude of shock absorbent 
tire combination more pronounced. The angle 7 
itself must remain small because of static stability 
reasons. 


EXAMPLE. THE SIZE OF THE LATERAL TIRE 
AMPLITUDE, AND THE MAGNITUDE OF THE ROTOR 
UNBALANCE Force AT WHICH SHOCK-STRUT ACTION 
SETS IN. SELECTION OF BLADE DAMPERS 


Consider a twin rotor helicopter with the following 
data: W = 6,950 lbs., J, = 31,800 slugs ft.*, H = 7 ft., 
H/l, = 0.8, h/H = 16; ¢€ = 0.1,6 = 0.01; K, = 
12,000 Ibs. per in., Ky = 6,000 Ibs. per in., Aig, = 5 
per cent, Aye, = 5 per cent (tires only). 

The following derived constants are obtained im- 
mediately : 


Tire Action Only 


Mass = M = W/g = 216 slugs = 18 lbs. sec.? 
per in. 
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Rolling moment of inertia, J, = 31,800 + 
10,600 = 42,400 slugs ft.” (5) 
Coefficient of coupling, ¢ = 10,600/42,400 = 


25 per cent (18b) 
Vertical vibration frequency, vy, = 25.8 rad. 

per sec. (14a) 
Lateral vibration frequency, vy = 18.24 rad. 

per sec. (12a) 
Vertical frequency ratio, g = 0.625 (22b) 
Lateral frequency ratio, a = 0.883 (22c) 
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Rolling frequency, v, = 16.12 rad. per sec. (22a) 
Percentage critical, pure lateral mode, Ayg, = 

5 per cent; Ay = 1.824/sec. (15) 
Percentage critical, pure rolling mode, \,q = 


3.1 per cent; A, = 1.00/sec. (22b) 


By Fig. 2 (interpolating for a = 0.883 at ¢ = 25 per 
cent) 


vy, = 0.763rvy = 13.9 rad. per sec. 
1.34vy = 24.4 rad. per sec. 


v2 
are the natural frequencies of the coupled system, 
and 

w, = 194 r.p.m., w. = 340 r.p.m. (26) 


are the critical speeds. 
By Figs. 6a and 6b 


(Hn/Y),, = 0.72, (Hn/Y),, = —0.44 


are the amplitude ratios in the natural modes. At the 
first resonance the vertical deflection of the tire is mag- 
nified by the factor 


hon(v1) /hoMstat. = 19 (42b) 


The apparent extinction coefficient of the landing 
gear in the lower mode has the value 


A, = 0.365 X 1.824 + 0.775 XK 1.00 = 1.44/sec. (29) 
(This is somewhat worse than Ay.) One therefore finds 
that a blade damping of ‘ 


¢,/mp? = 7.2/1.44 = 5.0/sec. (27, 28) 


is necessary to ensure stability in the absence of shock- 
strut action. 


The Least Excitation Which Ensures Shock-Strut Action 


Assume now that the coefficient of packing friction 
of the strut is f = 0.40. Then the vertical force per 
tire set must reach the value 


V* = fW = 2,780 Ibs. (45) 


and the vertical tire deflection must reach the 
value 


Lon* = fW/K, = 0.232 in. (46a) 


before shock-strut action can set in. In the lower 
mode the corresponding lateral amplitude is 


Y* = (An/Y)—'(A/ly)hon* = 0.26 in.  (46b) 


and the lateral unbalance force that produces this 
deflection is 


£* = [(Kahe/h)lwnsear./bon(rs)]on* = 116 Ibs. (47) 


(Note that if the frequency ratio a is large (say, a > 
4/2), then the lateral amplitude may become so large 
(say Y > 2 in.) that the helicopter overturns before 
shock-strut action sets in.) 


1946 
The Required Blade Damping When the Shock Struts 
rate 


For further analysis the spring and damping constants 
of the shock struts must be specified. Assume these to 
be 


K,; = 6,000 Ibs. per in., A,¢, = 60 per cent 
per pair, then for combined tire-shock-strut action 


18.24 rad. per sec., a = 0.525, 1 = 9.12 rad. 
per sec. 


Vy 


Ayy = 5 per cent, A,% = 19 per cent 
and therefore 
Ae = 0.14 XK 1.824 + 4.25 X 3.64 = 15.7/sec. 


is the apparent landing gear extinction coefficient. 
This yields 


Nh = 7.2/15.7 = 0.46 sec. 


for the extinction coefficient of the blade dampers. 
This is ten times smaller than the previously obtained 
value 5.0/sec. for the case of no shock-strut action. 

Evidently an adequate knowledge of shock absorber 
characteristics is of primary importance in ground vi- 
bration analysis. 


Comments on the Values of the Spring and Damping Con- 

stants of Helicopter Tires 

R. E. Crowe and the writer have tested five pairs of 
high-pressure helicopter tires for spring and damping 
constants. The axle of the wheels supported a tray 
loaded with weights varying from 4,000 to 11,000 Ibs. and 
measured the natural frequency and the decrement of 
the free oscillations. They found that the vertical 
spring constant K, of the tires is about one or two times 
the static spring constant (as obtained from static de- 
flection tests), the lateral spring constant Ky is about 
one-third to one-half of K,, and the fore-and-aft spring 
constant Ky (only one set of tires was tested for Ky) 
is about the same as K, or somewhat higher. The tire 
damping constants Cx, Cy, Cz range from 4 to 7 per 
cent critical for the various tire sets. Cy seems to be 
somewhat larger than Cz. The spring constants in- 
crease with increasing tire pressure and decrease with 
increasing amplitude of vibration (the latter variation 
is pronounced). The damping constants decrease 
somewhat with increasing tire pressure. Highest K, 
(desirable, it makes K,/K, large), K,/Ky (undesirable, 
it makes a large), and Cx, Cy, Cz (very desirable) were 
found for channel tread tires. The above numerical 
values refer to tires bouncing on steel. Using wood as 
foundation little change was noted, but on dirt K, (and 
this alone) was invariably higher. 


COMBINED TIRE-SHOCK-STRUT ACTION 


When the vertical vibratory force V cos Qt acting on 


the landing gear is sufficiently large, then the condi- 
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tion V > V* = fW necessary for shock-strut action, 
Eqs. (45) and (61c), is satisfied, and the vertical spring 
and damping ‘“‘constants’” Kz and Cz of the landing 
gear are expressions in the strut constants K;,, C,, in 
the tire constants K,, C,; and in the frequency 2. The 
purpose of this section is to establish these relations. 

Roughly speaking, oil and air in the strut are subject 
to the same displacement; the same holds true for the 
spring action and the damping action of the tires. The 
two systems are excited by the same vertical force. 
Correspondingly, by the well-known electrical analogy 
method, it can be said that the strut impedance K, + 
i2C, is connected in parallel with the tire impedance 
K,+ i2C;. Jointly they yield a resultant landing gear 
impedance of 


Kz + i0C, = [1/(K, + i8C,) + 1/(K, + 72C,)]— (48) 
Thus it is found that 
al K.K (Ki t+ K,) + (K,C? + K,C,?)Q? 
(K, + K,)? + (C, + C,)?2? 
— KPC. + KPC, + CCC, + C,)2? 
(K, + Ky)? + (C, + C,)?2? 
are the spring and damping constants of the landing 


gear. 
If the struts jam, then K, = © and, therefore, Kz = 
K,, Cz = C;. For very weakly damped struts 


K; and K,? > > C,2? >> C,*0? 





(49a) 


(49b) 





Cz 


are true. They lead to the familiar relations 
Kz ~ [(1/K,) + (1/K;)]~ (50a) 
Cs stead C.K ?/(K; oa K,)* (50b) 


Thus, spring action of tire and very weakly damped 
strut is that of two springs K, and K, connected in 
tandem; damping action is that of damper C,, reduced 
by the factor K,?/(K, + K,)*. In order to produce 
large strut damping it is necessary to have a large value 
for the ratio K,/Ks. 

In Figs. 7a and 7b the ratios Kz/K, and C,Q/2K, 
are plotted. The curves in Fig. 7a show that the land- 
ing gear spring constant varies between [(1/K,) + 
(1/K,)] — for very small strut damping and K, for very 
large strut damping. (C,2 can always be considered 
as small.) 

For C,2. = 0 the numerator C,0 of the ratio C,Q/2Kz 
reaches a peak 


{C22} maz, ims 1/o(K /Ks)((1/K,) + (1/K,)] Tan (51a) 
at 
CR = Ky + K, (51b) 


and then sinks to C;,Q = Oas CQ — ~. For CQ = 
small ~ 0, the curve C,0 is similar, except that now 
C2 CQas CQ — @. 

For 2 = vz the ratio C,2/2Kz represents the frac- 
tion of critical damping that is available in the vertical 
mode. C,2/2Kz reaches its greatest values in the range 
C,2/K; = 1 to 3. For conventional landing gears this 
rarely exceeds 35 per cent. Fig. 7b clearly indicates 
the importance of having a large K,/K, ratio, or else, 
of having a very large value for the vertical tire damp- 
ing. If CQ becomes excessive, then the shock struts 
lock and the available landing gear damping is that of 
the tires alone (asymptotes, Fig. 7b). 


Example 


Determine the combined tire-strut spring constant 
Kz and damping constant Cz of the helicopter of the 


> eee ° 
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previous section, at the lower natural frequency 1. 
Given: K, = 6,000 Ibs. per in., As, = 60 per cent. 
(Note that v, = 18.2 rad. per sec., and C,v,/K, = 1.2.) 
As first approximation take 2 = 0. Then from Fig. 
a, Kz = 0.666K, = 4,000 Ibs. per in. Therefore 
= 0.51, and by Fig. 2a, by interpolation, »; = 0.49vy 
8.95 rad. per sec. (first approximation). , 
In second approximation take 2 = 1, approx.- 


“J 


y 8 


Then 
C,0/K, = 2(ds9,/100)(r1/v,) = 1.2 X 0.49 = 0.59 (52a) 


C/K, = 2(A1m,/100)(K i/Ks)(v1/v) = 0.10 , 4 


2 X 8.95/25.8 = 0.07 (52b) 


From Fig. 7a the corresponding vertical spring constant 
is Kz = 0.71K, = 4,240 lbs. per in. Therefore vz = 
15.32 rad. per sec., v, = 9.59 rad. per sec., a = 0.525, 
and, by Fig. 2a, 1 = 0.50vy = 9.12 rad. per sec. (second 
approximation). From Fig. 7b one now finds, for the 
curve K,/K, = 2, C/K, = 0.07, at abscissa C,»/K, 
= 1.2 X 0.50 = 0.60 that Cz,/2Kz ~ 18 per cent 
(referred to the »; mode). Consequently, 


Ay% = 100 (Czn/2Kz)(»,/n) = 18 X 


9.59/9.12 = 19 per cent. (52c) 


is the percentage critical of the tire strut combination in 
the rolling mode, using the Kz and the Cz values of the 
vy, mode. This is the quantity that enters \, of Eq. 
(28) (see the example of the preceding section). 


CONSTANTS OF THE SHOCK ABSORBER 


The shock absorber serves a dual role in ground vi- 
brations. It acts both as a spring and as a damper. 
Thus, both the spring constant and the damping con- 
stant of the shock absorber are of great interest. There 
is, in addition, a third constant, no less important than 
the other two, which must be considered. This is the 
coefficient of packing friction. Evidently, shock ab- 
sorber action can start only when packing friction has 
already been overcome. 


Air-Spring Action 

Let D, denote the static displacement of the shock ab- 
sorber piston from fully extended position (Dy) = 0) 
to the static position assumed under load W,, let D, 
denote the total displacement (stroke) of the piston from 
Dy to fully compressed position, and let D, denote the 
virtual displacement of the piston—i.e., the displacement 


D, = D, + Dray, 


where D,., is the displacement that would be neces- 
sary to eliminate also the small residual air pocket that 
remains in fully compressed condition. Let po, ps, p:, 
bd, = © denote the air pressures of the strut, and 
Vo, Vs, Vt, Vp = O the air volumes when the piston is in 
positions D) = 0, D,, D;, D,. In the static position 


p; = W,/A,v = A(D, — D,) (53) 


1946 


where A is the piston area. In any other position, 
reached in the course of a vibration, the adiabatic rela- 
tion holds 


b = p,(v,/v)*4, v = A(D, — D) (54) 


Then the spring constant of the strut in position D is 
given by 


“n= Sta 2 (#)(*) i: 14a? (62) x 
dD dv /\dD Vs 


2-4 2 
(=n) os 
v Vv 


k, = k(D,) = 1.4W,/(D, — Ds) (56) 


where 


is the dynamic spring constant of the strut in the static 
position (as contrasted with the static spring constant 
in the static position 


kg, stat. — W,/(D, ba? D,) (57) 
which is of no interest). 


Example 


Consider a helicopter of weight W = 2W,; = 7,000 
Ibs., with two shock struts of 3-in. piston diameter (A 
= 7:07 sq.in.), 12-in. stroke (D, = 12 in.), 1-in. taxiing 
stroke (D, = 11 in.), and a residual airpocket of 2.65 
cu.in. (D, = 12.375 in.). Then 


k, = 1.4 X 3,500/1.375 = 3,560 Ibs. per in., 
bs = 3,500/7.07 = 495 Ibs. per sq.in. 


and 


Po = pys/¥o = 495 X 1.375/12.375 = 
55 Ibs. per sq.in. (58) 


Here fp» is the air pressure in the completely extended 
position of the strut. 

The stiffness of the air spring is k, for vibrations of 
small amplitudes. When the amplitudes are large, 
then the motion becomes asymmetric about the posi- 
tion D,, because of the nonlinearity of the air spring. 
If the strut is just bottomed (Dnar. = D; = 12 in.) in 
the course of a vibration about’ D,, then the other ex- 
treme position is not Dyn, = 10 in., but rather Dmin. 
= 8.5 in. This can be seen from Sussholz’s‘* Fig. 3b 
or Fig. 4. To a downward amplitude of 1 in. = D, — 
D, = 0.73 (D, — D,) from position D, there corresponds 
a rebounce to an upward amplitude of 1.85 (D, — D,) 
= 2.5 in. The increase in amplitude is associated 
with an increase in the period of vibrations. Instead 
of the value v, = +1/ W,/k,, which holds for infinitesimal 
oscillations, the frequency is now, by Sussholz’s Fig. 5, 
0.88 v,. Correspondingly, the effective spring constant 
of the air spring has the value 





keg = 0.88? XK Rk, = 2,750 lbs. per in. 
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for oscillations between Dn = 8.5 in. and Dngz = 
12 in. Evidently, k.g decreases as the amplitude of 


vibration increases. The reason for this is that the pis- 


ton spends a good deal more time in the region of small 
stiffness than in the region of large stiffness. (Note 
that the quantity 


an al tied tt kdD/(Dymaz. 


5,950 Ibs. per in. in the present case, is of no physical 
significance.) 

It is interesting to note that if the load on the strut 
is cut in half, then the period of vibrations P is increased 


by a factor of »/2, because 
P.’ se 1/2Wi/Rsn a \/2s0/Vs0 et 4 = 2 (60a) 
P,? W,/Rso Pen/Ven m4 Du" 
Here the subscripts and 0 refer to new and old. Like- 
wise, if the air pressure in the extended position is doubled, 
then the period of vibrations is increased by a factor of 
1/2, because again 
P,? ois W1/Rsn BA Ps0/Vs0 —_. - 
P,? Wi/Rso = Ps0/Vsn F 


These facts are not always appreciated. 


a Din.) (59) 


(60b) 





Damper Action 


A well-designed shock absorber should be approxi- 
mately critically damped. This means that for 
vibrations at the natural frequency », = Vk,g/W, of 
the shock strut the damping force should be about twice 
the spring force, c,v,/k, ~ 2. Since the spring constant 
is a rapidly increasing function of the piston displace- 
ment D, the damping constant must be an equally 
rapidly increasing function of D in order to produce 
uniform damping throughout the stroke. A suitable 
choice of the taper of the metering pin permits an 
approximation to these conditions. 

The damping constant of a shock strut can be de- 
termined, for instance, from an oscillogram that shows 
the relative displacement of the piston and the wheel 
axle in the course of a drop test. This, however, yields 
only a certain average for c, because the motion is 
nonlinear. For small oscillations about the position 
D, c, is a constant and can be calculated from the prop- 
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erties of the liquid, the pressure, the speed of flow, and 
the nature of the orifice. 


Packing Friction 

Shock absorber action does not set in until pack- 
ing friction has been overcome. If a steady static 
force 


Vi => fiWi (61a) 
is just sufficient to overcome packing friction, then an 
alternating force V; cos 2 is evidently insufficient. An 
alternating force V;' cos Qt where V;’ > Vi, exceeds V; 
only in part of the cycle. The higher the frequency of 
vibrations, the greater must be the ratio 

Vi'/V = | + fe (61b) 
before strut action can set in. A further increase of 
the force, from V;’ to V,*, may be necessary to ensure 
an approximately sinusoidal*response. Denoting the 
factor of increase by 1 + fs, 


V* = (1+ ff) + fo)fims =fMm 


The factor f can be called the dynamic coefficient of pack- 
ing friction. Although no estimates for f are available, 
it is believed that at frequencies of abuut 4 cycles per 
sec., the value 


(61c) 


f~ 04 (62) 
used earlier in the text is not unreasonable. An ex- 
perimental determination of the coefficient f; has been 
carried out by MacLean.’ He found that under favor- 


able conditions f ~ 0.05 to 0.07. 
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